-25/8 - (109*x)/8 + nl*(-1/4 - 61/(36*x) - x/4 + (109*x^2)/36) + 
 Pi^2*(3 + (3*x^2)/2 - 3/(4*(1 + x))) + 
 (-75/4 + nl*(-3/2 - (13*x)/6) + (33*x)/4 - 33*x^2)*HPL[{0}, x] + 
 (-18 - 9/(-1 + x) + 9/x + 9*x - 9*x^2)*HPL[{2}, x] + 
 (18 + 9/x + 9*x + 9*x^2 - 9/(1 + x))*HPL[{-1, 0}, x] + 
 (nl*(-2/3 - (2*x)/3) - 9/(2*(-1 + x)) + 18*x - 9*x^2 + 9/(2*(1 + x)))*
  HPL[{0, 0}, x] + (-18 - 9/(-1 + x) + 9/x + 9*x - 9*x^2)*HPL[{1, 0}, x] + 
 (67/4 - (5*nl)/6 - (3*Pi^2)/4)*PlusDistribution[0, 1 - x] + 
 PlusDistribution[-1, 1 - x]*(6 - (2*nl)/3 + (27*Zeta[3])/4) + 
 ep*(-25/8 - (109*x)/8 + nl*(-1/4 - 61/(36*x) - x/4 + (109*x^2)/36) + 
   Pi^2*(3 + (3*x^2)/2 - 3/(4*(1 + x))) + 
   (-75/4 + nl*(-3/2 - (13*x)/6) + (33*x)/4 - 33*x^2)*HPL[{0}, x] + 
   (-18 - 9/(-1 + x) + 9/x + 9*x - 9*x^2)*HPL[{2}, x] + 
   (18 + 9/x + 9*x + 9*x^2 - 9/(1 + x))*HPL[{-1, 0}, x] + 
   (nl*(-2/3 - (2*x)/3) - 9/(2*(-1 + x)) + 18*x - 9*x^2 + 9/(2*(1 + x)))*
    HPL[{0, 0}, x] + (-18 - 9/(-1 + x) + 9/x + 9*x - 9*x^2)*HPL[{1, 0}, x] + 
   (67/4 - (5*nl)/6 - (3*Pi^2)/4)*PlusDistribution[0, 1 - x] + 
   PlusDistribution[-1, 1 - x]*(6 - (2*nl)/3 + (27*Zeta[3])/4)) + 
 ep^2*(-25/8 - (109*x)/8 + nl*(-1/4 - 61/(36*x) - x/4 + (109*x^2)/36) + 
   Pi^2*(3 + (3*x^2)/2 - 3/(4*(1 + x))) + 
   (-75/4 + nl*(-3/2 - (13*x)/6) + (33*x)/4 - 33*x^2)*HPL[{0}, x] + 
   (-18 - 9/(-1 + x) + 9/x + 9*x - 9*x^2)*HPL[{2}, x] + 
   (18 + 9/x + 9*x + 9*x^2 - 9/(1 + x))*HPL[{-1, 0}, x] + 
   (nl*(-2/3 - (2*x)/3) - 9/(2*(-1 + x)) + 18*x - 9*x^2 + 9/(2*(1 + x)))*
    HPL[{0, 0}, x] + (-18 - 9/(-1 + x) + 9/x + 9*x - 9*x^2)*HPL[{1, 0}, x] + 
   (67/4 - (5*nl)/6 - (3*Pi^2)/4)*PlusDistribution[0, 1 - x] + 
   PlusDistribution[-1, 1 - x]*(6 - (2*nl)/3 + (27*Zeta[3])/4)) + 
 ep^3*(-25/8 - (109*x)/8 + nl*(-1/4 - 61/(36*x) - x/4 + (109*x^2)/36) + 
   Pi^2*(3 + (3*x^2)/2 - 3/(4*(1 + x))) + 
   (-75/4 + nl*(-3/2 - (13*x)/6) + (33*x)/4 - 33*x^2)*HPL[{0}, x] + 
   (-18 - 9/(-1 + x) + 9/x + 9*x - 9*x^2)*HPL[{2}, x] + 
   (18 + 9/x + 9*x + 9*x^2 - 9/(1 + x))*HPL[{-1, 0}, x] + 
   (nl*(-2/3 - (2*x)/3) - 9/(2*(-1 + x)) + 18*x - 9*x^2 + 9/(2*(1 + x)))*
    HPL[{0, 0}, x] + (-18 - 9/(-1 + x) + 9/x + 9*x - 9*x^2)*HPL[{1, 0}, x] + 
   (67/4 - (5*nl)/6 - (3*Pi^2)/4)*PlusDistribution[0, 1 - x] + 
   PlusDistribution[-1, 1 - x]*(6 - (2*nl)/3 + (27*Zeta[3])/4))
