-8/27 + Pi^2*(-32/81 - 32/(81*x) - (8*x)/81) + 40/(27*x) - (32*x)/27 + 
 (-16/9 - (20*x)/27)*HPL[{0}, x] + (-64/27 + 32/(9*x) - (32*x)/27)*
  HPL[{1}, x] + (64/27 + 64/(27*x) + (16*x)/27)*HPL[{2}, x] + 
 (32/27 + (8*x)/27)*HPL[{0, 0}, x] + 
 ep*(-8/27 + Pi^2*(-32/81 - 32/(81*x) - (8*x)/81) + 40/(27*x) - (32*x)/27 + 
   (-16/9 - (20*x)/27)*HPL[{0}, x] + (-64/27 + 32/(9*x) - (32*x)/27)*
    HPL[{1}, x] + (64/27 + 64/(27*x) + (16*x)/27)*HPL[{2}, x] + 
   (32/27 + (8*x)/27)*HPL[{0, 0}, x]) + 
 ep^2*(-8/27 + Pi^2*(-32/81 - 32/(81*x) - (8*x)/81) + 40/(27*x) - (32*x)/27 + 
   (-16/9 - (20*x)/27)*HPL[{0}, x] + (-64/27 + 32/(9*x) - (32*x)/27)*
    HPL[{1}, x] + (64/27 + 64/(27*x) + (16*x)/27)*HPL[{2}, x] + 
   (32/27 + (8*x)/27)*HPL[{0, 0}, x]) + 
 ep^3*(-8/27 + Pi^2*(-32/81 - 32/(81*x) - (8*x)/81) + 40/(27*x) - (32*x)/27 + 
   (-16/9 - (20*x)/27)*HPL[{0}, x] + (-64/27 + 32/(9*x) - (32*x)/27)*
    HPL[{1}, x] + (64/27 + 64/(27*x) + (16*x)/27)*HPL[{2}, x] + 
   (32/27 + (8*x)/27)*HPL[{0, 0}, x])
