-6743/32 + 1643/(8*x) + (4693*x)/32 - (1643*x^2)/8 + 
 nl^2*(1/24 + 61/(216*x) + x/24 - (109*x^2)/216) + 
 Pi^2*(-63/8 - 33/(4*x) - (21*x)/8 - (165*x^2)/8 - 33/(16*(1 + x))) + 
 Pi^2*(18 + 9/x + 9*x + 9*x^2 - 9/(1 + x))*HPL[{-1}, x] + 
 (4143/16 + nl^2*(1/4 + (13*x)/36) + 201/(4*(-1 + x)) - (519*x)/16 + 
   (483*x^2)/2 + nl*(-8 + Pi^2*(-1/3 - x/3) - 5/(2*(-1 + x)) + 61/(12*x) + 
     (41*x)/12 + 3*x^2) + Pi^2*(-9/(2*(-1 + x)) + 9*x - 9*x^2 + 
     9/(2*(1 + x))))*HPL[{0}, x] + 
 (879/8 - 201/(4*x) - (75*x)/8 + (201*x^2)/4 + 
   nl*(-17/4 + 91/(12*x) + (13*x)/4 - (139*x^2)/12) + 
   Pi^2*(-9 + 9/(2*x) + (9*x)/2 - (9*x^2)/2))*HPL[{1}, x] + 
 (189/4 - 99/(4*(-1 + x)) - 99/(4*x) - (81*x)/2 + (495*x^2)/4 + 
   nl*(15/2 + 3/(2*(-1 + x)) - 3/(2*x) + 5*x + (3*x^2)/2))*HPL[{2}, x] + 
 (81/(2*(-1 + x)) - 27/x - 135*x + 54*x^2 - 27/(2*(1 + x)) + nl*(2 + 2*x))*
  HPL[{3}, x] + (-27/x - 81*x - 27*x^2 + 27/(1 + x))*HPL[{-2, 0}, x] + 
 (-63/2 - 297/(4*x) - (225*x)/4 - (297*x^2)/4 - 99/(4*(1 + x)) + 
   nl*(-3 - 3/(2*x) - (3*x)/2 - (3*x^2)/2 + 3/(2*(1 + x))))*HPL[{-1, 0}, x] + 
 (-54 - 27/x - 27*x - 27*x^2 + 27/(1 + x))*HPL[{-1, 2}, x] + 
 (261/4 + nl^2*(1/9 + x/9) - 99/(8*(-1 + x)) - 27*x + (495*x^2)/4 + 
   99/(8*(1 + x)) + nl*(-16/3 + 3/(4*(-1 + x)) + (55*x)/6 + (3*x^2)/2 - 
     3/(4*(1 + x))))*HPL[{0, 0}, x] + 
 (63/2 - 99/(4*(-1 + x)) - 297/(4*x) - (225*x)/4 + (297*x^2)/4 + 
   nl*(3 + 3/(2*(-1 + x)) - 3/(2*x) - (3*x)/2 + (3*x^2)/2))*HPL[{1, 0}, x] + 
 (108 + 54/(-1 + x) - 54/x - 54*x + 54*x^2)*HPL[{1, 2}, x] + 
 (54 + 54/(-1 + x) - 54/x - 108*x + 54*x^2)*HPL[{2, 0}, x] + 
 (108 + 54/(-1 + x) - 54/x - 54*x + 54*x^2)*HPL[{2, 1}, x] + 
 (108 + 54/x + 54*x + 54*x^2 - 54/(1 + x))*HPL[{-1, -1, 0}, x] + 
 (-81 - 81/(2*x) - (81*x)/2 - (81*x^2)/2 + 81/(2*(1 + x)))*
  HPL[{-1, 0, 0}, x] + (54 + 27/(2*(-1 + x)) - 108*x + 27*x^2 - 
   27/(2*(1 + x)) + nl*(-2 + 4*x))*HPL[{0, 0, 0}, x] + 
 (81 + 81/(2*(-1 + x)) - 81/(2*x) - (81*x)/2 + (81*x^2)/2)*
  HPL[{1, 0, 0}, x] + (108 + 54/(-1 + x) - 54/x - 54*x + 54*x^2)*
  HPL[{1, 1, 0}, x] + (201/2 - 5*nl - (9*Pi^2)/2)*
  PlusDistribution[1, 1 - x] + (27/(-1 + x) - 54/x - 135*x - (81*x^2)/2 + 
   189/(4*(1 + x)))*Zeta[3] + PlusDistribution[-1, 1 - x]*
  (33/2 + nl^2/9 - (67*Pi^2)/8 + (3*Pi^4)/8 + 
   nl*(-17/6 + (5*Pi^2)/12 - (9*Zeta[3])/8) + (297*Zeta[3])/16) + 
 PlusDistribution[0, 1 - x]*(1025/16 + (5*nl^2)/36 - (33*Pi^2)/16 + 
   nl*(-85/12 + Pi^2/8) + (81*Zeta[3])/4) + 
 nl*(-185/8 + 135/(16*x) + (1069*x)/24 - (1093*x^2)/48 + 
   Pi^2*(-5/4 - (13*x)/12 - x^2/4 + 1/(8*(1 + x))) + (-2 - 2*x)*Zeta[3]) + 
 ep*(-6743/32 + 1643/(8*x) + (4693*x)/32 - (1643*x^2)/8 + 
   nl^2*(1/24 + 61/(216*x) + x/24 - (109*x^2)/216) + 
   Pi^2*(-63/8 - 33/(4*x) - (21*x)/8 - (165*x^2)/8 - 33/(16*(1 + x))) + 
   Pi^2*(18 + 9/x + 9*x + 9*x^2 - 9/(1 + x))*HPL[{-1}, x] + 
   (4143/16 + nl^2*(1/4 + (13*x)/36) + 201/(4*(-1 + x)) - (519*x)/16 + 
     (483*x^2)/2 + nl*(-8 + Pi^2*(-1/3 - x/3) - 5/(2*(-1 + x)) + 61/(12*x) + 
       (41*x)/12 + 3*x^2) + Pi^2*(-9/(2*(-1 + x)) + 9*x - 9*x^2 + 
       9/(2*(1 + x))))*HPL[{0}, x] + (879/8 - 201/(4*x) - (75*x)/8 + 
     (201*x^2)/4 + nl*(-17/4 + 91/(12*x) + (13*x)/4 - (139*x^2)/12) + 
     Pi^2*(-9 + 9/(2*x) + (9*x)/2 - (9*x^2)/2))*HPL[{1}, x] + 
   (189/4 - 99/(4*(-1 + x)) - 99/(4*x) - (81*x)/2 + (495*x^2)/4 + 
     nl*(15/2 + 3/(2*(-1 + x)) - 3/(2*x) + 5*x + (3*x^2)/2))*HPL[{2}, x] + 
   (81/(2*(-1 + x)) - 27/x - 135*x + 54*x^2 - 27/(2*(1 + x)) + nl*(2 + 2*x))*
    HPL[{3}, x] + (-27/x - 81*x - 27*x^2 + 27/(1 + x))*HPL[{-2, 0}, x] + 
   (-63/2 - 297/(4*x) - (225*x)/4 - (297*x^2)/4 - 99/(4*(1 + x)) + 
     nl*(-3 - 3/(2*x) - (3*x)/2 - (3*x^2)/2 + 3/(2*(1 + x))))*
    HPL[{-1, 0}, x] + (-54 - 27/x - 27*x - 27*x^2 + 27/(1 + x))*
    HPL[{-1, 2}, x] + (261/4 + nl^2*(1/9 + x/9) - 99/(8*(-1 + x)) - 27*x + 
     (495*x^2)/4 + 99/(8*(1 + x)) + nl*(-16/3 + 3/(4*(-1 + x)) + (55*x)/6 + 
       (3*x^2)/2 - 3/(4*(1 + x))))*HPL[{0, 0}, x] + 
   (63/2 - 99/(4*(-1 + x)) - 297/(4*x) - (225*x)/4 + (297*x^2)/4 + 
     nl*(3 + 3/(2*(-1 + x)) - 3/(2*x) - (3*x)/2 + (3*x^2)/2))*
    HPL[{1, 0}, x] + (108 + 54/(-1 + x) - 54/x - 54*x + 54*x^2)*
    HPL[{1, 2}, x] + (54 + 54/(-1 + x) - 54/x - 108*x + 54*x^2)*
    HPL[{2, 0}, x] + (108 + 54/(-1 + x) - 54/x - 54*x + 54*x^2)*
    HPL[{2, 1}, x] + (108 + 54/x + 54*x + 54*x^2 - 54/(1 + x))*
    HPL[{-1, -1, 0}, x] + (-81 - 81/(2*x) - (81*x)/2 - (81*x^2)/2 + 
     81/(2*(1 + x)))*HPL[{-1, 0, 0}, x] + 
   (54 + 27/(2*(-1 + x)) - 108*x + 27*x^2 - 27/(2*(1 + x)) + nl*(-2 + 4*x))*
    HPL[{0, 0, 0}, x] + (81 + 81/(2*(-1 + x)) - 81/(2*x) - (81*x)/2 + 
     (81*x^2)/2)*HPL[{1, 0, 0}, x] + (108 + 54/(-1 + x) - 54/x - 54*x + 
     54*x^2)*HPL[{1, 1, 0}, x] + (201/2 - 5*nl - (9*Pi^2)/2)*
    PlusDistribution[1, 1 - x] + (27/(-1 + x) - 54/x - 135*x - (81*x^2)/2 + 
     189/(4*(1 + x)))*Zeta[3] + PlusDistribution[-1, 1 - x]*
    (33/2 + nl^2/9 - (67*Pi^2)/8 + (3*Pi^4)/8 + 
     nl*(-17/6 + (5*Pi^2)/12 - (9*Zeta[3])/8) + (297*Zeta[3])/16) + 
   PlusDistribution[0, 1 - x]*(1025/16 + (5*nl^2)/36 - (33*Pi^2)/16 + 
     nl*(-85/12 + Pi^2/8) + (81*Zeta[3])/4) + 
   nl*(-185/8 + 135/(16*x) + (1069*x)/24 - (1093*x^2)/48 + 
     Pi^2*(-5/4 - (13*x)/12 - x^2/4 + 1/(8*(1 + x))) + (-2 - 2*x)*Zeta[3])) + 
 ep^2*(-6743/32 + 1643/(8*x) + (4693*x)/32 - (1643*x^2)/8 + 
   nl^2*(1/24 + 61/(216*x) + x/24 - (109*x^2)/216) + 
   Pi^2*(-63/8 - 33/(4*x) - (21*x)/8 - (165*x^2)/8 - 33/(16*(1 + x))) + 
   Pi^2*(18 + 9/x + 9*x + 9*x^2 - 9/(1 + x))*HPL[{-1}, x] + 
   (4143/16 + nl^2*(1/4 + (13*x)/36) + 201/(4*(-1 + x)) - (519*x)/16 + 
     (483*x^2)/2 + nl*(-8 + Pi^2*(-1/3 - x/3) - 5/(2*(-1 + x)) + 61/(12*x) + 
       (41*x)/12 + 3*x^2) + Pi^2*(-9/(2*(-1 + x)) + 9*x - 9*x^2 + 
       9/(2*(1 + x))))*HPL[{0}, x] + (879/8 - 201/(4*x) - (75*x)/8 + 
     (201*x^2)/4 + nl*(-17/4 + 91/(12*x) + (13*x)/4 - (139*x^2)/12) + 
     Pi^2*(-9 + 9/(2*x) + (9*x)/2 - (9*x^2)/2))*HPL[{1}, x] + 
   (189/4 - 99/(4*(-1 + x)) - 99/(4*x) - (81*x)/2 + (495*x^2)/4 + 
     nl*(15/2 + 3/(2*(-1 + x)) - 3/(2*x) + 5*x + (3*x^2)/2))*HPL[{2}, x] + 
   (81/(2*(-1 + x)) - 27/x - 135*x + 54*x^2 - 27/(2*(1 + x)) + nl*(2 + 2*x))*
    HPL[{3}, x] + (-27/x - 81*x - 27*x^2 + 27/(1 + x))*HPL[{-2, 0}, x] + 
   (-63/2 - 297/(4*x) - (225*x)/4 - (297*x^2)/4 - 99/(4*(1 + x)) + 
     nl*(-3 - 3/(2*x) - (3*x)/2 - (3*x^2)/2 + 3/(2*(1 + x))))*
    HPL[{-1, 0}, x] + (-54 - 27/x - 27*x - 27*x^2 + 27/(1 + x))*
    HPL[{-1, 2}, x] + (261/4 + nl^2*(1/9 + x/9) - 99/(8*(-1 + x)) - 27*x + 
     (495*x^2)/4 + 99/(8*(1 + x)) + nl*(-16/3 + 3/(4*(-1 + x)) + (55*x)/6 + 
       (3*x^2)/2 - 3/(4*(1 + x))))*HPL[{0, 0}, x] + 
   (63/2 - 99/(4*(-1 + x)) - 297/(4*x) - (225*x)/4 + (297*x^2)/4 + 
     nl*(3 + 3/(2*(-1 + x)) - 3/(2*x) - (3*x)/2 + (3*x^2)/2))*
    HPL[{1, 0}, x] + (108 + 54/(-1 + x) - 54/x - 54*x + 54*x^2)*
    HPL[{1, 2}, x] + (54 + 54/(-1 + x) - 54/x - 108*x + 54*x^2)*
    HPL[{2, 0}, x] + (108 + 54/(-1 + x) - 54/x - 54*x + 54*x^2)*
    HPL[{2, 1}, x] + (108 + 54/x + 54*x + 54*x^2 - 54/(1 + x))*
    HPL[{-1, -1, 0}, x] + (-81 - 81/(2*x) - (81*x)/2 - (81*x^2)/2 + 
     81/(2*(1 + x)))*HPL[{-1, 0, 0}, x] + 
   (54 + 27/(2*(-1 + x)) - 108*x + 27*x^2 - 27/(2*(1 + x)) + nl*(-2 + 4*x))*
    HPL[{0, 0, 0}, x] + (81 + 81/(2*(-1 + x)) - 81/(2*x) - (81*x)/2 + 
     (81*x^2)/2)*HPL[{1, 0, 0}, x] + (108 + 54/(-1 + x) - 54/x - 54*x + 
     54*x^2)*HPL[{1, 1, 0}, x] + (201/2 - 5*nl - (9*Pi^2)/2)*
    PlusDistribution[1, 1 - x] + (27/(-1 + x) - 54/x - 135*x - (81*x^2)/2 + 
     189/(4*(1 + x)))*Zeta[3] + PlusDistribution[-1, 1 - x]*
    (33/2 + nl^2/9 - (67*Pi^2)/8 + (3*Pi^4)/8 + 
     nl*(-17/6 + (5*Pi^2)/12 - (9*Zeta[3])/8) + (297*Zeta[3])/16) + 
   PlusDistribution[0, 1 - x]*(1025/16 + (5*nl^2)/36 - (33*Pi^2)/16 + 
     nl*(-85/12 + Pi^2/8) + (81*Zeta[3])/4) + 
   nl*(-185/8 + 135/(16*x) + (1069*x)/24 - (1093*x^2)/48 + 
     Pi^2*(-5/4 - (13*x)/12 - x^2/4 + 1/(8*(1 + x))) + (-2 - 2*x)*Zeta[3])) + 
 ep^3*(-6743/32 + 1643/(8*x) + (4693*x)/32 - (1643*x^2)/8 + 
   nl^2*(1/24 + 61/(216*x) + x/24 - (109*x^2)/216) + 
   Pi^2*(-63/8 - 33/(4*x) - (21*x)/8 - (165*x^2)/8 - 33/(16*(1 + x))) + 
   Pi^2*(18 + 9/x + 9*x + 9*x^2 - 9/(1 + x))*HPL[{-1}, x] + 
   (4143/16 + nl^2*(1/4 + (13*x)/36) + 201/(4*(-1 + x)) - (519*x)/16 + 
     (483*x^2)/2 + nl*(-8 + Pi^2*(-1/3 - x/3) - 5/(2*(-1 + x)) + 61/(12*x) + 
       (41*x)/12 + 3*x^2) + Pi^2*(-9/(2*(-1 + x)) + 9*x - 9*x^2 + 
       9/(2*(1 + x))))*HPL[{0}, x] + (879/8 - 201/(4*x) - (75*x)/8 + 
     (201*x^2)/4 + nl*(-17/4 + 91/(12*x) + (13*x)/4 - (139*x^2)/12) + 
     Pi^2*(-9 + 9/(2*x) + (9*x)/2 - (9*x^2)/2))*HPL[{1}, x] + 
   (189/4 - 99/(4*(-1 + x)) - 99/(4*x) - (81*x)/2 + (495*x^2)/4 + 
     nl*(15/2 + 3/(2*(-1 + x)) - 3/(2*x) + 5*x + (3*x^2)/2))*HPL[{2}, x] + 
   (81/(2*(-1 + x)) - 27/x - 135*x + 54*x^2 - 27/(2*(1 + x)) + nl*(2 + 2*x))*
    HPL[{3}, x] + (-27/x - 81*x - 27*x^2 + 27/(1 + x))*HPL[{-2, 0}, x] + 
   (-63/2 - 297/(4*x) - (225*x)/4 - (297*x^2)/4 - 99/(4*(1 + x)) + 
     nl*(-3 - 3/(2*x) - (3*x)/2 - (3*x^2)/2 + 3/(2*(1 + x))))*
    HPL[{-1, 0}, x] + (-54 - 27/x - 27*x - 27*x^2 + 27/(1 + x))*
    HPL[{-1, 2}, x] + (261/4 + nl^2*(1/9 + x/9) - 99/(8*(-1 + x)) - 27*x + 
     (495*x^2)/4 + 99/(8*(1 + x)) + nl*(-16/3 + 3/(4*(-1 + x)) + (55*x)/6 + 
       (3*x^2)/2 - 3/(4*(1 + x))))*HPL[{0, 0}, x] + 
   (63/2 - 99/(4*(-1 + x)) - 297/(4*x) - (225*x)/4 + (297*x^2)/4 + 
     nl*(3 + 3/(2*(-1 + x)) - 3/(2*x) - (3*x)/2 + (3*x^2)/2))*
    HPL[{1, 0}, x] + (108 + 54/(-1 + x) - 54/x - 54*x + 54*x^2)*
    HPL[{1, 2}, x] + (54 + 54/(-1 + x) - 54/x - 108*x + 54*x^2)*
    HPL[{2, 0}, x] + (108 + 54/(-1 + x) - 54/x - 54*x + 54*x^2)*
    HPL[{2, 1}, x] + (108 + 54/x + 54*x + 54*x^2 - 54/(1 + x))*
    HPL[{-1, -1, 0}, x] + (-81 - 81/(2*x) - (81*x)/2 - (81*x^2)/2 + 
     81/(2*(1 + x)))*HPL[{-1, 0, 0}, x] + 
   (54 + 27/(2*(-1 + x)) - 108*x + 27*x^2 - 27/(2*(1 + x)) + nl*(-2 + 4*x))*
    HPL[{0, 0, 0}, x] + (81 + 81/(2*(-1 + x)) - 81/(2*x) - (81*x)/2 + 
     (81*x^2)/2)*HPL[{1, 0, 0}, x] + (108 + 54/(-1 + x) - 54/x - 54*x + 
     54*x^2)*HPL[{1, 1, 0}, x] + (201/2 - 5*nl - (9*Pi^2)/2)*
    PlusDistribution[1, 1 - x] + (27/(-1 + x) - 54/x - 135*x - (81*x^2)/2 + 
     189/(4*(1 + x)))*Zeta[3] + PlusDistribution[-1, 1 - x]*
    (33/2 + nl^2/9 - (67*Pi^2)/8 + (3*Pi^4)/8 + 
     nl*(-17/6 + (5*Pi^2)/12 - (9*Zeta[3])/8) + (297*Zeta[3])/16) + 
   PlusDistribution[0, 1 - x]*(1025/16 + (5*nl^2)/36 - (33*Pi^2)/16 + 
     nl*(-85/12 + Pi^2/8) + (81*Zeta[3])/4) + 
   nl*(-185/8 + 135/(16*x) + (1069*x)/24 - (1093*x^2)/48 + 
     Pi^2*(-5/4 - (13*x)/12 - x^2/4 + 1/(8*(1 + x))) + (-2 - 2*x)*Zeta[3]))
