U1[x,ep] = (
  + 1/2 - x^2/2 + x*HPL[{0}, x]
  + ep*(13/4 - (2*Pi^2*x)/3 - (13*x^2)/4 + 3*x*HPL[{0}, x] - (x^2*HPL[{0}, x])/2 + 2*HPL[{1}, x] - 2*x^2*HPL[{1}, x] + 4*x*HPL[{2}, x] - x*HPL[{0, 0}, x])
  + ep^2*(115/8 - Pi^2/2 - 2*Pi^2*x - (115*x^2)/8 + (5*Pi^2*x^2)/6 + 6*x*HPL[{0}, x] - (Pi^2*x*HPL[{0}, x])/3 - (13*x^2*HPL[{0}, x])/4 + 13*HPL[{1}, x] - 13*x^2*HPL[{1}, x] + 12*x*HPL[{2}, x] - 2*x^2*HPL[{2}, x] - 4*x*HPL[{3}, x] - 3*x*HPL[{0, 0}, x] + (x^2*HPL[{0, 0}, x])/2 + HPL[{1, 0}, x] - x^2*HPL[{1, 0}, x] + 8*HPL[{1, 1}, x] - 8*x^2*HPL[{1, 1}, x] + 2*x*HPL[{2, 0}, x] + 16*x*HPL[{2, 1}, x] + x*HPL[{0, 0, 0}, x] - 8*x*Zeta[3])
  + ep^3*(865/16 - (13*Pi^2)/4 - 4*Pi^2*x + (2*Pi^4*x)/9 - (865*x^2)/16 + (65*Pi^2*x^2)/12 + (15*x*HPL[{0}, x])/2 - Pi^2*x*HPL[{0}, x] - (115*x^2*HPL[{0}, x])/8 + (Pi^2*x^2*HPL[{0}, x])/6 + (115*HPL[{1}, x])/2 - (8*Pi^2*HPL[{1}, x])/3 - (115*x^2*HPL[{1}, x])/2 + (8*Pi^2*x^2*HPL[{1}, x])/3 + 24*x*HPL[{2}, x] - (16*Pi^2*x*HPL[{2}, x])/3 - 13*x^2*HPL[{2}, x] - 12*x*HPL[{3}, x] + 2*x^2*HPL[{3}, x] + 4*x*HPL[{4}, x] - 6*x*HPL[{0, 0}, x] + (Pi^2*x*HPL[{0, 0}, x])/3 + (13*x^2*HPL[{0, 0}, x])/4 + (13*HPL[{1, 0}, x])/2 - (13*x^2*HPL[{1, 0}, x])/2 + 52*HPL[{1, 1}, x] - 52*x^2*HPL[{1, 1}, x] + 4*HPL[{1, 2}, x] - 4*x^2*HPL[{1, 2}, x] + 6*x*HPL[{2, 0}, x] - x^2*HPL[{2, 0}, x] + 48*x*HPL[{2, 1}, x] - 8*x^2*HPL[{2, 1}, x] + 8*x*HPL[{2, 2}, x] - 2*x*HPL[{3, 0}, x] - 16*x*HPL[{3, 1}, x] + 3*x*HPL[{0, 0, 0}, x] - (x^2*HPL[{0, 0, 0}, x])/2 - HPL[{1, 0, 0}, x] + x^2*HPL[{1, 0, 0}, x] + 4*HPL[{1, 1, 0}, x] - 4*x^2*HPL[{1, 1, 0}, x] + 32*HPL[{1, 1, 1}, x] - 32*x^2*HPL[{1, 1, 1}, x] - 2*x*HPL[{2, 0, 0}, x] + 8*x*HPL[{2, 1, 0}, x] + 64*x*HPL[{2, 1, 1}, x] - x*HPL[{0, 0, 0, 0}, x] - 5*Zeta[3] - 24*x*Zeta[3] + 9*x^2*Zeta[3] - 2*x*HPL[{0}, x]*Zeta[3])
  + ep^4*(5971/32 - (115*Pi^2)/8 + (7*Pi^4)/60 - 5*Pi^2*x + (2*Pi^4*x)/3 - (5971*x^2)/32 + (575*Pi^2*x^2)/24 - (41*Pi^4*x^2)/180 - (21*x*HPL[{0}, x])/4 - 2*Pi^2*x*HPL[{0}, x] + (Pi^4*x*HPL[{0}, x])/90 - (865*x^2*HPL[{0}, x])/16 + (13*Pi^2*x^2*HPL[{0}, x])/12 + (865*HPL[{1}, x])/4 - (52*Pi^2*HPL[{1}, x])/3 - (865*x^2*HPL[{1}, x])/4 + (52*Pi^2*x^2*HPL[{1}, x])/3 + 30*x*HPL[{2}, x] - 16*Pi^2*x*HPL[{2}, x] - (115*x^2*HPL[{2}, x])/2 + (8*Pi^2*x^2*HPL[{2}, x])/3 - 24*x*HPL[{3}, x] + (16*Pi^2*x*HPL[{3}, x])/3 + 13*x^2*HPL[{3}, x] + 12*x*HPL[{4}, x] - 2*x^2*HPL[{4}, x] - 4*x*HPL[{5}, x] - (15*x*HPL[{0, 0}, x])/2 + Pi^2*x*HPL[{0, 0}, x] + (115*x^2*HPL[{0, 0}, x])/8 - (Pi^2*x^2*HPL[{0, 0}, x])/6 + (115*HPL[{1, 0}, x])/4 - (Pi^2*HPL[{1, 0}, x])/3 - (115*x^2*HPL[{1, 0}, x])/4 + (Pi^2*x^2*HPL[{1, 0}, x])/3 + 230*HPL[{1, 1}, x] - (32*Pi^2*HPL[{1, 1}, x])/3 - 230*x^2*HPL[{1, 1}, x] + (32*Pi^2*x^2*HPL[{1, 1}, x])/3 + 26*HPL[{1, 2}, x] - 26*x^2*HPL[{1, 2}, x] - 4*HPL[{1, 3}, x] + 4*x^2*HPL[{1, 3}, x] + 12*x*HPL[{2, 0}, x] - (2*Pi^2*x*HPL[{2, 0}, x])/3 - (13*x^2*HPL[{2, 0}, x])/2 + 96*x*HPL[{2, 1}, x] - (64*Pi^2*x*HPL[{2, 1}, x])/3 - 52*x^2*HPL[{2, 1}, x] + 24*x*HPL[{2, 2}, x] - 4*x^2*HPL[{2, 2}, x] - 8*x*HPL[{2, 3}, x] - 6*x*HPL[{3, 0}, x] + x^2*HPL[{3, 0}, x] - 48*x*HPL[{3, 1}, x] + 8*x^2*HPL[{3, 1}, x] - 8*x*HPL[{3, 2}, x] + 2*x*HPL[{4, 0}, x] + 16*x*HPL[{4, 1}, x] + 6*x*HPL[{0, 0, 0}, x] - (Pi^2*x*HPL[{0, 0, 0}, x])/3 - (13*x^2*HPL[{0, 0, 0}, x])/4 - (13*HPL[{1, 0, 0}, x])/2 + (13*x^2*HPL[{1, 0, 0}, x])/2 + 26*HPL[{1, 1, 0}, x] - 26*x^2*HPL[{1, 1, 0}, x] + 208*HPL[{1, 1, 1}, x] - 208*x^2*HPL[{1, 1, 1}, x] + 16*HPL[{1, 1, 2}, x] - 16*x^2*HPL[{1, 1, 2}, x] + 2*HPL[{1, 2, 0}, x] - 2*x^2*HPL[{1, 2, 0}, x] + 16*HPL[{1, 2, 1}, x] - 16*x^2*HPL[{1, 2, 1}, x] - 6*x*HPL[{2, 0, 0}, x] + x^2*HPL[{2, 0, 0}, x] + 24*x*HPL[{2, 1, 0}, x] - 4*x^2*HPL[{2, 1, 0}, x] + 192*x*HPL[{2, 1, 1}, x] - 32*x^2*HPL[{2, 1, 1}, x] + 32*x*HPL[{2, 1, 2}, x] + 4*x*HPL[{2, 2, 0}, x] + 32*x*HPL[{2, 2, 1}, x] + 2*x*HPL[{3, 0, 0}, x] - 8*x*HPL[{3, 1, 0}, x] - 64*x*HPL[{3, 1, 1}, x] - 3*x*HPL[{0, 0, 0, 0}, x] + (x^2*HPL[{0, 0, 0, 0}, x])/2 + HPL[{1, 0, 0, 0}, x] - x^2*HPL[{1, 0, 0, 0}, x] - 4*HPL[{1, 1, 0, 0}, x] + 4*x^2*HPL[{1, 1, 0, 0}, x] + 16*HPL[{1, 1, 1, 0}, x] - 16*x^2*HPL[{1, 1, 1, 0}, x] + 128*HPL[{1, 1, 1, 1}, x] - 128*x^2*HPL[{1, 1, 1, 1}, x] + 2*x*HPL[{2, 0, 0, 0}, x] - 8*x*HPL[{2, 1, 0, 0}, x] + 32*x*HPL[{2, 1, 1, 0}, x] + 256*x*HPL[{2, 1, 1, 1}, x] + x*HPL[{0, 0, 0, 0, 0}, x] - (65*Zeta[3])/2 - 48*x*Zeta[3] + (28*Pi^2*x*Zeta[3])/3 + (117*x^2*Zeta[3])/2 - 6*x*HPL[{0}, x]*Zeta[3] + x^2*HPL[{0}, x]*Zeta[3] - 28*HPL[{1}, x]*Zeta[3] + 28*x^2*HPL[{1}, x]*Zeta[3] - 56*x*HPL[{2}, x]*Zeta[3] + 2*x*HPL[{0, 0}, x]*Zeta[3] - 48*x*Zeta[5])
  + ep^5*(39193/64 - (865*Pi^2)/16 + (91*Pi^4)/120 + (7*Pi^2*x)/2 + (4*Pi^4*x)/3 - (Pi^6*x)/27 - (39193*x^2)/64 + (4325*Pi^2*x^2)/48 - (533*Pi^4*x^2)/360 - (651*x*HPL[{0}, x])/8 - (5*Pi^2*x*HPL[{0}, x])/2 + (Pi^4*x*HPL[{0}, x])/30 - (5971*x^2*HPL[{0}, x])/32 + (115*Pi^2*x^2*HPL[{0}, x])/24 - (Pi^4*x^2*HPL[{0}, x])/180 + (5971*HPL[{1}, x])/8 - (230*Pi^2*HPL[{1}, x])/3 + (31*Pi^4*HPL[{1}, x])/45 - (5971*x^2*HPL[{1}, x])/8 + (230*Pi^2*x^2*HPL[{1}, x])/3 - (31*Pi^4*x^2*HPL[{1}, x])/45 - 21*x*HPL[{2}, x] - 32*Pi^2*x*HPL[{2}, x] + (62*Pi^4*x*HPL[{2}, x])/45 - (865*x^2*HPL[{2}, x])/4 + (52*Pi^2*x^2*HPL[{2}, x])/3 - 30*x*HPL[{3}, x] + 16*Pi^2*x*HPL[{3}, x] + (115*x^2*HPL[{3}, x])/2 - (8*Pi^2*x^2*HPL[{3}, x])/3 + 24*x*HPL[{4}, x] - (16*Pi^2*x*HPL[{4}, x])/3 - 13*x^2*HPL[{4}, x] - 12*x*HPL[{5}, x] + 2*x^2*HPL[{5}, x] + 4*x*HPL[{6}, x] + (21*x*HPL[{0, 0}, x])/4 + 2*Pi^2*x*HPL[{0, 0}, x] - (Pi^4*x*HPL[{0, 0}, x])/90 + (865*x^2*HPL[{0, 0}, x])/16 - (13*Pi^2*x^2*HPL[{0, 0}, x])/12 + (865*HPL[{1, 0}, x])/8 - (13*Pi^2*HPL[{1, 0}, x])/6 - (865*x^2*HPL[{1, 0}, x])/8 + (13*Pi^2*x^2*HPL[{1, 0}, x])/6 + 865*HPL[{1, 1}, x] - (208*Pi^2*HPL[{1, 1}, x])/3 - 865*x^2*HPL[{1, 1}, x] + (208*Pi^2*x^2*HPL[{1, 1}, x])/3 + 115*HPL[{1, 2}, x] - (16*Pi^2*HPL[{1, 2}, x])/3 - 115*x^2*HPL[{1, 2}, x] + (16*Pi^2*x^2*HPL[{1, 2}, x])/3 - 26*HPL[{1, 3}, x] + 26*x^2*HPL[{1, 3}, x] + 4*HPL[{1, 4}, x] - 4*x^2*HPL[{1, 4}, x] + 15*x*HPL[{2, 0}, x] - 2*Pi^2*x*HPL[{2, 0}, x] - (115*x^2*HPL[{2, 0}, x])/4 + (Pi^2*x^2*HPL[{2, 0}, x])/3 + 120*x*HPL[{2, 1}, x] - 64*Pi^2*x*HPL[{2, 1}, x] - 230*x^2*HPL[{2, 1}, x] + (32*Pi^2*x^2*HPL[{2, 1}, x])/3 + 48*x*HPL[{2, 2}, x] - (32*Pi^2*x*HPL[{2, 2}, x])/3 - 26*x^2*HPL[{2, 2}, x] - 24*x*HPL[{2, 3}, x] + 4*x^2*HPL[{2, 3}, x] + 8*x*HPL[{2, 4}, x] - 12*x*HPL[{3, 0}, x] + (2*Pi^2*x*HPL[{3, 0}, x])/3 + (13*x^2*HPL[{3, 0}, x])/2 - 96*x*HPL[{3, 1}, x] + (64*Pi^2*x*HPL[{3, 1}, x])/3 + 52*x^2*HPL[{3, 1}, x] - 24*x*HPL[{3, 2}, x] + 4*x^2*HPL[{3, 2}, x] + 8*x*HPL[{3, 3}, x] + 6*x*HPL[{4, 0}, x] - x^2*HPL[{4, 0}, x] + 48*x*HPL[{4, 1}, x] - 8*x^2*HPL[{4, 1}, x] + 8*x*HPL[{4, 2}, x] - 2*x*HPL[{5, 0}, x] - 16*x*HPL[{5, 1}, x] + (15*x*HPL[{0, 0, 0}, x])/2 - Pi^2*x*HPL[{0, 0, 0}, x] - (115*x^2*HPL[{0, 0, 0}, x])/8 + (Pi^2*x^2*HPL[{0, 0, 0}, x])/6 - (115*HPL[{1, 0, 0}, x])/4 + (Pi^2*HPL[{1, 0, 0}, x])/3 + (115*x^2*HPL[{1, 0, 0}, x])/4 - (Pi^2*x^2*HPL[{1, 0, 0}, x])/3 + 115*HPL[{1, 1, 0}, x] - (4*Pi^2*HPL[{1, 1, 0}, x])/3 - 115*x^2*HPL[{1, 1, 0}, x] + (4*Pi^2*x^2*HPL[{1, 1, 0}, x])/3 + 920*HPL[{1, 1, 1}, x] - (128*Pi^2*HPL[{1, 1, 1}, x])/3 - 920*x^2*HPL[{1, 1, 1}, x] + (128*Pi^2*x^2*HPL[{1, 1, 1}, x])/3 + 104*HPL[{1, 1, 2}, x] - 104*x^2*HPL[{1, 1, 2}, x] - 16*HPL[{1, 1, 3}, x] + 16*x^2*HPL[{1, 1, 3}, x] + 13*HPL[{1, 2, 0}, x] - 13*x^2*HPL[{1, 2, 0}, x] + 104*HPL[{1, 2, 1}, x] - 104*x^2*HPL[{1, 2, 1}, x] + 8*HPL[{1, 2, 2}, x] - 8*x^2*HPL[{1, 2, 2}, x] - 2*HPL[{1, 3, 0}, x] + 2*x^2*HPL[{1, 3, 0}, x] - 16*HPL[{1, 3, 1}, x] + 16*x^2*HPL[{1, 3, 1}, x] - 12*x*HPL[{2, 0, 0}, x] + (2*Pi^2*x*HPL[{2, 0, 0}, x])/3 + (13*x^2*HPL[{2, 0, 0}, x])/2 + 48*x*HPL[{2, 1, 0}, x] - (8*Pi^2*x*HPL[{2, 1, 0}, x])/3 - 26*x^2*HPL[{2, 1, 0}, x] + 384*x*HPL[{2, 1, 1}, x] - (256*Pi^2*x*HPL[{2, 1, 1}, x])/3 - 208*x^2*HPL[{2, 1, 1}, x] + 96*x*HPL[{2, 1, 2}, x] - 16*x^2*HPL[{2, 1, 2}, x] - 32*x*HPL[{2, 1, 3}, x] + 12*x*HPL[{2, 2, 0}, x] - 2*x^2*HPL[{2, 2, 0}, x] + 96*x*HPL[{2, 2, 1}, x] - 16*x^2*HPL[{2, 2, 1}, x] + 16*x*HPL[{2, 2, 2}, x] - 4*x*HPL[{2, 3, 0}, x] - 32*x*HPL[{2, 3, 1}, x] + 6*x*HPL[{3, 0, 0}, x] - x^2*HPL[{3, 0, 0}, x] - 24*x*HPL[{3, 1, 0}, x] + 4*x^2*HPL[{3, 1, 0}, x] - 192*x*HPL[{3, 1, 1}, x] + 32*x^2*HPL[{3, 1, 1}, x] - 32*x*HPL[{3, 1, 2}, x] - 4*x*HPL[{3, 2, 0}, x] - 32*x*HPL[{3, 2, 1}, x] - 2*x*HPL[{4, 0, 0}, x] + 8*x*HPL[{4, 1, 0}, x] + 64*x*HPL[{4, 1, 1}, x] - 6*x*HPL[{0, 0, 0, 0}, x] + (Pi^2*x*HPL[{0, 0, 0, 0}, x])/3 + (13*x^2*HPL[{0, 0, 0, 0}, x])/4 + (13*HPL[{1, 0, 0, 0}, x])/2 - (13*x^2*HPL[{1, 0, 0, 0}, x])/2 - 26*HPL[{1, 1, 0, 0}, x] + 26*x^2*HPL[{1, 1, 0, 0}, x] + 104*HPL[{1, 1, 1, 0}, x] - 104*x^2*HPL[{1, 1, 1, 0}, x] + 832*HPL[{1, 1, 1, 1}, x] - 832*x^2*HPL[{1, 1, 1, 1}, x] + 64*HPL[{1, 1, 1, 2}, x] - 64*x^2*HPL[{1, 1, 1, 2}, x] + 8*HPL[{1, 1, 2, 0}, x] - 8*x^2*HPL[{1, 1, 2, 0}, x] + 64*HPL[{1, 1, 2, 1}, x] - 64*x^2*HPL[{1, 1, 2, 1}, x] - 2*HPL[{1, 2, 0, 0}, x] + 2*x^2*HPL[{1, 2, 0, 0}, x] + 8*HPL[{1, 2, 1, 0}, x] - 8*x^2*HPL[{1, 2, 1, 0}, x] + 64*HPL[{1, 2, 1, 1}, x] - 64*x^2*HPL[{1, 2, 1, 1}, x] + 6*x*HPL[{2, 0, 0, 0}, x] - x^2*HPL[{2, 0, 0, 0}, x] - 24*x*HPL[{2, 1, 0, 0}, x] + 4*x^2*HPL[{2, 1, 0, 0}, x] + 96*x*HPL[{2, 1, 1, 0}, x] - 16*x^2*HPL[{2, 1, 1, 0}, x] + 768*x*HPL[{2, 1, 1, 1}, x] - 128*x^2*HPL[{2, 1, 1, 1}, x] + 128*x*HPL[{2, 1, 1, 2}, x] + 16*x*HPL[{2, 1, 2, 0}, x] + 128*x*HPL[{2, 1, 2, 1}, x] - 4*x*HPL[{2, 2, 0, 0}, x] + 16*x*HPL[{2, 2, 1, 0}, x] + 128*x*HPL[{2, 2, 1, 1}, x] - 2*x*HPL[{3, 0, 0, 0}, x] + 8*x*HPL[{3, 1, 0, 0}, x] - 32*x*HPL[{3, 1, 1, 0}, x] - 256*x*HPL[{3, 1, 1, 1}, x] + 3*x*HPL[{0, 0, 0, 0, 0}, x] - (x^2*HPL[{0, 0, 0, 0, 0}, x])/2 - HPL[{1, 0, 0, 0, 0}, x] + x^2*HPL[{1, 0, 0, 0, 0}, x] + 4*HPL[{1, 1, 0, 0, 0}, x] - 4*x^2*HPL[{1, 1, 0, 0, 0}, x] - 16*HPL[{1, 1, 1, 0, 0}, x] + 16*x^2*HPL[{1, 1, 1, 0, 0}, x] + 64*HPL[{1, 1, 1, 1, 0}, x] - 64*x^2*HPL[{1, 1, 1, 1, 0}, x] + 512*HPL[{1, 1, 1, 1, 1}, x] - 512*x^2*HPL[{1, 1, 1, 1, 1}, x] - 2*x*HPL[{2, 0, 0, 0, 0}, x] + 8*x*HPL[{2, 1, 0, 0, 0}, x] - 32*x*HPL[{2, 1, 1, 0, 0}, x] + 128*x*HPL[{2, 1, 1, 1, 0}, x] + 1024*x*HPL[{2, 1, 1, 1, 1}, x] - x*HPL[{0, 0, 0, 0, 0, 0}, x] - (575*Zeta[3])/4 + 5*Pi^2*Zeta[3] - 60*x*Zeta[3] + 28*Pi^2*x*Zeta[3] + (1035*x^2*Zeta[3])/4 - (29*Pi^2*x^2*Zeta[3])/3 - 12*x*HPL[{0}, x]*Zeta[3] + (2*Pi^2*x*HPL[{0}, x]*Zeta[3])/3 + (13*x^2*HPL[{0}, x]*Zeta[3])/2 - 182*HPL[{1}, x]*Zeta[3] + 182*x^2*HPL[{1}, x]*Zeta[3] - 168*x*HPL[{2}, x]*Zeta[3] + 28*x^2*HPL[{2}, x]*Zeta[3] + 56*x*HPL[{3}, x]*Zeta[3] + 6*x*HPL[{0, 0}, x]*Zeta[3] - x^2*HPL[{0, 0}, x]*Zeta[3] - 2*HPL[{1, 0}, x]*Zeta[3] + 2*x^2*HPL[{1, 0}, x]*Zeta[3] - 112*HPL[{1, 1}, x]*Zeta[3] + 112*x^2*HPL[{1, 1}, x]*Zeta[3] - 4*x*HPL[{2, 0}, x]*Zeta[3] - 224*x*HPL[{2, 1}, x]*Zeta[3] - 2*x*HPL[{0, 0, 0}, x]*Zeta[3] + 48*x*Zeta[3]^2 - 27*Zeta[5] - 144*x*Zeta[5] + 51*x^2*Zeta[5] - 6*x*HPL[{0}, x]*Zeta[5])
  + ALARM*ep^6 );
U1[x -> 1,ep] = (1-x)^(3 - 4*ep) * (
  + 1/6 + (11*ep)/9 + ep^2*(170/27 - (2*Pi^2)/9) + ep^3*(2300/81 - (44*Pi^2)/27 - (10*Zeta[3])/3) + ep^4*(29288/243 - (680*Pi^2)/81 + (4*Pi^4)/135 - (220*Zeta[3])/9) + ep^5*(362096/729 - (9200*Pi^2)/243 + (88*Pi^4)/405 - (3400*Zeta[3])/27 + (40*Pi^2*Zeta[3])/9 - 34*Zeta[5])
  + ALARM*ep^6 );

U1a[x,ep] = (
  + 1/6 + x/4 - x^2/2 + x^3/12 + (x*HPL[{0}, x])/2
  + ep*(13/12 + (35*x)/24 - (Pi^2*x)/3 - (37*x^2)/12 + (13*x^3)/24 + (3*x*HPL[{0}, x])/2 - (x^2*HPL[{0}, x])/2 + (x^3*HPL[{0}, x])/12 + (2*HPL[{1}, x])/3 + x*HPL[{1}, x] - 2*x^2*HPL[{1}, x] + (x^3*HPL[{1}, x])/3 + 2*x*HPL[{2}, x] - (x*HPL[{0, 0}, x])/2)
  + ep^2*(115/24 - Pi^2/6 + (95*x)/16 - (5*Pi^2*x)/4 - (105*x^2)/8 + (5*Pi^2*x^2)/6 + (115*x^3)/48 - (5*Pi^2*x^3)/36 + (19*x*HPL[{0}, x])/6 - (Pi^2*x*HPL[{0}, x])/6 - (35*x^2*HPL[{0}, x])/12 + (13*x^3*HPL[{0}, x])/24 + (13*HPL[{1}, x])/3 + (35*x*HPL[{1}, x])/6 - (37*x^2*HPL[{1}, x])/3 + (13*x^3*HPL[{1}, x])/6 + 6*x*HPL[{2}, x] - 2*x^2*HPL[{2}, x] + (x^3*HPL[{2}, x])/3 - 2*x*HPL[{3}, x] - (3*x*HPL[{0, 0}, x])/2 + (x^2*HPL[{0, 0}, x])/2 - (x^3*HPL[{0, 0}, x])/12 + HPL[{1, 0}, x]/3 + (x*HPL[{1, 0}, x])/2 - x^2*HPL[{1, 0}, x] + (x^3*HPL[{1, 0}, x])/6 + (8*HPL[{1, 1}, x])/3 + 4*x*HPL[{1, 1}, x] - 8*x^2*HPL[{1, 1}, x] + (4*x^3*HPL[{1, 1}, x])/3 + x*HPL[{2, 0}, x] + 8*x*HPL[{2, 1}, x] + (x*HPL[{0, 0, 0}, x])/2 - 4*x*Zeta[3])
  + ep^3*(865/48 - (13*Pi^2)/12 + (2015*x)/96 - (257*Pi^2*x)/72 + (Pi^4*x)/9 - (2305*x^2)/48 + (181*Pi^2*x^2)/36 + (865*x^3)/96 - (65*Pi^2*x^3)/72 + 5*x*HPL[{0}, x] - (Pi^2*x*HPL[{0}, x])/2 - (95*x^2*HPL[{0}, x])/8 + (Pi^2*x^2*HPL[{0}, x])/6 + (115*x^3*HPL[{0}, x])/48 - (Pi^2*x^3*HPL[{0}, x])/36 + (115*HPL[{1}, x])/6 - (8*Pi^2*HPL[{1}, x])/9 + (95*x*HPL[{1}, x])/4 - (4*Pi^2*x*HPL[{1}, x])/3 - (105*x^2*HPL[{1}, x])/2 + (8*Pi^2*x^2*HPL[{1}, x])/3 + (115*x^3*HPL[{1}, x])/12 - (4*Pi^2*x^3*HPL[{1}, x])/9 + (38*x*HPL[{2}, x])/3 - (8*Pi^2*x*HPL[{2}, x])/3 - (35*x^2*HPL[{2}, x])/3 + (13*x^3*HPL[{2}, x])/6 - 6*x*HPL[{3}, x] + 2*x^2*HPL[{3}, x] - (x^3*HPL[{3}, x])/3 + 2*x*HPL[{4}, x] - (19*x*HPL[{0, 0}, x])/6 + (Pi^2*x*HPL[{0, 0}, x])/6 + (35*x^2*HPL[{0, 0}, x])/12 - (13*x^3*HPL[{0, 0}, x])/24 + (13*HPL[{1, 0}, x])/6 + (35*x*HPL[{1, 0}, x])/12 - (37*x^2*HPL[{1, 0}, x])/6 + (13*x^3*HPL[{1, 0}, x])/12 + (52*HPL[{1, 1}, x])/3 + (70*x*HPL[{1, 1}, x])/3 - (148*x^2*HPL[{1, 1}, x])/3 + (26*x^3*HPL[{1, 1}, x])/3 + (4*HPL[{1, 2}, x])/3 + 2*x*HPL[{1, 2}, x] - 4*x^2*HPL[{1, 2}, x] + (2*x^3*HPL[{1, 2}, x])/3 + 3*x*HPL[{2, 0}, x] - x^2*HPL[{2, 0}, x] + (x^3*HPL[{2, 0}, x])/6 + 24*x*HPL[{2, 1}, x] - 8*x^2*HPL[{2, 1}, x] + (4*x^3*HPL[{2, 1}, x])/3 + 4*x*HPL[{2, 2}, x] - x*HPL[{3, 0}, x] - 8*x*HPL[{3, 1}, x] + (3*x*HPL[{0, 0, 0}, x])/2 - (x^2*HPL[{0, 0, 0}, x])/2 + (x^3*HPL[{0, 0, 0}, x])/12 - HPL[{1, 0, 0}, x]/3 - (x*HPL[{1, 0, 0}, x])/2 + x^2*HPL[{1, 0, 0}, x] - (x^3*HPL[{1, 0, 0}, x])/6 + (4*HPL[{1, 1, 0}, x])/3 + 2*x*HPL[{1, 1, 0}, x] - 4*x^2*HPL[{1, 1, 0}, x] + (2*x^3*HPL[{1, 1, 0}, x])/3 + (32*HPL[{1, 1, 1}, x])/3 + 16*x*HPL[{1, 1, 1}, x] - 32*x^2*HPL[{1, 1, 1}, x] + (16*x^3*HPL[{1, 1, 1}, x])/3 - x*HPL[{2, 0, 0}, x] + 4*x*HPL[{2, 1, 0}, x] + 32*x*HPL[{2, 1, 1}, x] - (x*HPL[{0, 0, 0, 0}, x])/2 - (5*Zeta[3])/3 - (29*x*Zeta[3])/2 + 9*x^2*Zeta[3] - (3*x^3*Zeta[3])/2 - x*HPL[{0}, x]*Zeta[3])
  + ep^4*(5971/96 - (115*Pi^2)/24 + (7*Pi^4)/180 + (4431*x)/64 - (445*Pi^2*x)/48 + (47*Pi^4*x)/120 - (5201*x^2)/32 + (505*Pi^2*x^2)/24 - (41*Pi^4*x^2)/180 + (5971*x^3)/192 - (575*Pi^2*x^3)/144 + (41*Pi^4*x^3)/1080 + (41*x*HPL[{0}, x])/12 - (19*Pi^2*x*HPL[{0}, x])/18 + (Pi^4*x*HPL[{0}, x])/180 - (2015*x^2*HPL[{0}, x])/48 + (35*Pi^2*x^2*HPL[{0}, x])/36 + (865*x^3*HPL[{0}, x])/96 - (13*Pi^2*x^3*HPL[{0}, x])/72 + (865*HPL[{1}, x])/12 - (52*Pi^2*HPL[{1}, x])/9 + (2015*x*HPL[{1}, x])/24 - (70*Pi^2*x*HPL[{1}, x])/9 - (2305*x^2*HPL[{1}, x])/12 + (148*Pi^2*x^2*HPL[{1}, x])/9 + (865*x^3*HPL[{1}, x])/24 - (26*Pi^2*x^3*HPL[{1}, x])/9 + 20*x*HPL[{2}, x] - 8*Pi^2*x*HPL[{2}, x] - (95*x^2*HPL[{2}, x])/2 + (8*Pi^2*x^2*HPL[{2}, x])/3 + (115*x^3*HPL[{2}, x])/12 - (4*Pi^2*x^3*HPL[{2}, x])/9 - (38*x*HPL[{3}, x])/3 + (8*Pi^2*x*HPL[{3}, x])/3 + (35*x^2*HPL[{3}, x])/3 - (13*x^3*HPL[{3}, x])/6 + 6*x*HPL[{4}, x] - 2*x^2*HPL[{4}, x] + (x^3*HPL[{4}, x])/3 - 2*x*HPL[{5}, x] - 5*x*HPL[{0, 0}, x] + (Pi^2*x*HPL[{0, 0}, x])/2 + (95*x^2*HPL[{0, 0}, x])/8 - (Pi^2*x^2*HPL[{0, 0}, x])/6 - (115*x^3*HPL[{0, 0}, x])/48 + (Pi^2*x^3*HPL[{0, 0}, x])/36 + (115*HPL[{1, 0}, x])/12 - (Pi^2*HPL[{1, 0}, x])/9 + (95*x*HPL[{1, 0}, x])/8 - (Pi^2*x*HPL[{1, 0}, x])/6 - (105*x^2*HPL[{1, 0}, x])/4 + (Pi^2*x^2*HPL[{1, 0}, x])/3 + (115*x^3*HPL[{1, 0}, x])/24 - (Pi^2*x^3*HPL[{1, 0}, x])/18 + (230*HPL[{1, 1}, x])/3 - (32*Pi^2*HPL[{1, 1}, x])/9 + 95*x*HPL[{1, 1}, x] - (16*Pi^2*x*HPL[{1, 1}, x])/3 - 210*x^2*HPL[{1, 1}, x] + (32*Pi^2*x^2*HPL[{1, 1}, x])/3 + (115*x^3*HPL[{1, 1}, x])/3 - (16*Pi^2*x^3*HPL[{1, 1}, x])/9 + (26*HPL[{1, 2}, x])/3 + (35*x*HPL[{1, 2}, x])/3 - (74*x^2*HPL[{1, 2}, x])/3 + (13*x^3*HPL[{1, 2}, x])/3 - (4*HPL[{1, 3}, x])/3 - 2*x*HPL[{1, 3}, x] + 4*x^2*HPL[{1, 3}, x] - (2*x^3*HPL[{1, 3}, x])/3 + (19*x*HPL[{2, 0}, x])/3 - (Pi^2*x*HPL[{2, 0}, x])/3 - (35*x^2*HPL[{2, 0}, x])/6 + (13*x^3*HPL[{2, 0}, x])/12 + (152*x*HPL[{2, 1}, x])/3 - (32*Pi^2*x*HPL[{2, 1}, x])/3 - (140*x^2*HPL[{2, 1}, x])/3 + (26*x^3*HPL[{2, 1}, x])/3 + 12*x*HPL[{2, 2}, x] - 4*x^2*HPL[{2, 2}, x] + (2*x^3*HPL[{2, 2}, x])/3 - 4*x*HPL[{2, 3}, x] - 3*x*HPL[{3, 0}, x] + x^2*HPL[{3, 0}, x] - (x^3*HPL[{3, 0}, x])/6 - 24*x*HPL[{3, 1}, x] + 8*x^2*HPL[{3, 1}, x] - (4*x^3*HPL[{3, 1}, x])/3 - 4*x*HPL[{3, 2}, x] + x*HPL[{4, 0}, x] + 8*x*HPL[{4, 1}, x] + (19*x*HPL[{0, 0, 0}, x])/6 - (Pi^2*x*HPL[{0, 0, 0}, x])/6 - (35*x^2*HPL[{0, 0, 0}, x])/12 + (13*x^3*HPL[{0, 0, 0}, x])/24 - (13*HPL[{1, 0, 0}, x])/6 - (35*x*HPL[{1, 0, 0}, x])/12 + (37*x^2*HPL[{1, 0, 0}, x])/6 - (13*x^3*HPL[{1, 0, 0}, x])/12 + (26*HPL[{1, 1, 0}, x])/3 + (35*x*HPL[{1, 1, 0}, x])/3 - (74*x^2*HPL[{1, 1, 0}, x])/3 + (13*x^3*HPL[{1, 1, 0}, x])/3 + (208*HPL[{1, 1, 1}, x])/3 + (280*x*HPL[{1, 1, 1}, x])/3 - (592*x^2*HPL[{1, 1, 1}, x])/3 + (104*x^3*HPL[{1, 1, 1}, x])/3 + (16*HPL[{1, 1, 2}, x])/3 + 8*x*HPL[{1, 1, 2}, x] - 16*x^2*HPL[{1, 1, 2}, x] + (8*x^3*HPL[{1, 1, 2}, x])/3 + (2*HPL[{1, 2, 0}, x])/3 + x*HPL[{1, 2, 0}, x] - 2*x^2*HPL[{1, 2, 0}, x] + (x^3*HPL[{1, 2, 0}, x])/3 + (16*HPL[{1, 2, 1}, x])/3 + 8*x*HPL[{1, 2, 1}, x] - 16*x^2*HPL[{1, 2, 1}, x] + (8*x^3*HPL[{1, 2, 1}, x])/3 - 3*x*HPL[{2, 0, 0}, x] + x^2*HPL[{2, 0, 0}, x] - (x^3*HPL[{2, 0, 0}, x])/6 + 12*x*HPL[{2, 1, 0}, x] - 4*x^2*HPL[{2, 1, 0}, x] + (2*x^3*HPL[{2, 1, 0}, x])/3 + 96*x*HPL[{2, 1, 1}, x] - 32*x^2*HPL[{2, 1, 1}, x] + (16*x^3*HPL[{2, 1, 1}, x])/3 + 16*x*HPL[{2, 1, 2}, x] + 2*x*HPL[{2, 2, 0}, x] + 16*x*HPL[{2, 2, 1}, x] + x*HPL[{3, 0, 0}, x] - 4*x*HPL[{3, 1, 0}, x] - 32*x*HPL[{3, 1, 1}, x] - (3*x*HPL[{0, 0, 0, 0}, x])/2 + (x^2*HPL[{0, 0, 0, 0}, x])/2 - (x^3*HPL[{0, 0, 0, 0}, x])/12 + HPL[{1, 0, 0, 0}, x]/3 + (x*HPL[{1, 0, 0, 0}, x])/2 - x^2*HPL[{1, 0, 0, 0}, x] + (x^3*HPL[{1, 0, 0, 0}, x])/6 - (4*HPL[{1, 1, 0, 0}, x])/3 - 2*x*HPL[{1, 1, 0, 0}, x] + 4*x^2*HPL[{1, 1, 0, 0}, x] - (2*x^3*HPL[{1, 1, 0, 0}, x])/3 + (16*HPL[{1, 1, 1, 0}, x])/3 + 8*x*HPL[{1, 1, 1, 0}, x] - 16*x^2*HPL[{1, 1, 1, 0}, x] + (8*x^3*HPL[{1, 1, 1, 0}, x])/3 + (128*HPL[{1, 1, 1, 1}, x])/3 + 64*x*HPL[{1, 1, 1, 1}, x] - 128*x^2*HPL[{1, 1, 1, 1}, x] + (64*x^3*HPL[{1, 1, 1, 1}, x])/3 + x*HPL[{2, 0, 0, 0}, x] - 4*x*HPL[{2, 1, 0, 0}, x] + 16*x*HPL[{2, 1, 1, 0}, x] + 128*x*HPL[{2, 1, 1, 1}, x] + (x*HPL[{0, 0, 0, 0, 0}, x])/2 - (65*Zeta[3])/6 - (479*x*Zeta[3])/12 + (14*Pi^2*x*Zeta[3])/3 + (325*x^2*Zeta[3])/6 - (39*x^3*Zeta[3])/4 - 3*x*HPL[{0}, x]*Zeta[3] + x^2*HPL[{0}, x]*Zeta[3] - (x^3*HPL[{0}, x]*Zeta[3])/6 - (28*HPL[{1}, x]*Zeta[3])/3 - 14*x*HPL[{1}, x]*Zeta[3] + 28*x^2*HPL[{1}, x]*Zeta[3] - (14*x^3*HPL[{1}, x]*Zeta[3])/3 - 28*x*HPL[{2}, x]*Zeta[3] + x*HPL[{0, 0}, x]*Zeta[3] - 24*x*Zeta[5])
  + ep^5*(39193/192 - (865*Pi^2)/48 + (91*Pi^4)/360 + (84455*x)/384 - (6701*Pi^2*x)/288 + (451*Pi^4*x)/432 - (Pi^6*x)/54 - (101017*x^2)/192 + (10945*Pi^2*x^2)/144 - (1477*Pi^4*x^2)/1080 + (39193*x^3)/384 - (4325*Pi^2*x^3)/288 + (533*Pi^4*x^3)/2160 - (133*x*HPL[{0}, x])/8 - (5*Pi^2*x*HPL[{0}, x])/3 + (Pi^4*x*HPL[{0}, x])/60 - (4431*x^2*HPL[{0}, x])/32 + (95*Pi^2*x^2*HPL[{0}, x])/24 - (Pi^4*x^2*HPL[{0}, x])/180 + (5971*x^3*HPL[{0}, x])/192 - (115*Pi^2*x^3*HPL[{0}, x])/144 + (Pi^4*x^3*HPL[{0}, x])/1080 + (5971*HPL[{1}, x])/24 - (230*Pi^2*HPL[{1}, x])/9 + (31*Pi^4*HPL[{1}, x])/135 + (4431*x*HPL[{1}, x])/16 - (95*Pi^2*x*HPL[{1}, x])/3 + (31*Pi^4*x*HPL[{1}, x])/90 - (5201*x^2*HPL[{1}, x])/8 + 70*Pi^2*x^2*HPL[{1}, x] - (31*Pi^4*x^2*HPL[{1}, x])/45 + (5971*x^3*HPL[{1}, x])/48 - (115*Pi^2*x^3*HPL[{1}, x])/9 + (31*Pi^4*x^3*HPL[{1}, x])/270 + (41*x*HPL[{2}, x])/3 - (152*Pi^2*x*HPL[{2}, x])/9 + (31*Pi^4*x*HPL[{2}, x])/45 - (2015*x^2*HPL[{2}, x])/12 + (140*Pi^2*x^2*HPL[{2}, x])/9 + (865*x^3*HPL[{2}, x])/24 - (26*Pi^2*x^3*HPL[{2}, x])/9 - 20*x*HPL[{3}, x] + 8*Pi^2*x*HPL[{3}, x] + (95*x^2*HPL[{3}, x])/2 - (8*Pi^2*x^2*HPL[{3}, x])/3 - (115*x^3*HPL[{3}, x])/12 + (4*Pi^2*x^3*HPL[{3}, x])/9 + (38*x*HPL[{4}, x])/3 - (8*Pi^2*x*HPL[{4}, x])/3 - (35*x^2*HPL[{4}, x])/3 + (13*x^3*HPL[{4}, x])/6 - 6*x*HPL[{5}, x] + 2*x^2*HPL[{5}, x] - (x^3*HPL[{5}, x])/3 + 2*x*HPL[{6}, x] - (41*x*HPL[{0, 0}, x])/12 + (19*Pi^2*x*HPL[{0, 0}, x])/18 - (Pi^4*x*HPL[{0, 0}, x])/180 + (2015*x^2*HPL[{0, 0}, x])/48 - (35*Pi^2*x^2*HPL[{0, 0}, x])/36 - (865*x^3*HPL[{0, 0}, x])/96 + (13*Pi^2*x^3*HPL[{0, 0}, x])/72 + (865*HPL[{1, 0}, x])/24 - (13*Pi^2*HPL[{1, 0}, x])/18 + (2015*x*HPL[{1, 0}, x])/48 - (35*Pi^2*x*HPL[{1, 0}, x])/36 - (2305*x^2*HPL[{1, 0}, x])/24 + (37*Pi^2*x^2*HPL[{1, 0}, x])/18 + (865*x^3*HPL[{1, 0}, x])/48 - (13*Pi^2*x^3*HPL[{1, 0}, x])/36 + (865*HPL[{1, 1}, x])/3 - (208*Pi^2*HPL[{1, 1}, x])/9 + (2015*x*HPL[{1, 1}, x])/6 - (280*Pi^2*x*HPL[{1, 1}, x])/9 - (2305*x^2*HPL[{1, 1}, x])/3 + (592*Pi^2*x^2*HPL[{1, 1}, x])/9 + (865*x^3*HPL[{1, 1}, x])/6 - (104*Pi^2*x^3*HPL[{1, 1}, x])/9 + (115*HPL[{1, 2}, x])/3 - (16*Pi^2*HPL[{1, 2}, x])/9 + (95*x*HPL[{1, 2}, x])/2 - (8*Pi^2*x*HPL[{1, 2}, x])/3 - 105*x^2*HPL[{1, 2}, x] + (16*Pi^2*x^2*HPL[{1, 2}, x])/3 + (115*x^3*HPL[{1, 2}, x])/6 - (8*Pi^2*x^3*HPL[{1, 2}, x])/9 - (26*HPL[{1, 3}, x])/3 - (35*x*HPL[{1, 3}, x])/3 + (74*x^2*HPL[{1, 3}, x])/3 - (13*x^3*HPL[{1, 3}, x])/3 + (4*HPL[{1, 4}, x])/3 + 2*x*HPL[{1, 4}, x] - 4*x^2*HPL[{1, 4}, x] + (2*x^3*HPL[{1, 4}, x])/3 + 10*x*HPL[{2, 0}, x] - Pi^2*x*HPL[{2, 0}, x] - (95*x^2*HPL[{2, 0}, x])/4 + (Pi^2*x^2*HPL[{2, 0}, x])/3 + (115*x^3*HPL[{2, 0}, x])/24 - (Pi^2*x^3*HPL[{2, 0}, x])/18 + 80*x*HPL[{2, 1}, x] - 32*Pi^2*x*HPL[{2, 1}, x] - 190*x^2*HPL[{2, 1}, x] + (32*Pi^2*x^2*HPL[{2, 1}, x])/3 + (115*x^3*HPL[{2, 1}, x])/3 - (16*Pi^2*x^3*HPL[{2, 1}, x])/9 + (76*x*HPL[{2, 2}, x])/3 - (16*Pi^2*x*HPL[{2, 2}, x])/3 - (70*x^2*HPL[{2, 2}, x])/3 + (13*x^3*HPL[{2, 2}, x])/3 - 12*x*HPL[{2, 3}, x] + 4*x^2*HPL[{2, 3}, x] - (2*x^3*HPL[{2, 3}, x])/3 + 4*x*HPL[{2, 4}, x] - (19*x*HPL[{3, 0}, x])/3 + (Pi^2*x*HPL[{3, 0}, x])/3 + (35*x^2*HPL[{3, 0}, x])/6 - (13*x^3*HPL[{3, 0}, x])/12 - (152*x*HPL[{3, 1}, x])/3 + (32*Pi^2*x*HPL[{3, 1}, x])/3 + (140*x^2*HPL[{3, 1}, x])/3 - (26*x^3*HPL[{3, 1}, x])/3 - 12*x*HPL[{3, 2}, x] + 4*x^2*HPL[{3, 2}, x] - (2*x^3*HPL[{3, 2}, x])/3 + 4*x*HPL[{3, 3}, x] + 3*x*HPL[{4, 0}, x] - x^2*HPL[{4, 0}, x] + (x^3*HPL[{4, 0}, x])/6 + 24*x*HPL[{4, 1}, x] - 8*x^2*HPL[{4, 1}, x] + (4*x^3*HPL[{4, 1}, x])/3 + 4*x*HPL[{4, 2}, x] - x*HPL[{5, 0}, x] - 8*x*HPL[{5, 1}, x] + 5*x*HPL[{0, 0, 0}, x] - (Pi^2*x*HPL[{0, 0, 0}, x])/2 - (95*x^2*HPL[{0, 0, 0}, x])/8 + (Pi^2*x^2*HPL[{0, 0, 0}, x])/6 + (115*x^3*HPL[{0, 0, 0}, x])/48 - (Pi^2*x^3*HPL[{0, 0, 0}, x])/36 - (115*HPL[{1, 0, 0}, x])/12 + (Pi^2*HPL[{1, 0, 0}, x])/9 - (95*x*HPL[{1, 0, 0}, x])/8 + (Pi^2*x*HPL[{1, 0, 0}, x])/6 + (105*x^2*HPL[{1, 0, 0}, x])/4 - (Pi^2*x^2*HPL[{1, 0, 0}, x])/3 - (115*x^3*HPL[{1, 0, 0}, x])/24 + (Pi^2*x^3*HPL[{1, 0, 0}, x])/18 + (115*HPL[{1, 1, 0}, x])/3 - (4*Pi^2*HPL[{1, 1, 0}, x])/9 + (95*x*HPL[{1, 1, 0}, x])/2 - (2*Pi^2*x*HPL[{1, 1, 0}, x])/3 - 105*x^2*HPL[{1, 1, 0}, x] + (4*Pi^2*x^2*HPL[{1, 1, 0}, x])/3 + (115*x^3*HPL[{1, 1, 0}, x])/6 - (2*Pi^2*x^3*HPL[{1, 1, 0}, x])/9 + (920*HPL[{1, 1, 1}, x])/3 - (128*Pi^2*HPL[{1, 1, 1}, x])/9 + 380*x*HPL[{1, 1, 1}, x] - (64*Pi^2*x*HPL[{1, 1, 1}, x])/3 - 840*x^2*HPL[{1, 1, 1}, x] + (128*Pi^2*x^2*HPL[{1, 1, 1}, x])/3 + (460*x^3*HPL[{1, 1, 1}, x])/3 - (64*Pi^2*x^3*HPL[{1, 1, 1}, x])/9 + (104*HPL[{1, 1, 2}, x])/3 + (140*x*HPL[{1, 1, 2}, x])/3 - (296*x^2*HPL[{1, 1, 2}, x])/3 + (52*x^3*HPL[{1, 1, 2}, x])/3 - (16*HPL[{1, 1, 3}, x])/3 - 8*x*HPL[{1, 1, 3}, x] + 16*x^2*HPL[{1, 1, 3}, x] - (8*x^3*HPL[{1, 1, 3}, x])/3 + (13*HPL[{1, 2, 0}, x])/3 + (35*x*HPL[{1, 2, 0}, x])/6 - (37*x^2*HPL[{1, 2, 0}, x])/3 + (13*x^3*HPL[{1, 2, 0}, x])/6 + (104*HPL[{1, 2, 1}, x])/3 + (140*x*HPL[{1, 2, 1}, x])/3 - (296*x^2*HPL[{1, 2, 1}, x])/3 + (52*x^3*HPL[{1, 2, 1}, x])/3 + (8*HPL[{1, 2, 2}, x])/3 + 4*x*HPL[{1, 2, 2}, x] - 8*x^2*HPL[{1, 2, 2}, x] + (4*x^3*HPL[{1, 2, 2}, x])/3 - (2*HPL[{1, 3, 0}, x])/3 - x*HPL[{1, 3, 0}, x] + 2*x^2*HPL[{1, 3, 0}, x] - (x^3*HPL[{1, 3, 0}, x])/3 - (16*HPL[{1, 3, 1}, x])/3 - 8*x*HPL[{1, 3, 1}, x] + 16*x^2*HPL[{1, 3, 1}, x] - (8*x^3*HPL[{1, 3, 1}, x])/3 - (19*x*HPL[{2, 0, 0}, x])/3 + (Pi^2*x*HPL[{2, 0, 0}, x])/3 + (35*x^2*HPL[{2, 0, 0}, x])/6 - (13*x^3*HPL[{2, 0, 0}, x])/12 + (76*x*HPL[{2, 1, 0}, x])/3 - (4*Pi^2*x*HPL[{2, 1, 0}, x])/3 - (70*x^2*HPL[{2, 1, 0}, x])/3 + (13*x^3*HPL[{2, 1, 0}, x])/3 + (608*x*HPL[{2, 1, 1}, x])/3 - (128*Pi^2*x*HPL[{2, 1, 1}, x])/3 - (560*x^2*HPL[{2, 1, 1}, x])/3 + (104*x^3*HPL[{2, 1, 1}, x])/3 + 48*x*HPL[{2, 1, 2}, x] - 16*x^2*HPL[{2, 1, 2}, x] + (8*x^3*HPL[{2, 1, 2}, x])/3 - 16*x*HPL[{2, 1, 3}, x] + 6*x*HPL[{2, 2, 0}, x] - 2*x^2*HPL[{2, 2, 0}, x] + (x^3*HPL[{2, 2, 0}, x])/3 + 48*x*HPL[{2, 2, 1}, x] - 16*x^2*HPL[{2, 2, 1}, x] + (8*x^3*HPL[{2, 2, 1}, x])/3 + 8*x*HPL[{2, 2, 2}, x] - 2*x*HPL[{2, 3, 0}, x] - 16*x*HPL[{2, 3, 1}, x] + 3*x*HPL[{3, 0, 0}, x] - x^2*HPL[{3, 0, 0}, x] + (x^3*HPL[{3, 0, 0}, x])/6 - 12*x*HPL[{3, 1, 0}, x] + 4*x^2*HPL[{3, 1, 0}, x] - (2*x^3*HPL[{3, 1, 0}, x])/3 - 96*x*HPL[{3, 1, 1}, x] + 32*x^2*HPL[{3, 1, 1}, x] - (16*x^3*HPL[{3, 1, 1}, x])/3 - 16*x*HPL[{3, 1, 2}, x] - 2*x*HPL[{3, 2, 0}, x] - 16*x*HPL[{3, 2, 1}, x] - x*HPL[{4, 0, 0}, x] + 4*x*HPL[{4, 1, 0}, x] + 32*x*HPL[{4, 1, 1}, x] - (19*x*HPL[{0, 0, 0, 0}, x])/6 + (Pi^2*x*HPL[{0, 0, 0, 0}, x])/6 + (35*x^2*HPL[{0, 0, 0, 0}, x])/12 - (13*x^3*HPL[{0, 0, 0, 0}, x])/24 + (13*HPL[{1, 0, 0, 0}, x])/6 + (35*x*HPL[{1, 0, 0, 0}, x])/12 - (37*x^2*HPL[{1, 0, 0, 0}, x])/6 + (13*x^3*HPL[{1, 0, 0, 0}, x])/12 - (26*HPL[{1, 1, 0, 0}, x])/3 - (35*x*HPL[{1, 1, 0, 0}, x])/3 + (74*x^2*HPL[{1, 1, 0, 0}, x])/3 - (13*x^3*HPL[{1, 1, 0, 0}, x])/3 + (104*HPL[{1, 1, 1, 0}, x])/3 + (140*x*HPL[{1, 1, 1, 0}, x])/3 - (296*x^2*HPL[{1, 1, 1, 0}, x])/3 + (52*x^3*HPL[{1, 1, 1, 0}, x])/3 + (832*HPL[{1, 1, 1, 1}, x])/3 + (1120*x*HPL[{1, 1, 1, 1}, x])/3 - (2368*x^2*HPL[{1, 1, 1, 1}, x])/3 + (416*x^3*HPL[{1, 1, 1, 1}, x])/3 + (64*HPL[{1, 1, 1, 2}, x])/3 + 32*x*HPL[{1, 1, 1, 2}, x] - 64*x^2*HPL[{1, 1, 1, 2}, x] + (32*x^3*HPL[{1, 1, 1, 2}, x])/3 + (8*HPL[{1, 1, 2, 0}, x])/3 + 4*x*HPL[{1, 1, 2, 0}, x] - 8*x^2*HPL[{1, 1, 2, 0}, x] + (4*x^3*HPL[{1, 1, 2, 0}, x])/3 + (64*HPL[{1, 1, 2, 1}, x])/3 + 32*x*HPL[{1, 1, 2, 1}, x] - 64*x^2*HPL[{1, 1, 2, 1}, x] + (32*x^3*HPL[{1, 1, 2, 1}, x])/3 - (2*HPL[{1, 2, 0, 0}, x])/3 - x*HPL[{1, 2, 0, 0}, x] + 2*x^2*HPL[{1, 2, 0, 0}, x] - (x^3*HPL[{1, 2, 0, 0}, x])/3 + (8*HPL[{1, 2, 1, 0}, x])/3 + 4*x*HPL[{1, 2, 1, 0}, x] - 8*x^2*HPL[{1, 2, 1, 0}, x] + (4*x^3*HPL[{1, 2, 1, 0}, x])/3 + (64*HPL[{1, 2, 1, 1}, x])/3 + 32*x*HPL[{1, 2, 1, 1}, x] - 64*x^2*HPL[{1, 2, 1, 1}, x] + (32*x^3*HPL[{1, 2, 1, 1}, x])/3 + 3*x*HPL[{2, 0, 0, 0}, x] - x^2*HPL[{2, 0, 0, 0}, x] + (x^3*HPL[{2, 0, 0, 0}, x])/6 - 12*x*HPL[{2, 1, 0, 0}, x] + 4*x^2*HPL[{2, 1, 0, 0}, x] - (2*x^3*HPL[{2, 1, 0, 0}, x])/3 + 48*x*HPL[{2, 1, 1, 0}, x] - 16*x^2*HPL[{2, 1, 1, 0}, x] + (8*x^3*HPL[{2, 1, 1, 0}, x])/3 + 384*x*HPL[{2, 1, 1, 1}, x] - 128*x^2*HPL[{2, 1, 1, 1}, x] + (64*x^3*HPL[{2, 1, 1, 1}, x])/3 + 64*x*HPL[{2, 1, 1, 2}, x] + 8*x*HPL[{2, 1, 2, 0}, x] + 64*x*HPL[{2, 1, 2, 1}, x] - 2*x*HPL[{2, 2, 0, 0}, x] + 8*x*HPL[{2, 2, 1, 0}, x] + 64*x*HPL[{2, 2, 1, 1}, x] - x*HPL[{3, 0, 0, 0}, x] + 4*x*HPL[{3, 1, 0, 0}, x] - 16*x*HPL[{3, 1, 1, 0}, x] - 128*x*HPL[{3, 1, 1, 1}, x] + (3*x*HPL[{0, 0, 0, 0, 0}, x])/2 - (x^2*HPL[{0, 0, 0, 0, 0}, x])/2 + (x^3*HPL[{0, 0, 0, 0, 0}, x])/12 - HPL[{1, 0, 0, 0, 0}, x]/3 - (x*HPL[{1, 0, 0, 0, 0}, x])/2 + x^2*HPL[{1, 0, 0, 0, 0}, x] - (x^3*HPL[{1, 0, 0, 0, 0}, x])/6 + (4*HPL[{1, 1, 0, 0, 0}, x])/3 + 2*x*HPL[{1, 1, 0, 0, 0}, x] - 4*x^2*HPL[{1, 1, 0, 0, 0}, x] + (2*x^3*HPL[{1, 1, 0, 0, 0}, x])/3 - (16*HPL[{1, 1, 1, 0, 0}, x])/3 - 8*x*HPL[{1, 1, 1, 0, 0}, x] + 16*x^2*HPL[{1, 1, 1, 0, 0}, x] - (8*x^3*HPL[{1, 1, 1, 0, 0}, x])/3 + (64*HPL[{1, 1, 1, 1, 0}, x])/3 + 32*x*HPL[{1, 1, 1, 1, 0}, x] - 64*x^2*HPL[{1, 1, 1, 1, 0}, x] + (32*x^3*HPL[{1, 1, 1, 1, 0}, x])/3 + (512*HPL[{1, 1, 1, 1, 1}, x])/3 + 256*x*HPL[{1, 1, 1, 1, 1}, x] - 512*x^2*HPL[{1, 1, 1, 1, 1}, x] + (256*x^3*HPL[{1, 1, 1, 1, 1}, x])/3 - x*HPL[{2, 0, 0, 0, 0}, x] + 4*x*HPL[{2, 1, 0, 0, 0}, x] - 16*x*HPL[{2, 1, 1, 0, 0}, x] + 64*x*HPL[{2, 1, 1, 1, 0}, x] + 512*x*HPL[{2, 1, 1, 1, 1}, x] - (x*HPL[{0, 0, 0, 0, 0, 0}, x])/2 - (575*Zeta[3])/12 + (5*Pi^2*Zeta[3])/3 - (795*x*Zeta[3])/8 + (33*Pi^2*x*Zeta[3])/2 + (905*x^2*Zeta[3])/4 - (29*Pi^2*x^2*Zeta[3])/3 - (345*x^3*Zeta[3])/8 + (29*Pi^2*x^3*Zeta[3])/18 - (19*x*HPL[{0}, x]*Zeta[3])/3 + (Pi^2*x*HPL[{0}, x]*Zeta[3])/3 + (35*x^2*HPL[{0}, x]*Zeta[3])/6 - (13*x^3*HPL[{0}, x]*Zeta[3])/12 - (182*HPL[{1}, x]*Zeta[3])/3 - (245*x*HPL[{1}, x]*Zeta[3])/3 + (518*x^2*HPL[{1}, x]*Zeta[3])/3 - (91*x^3*HPL[{1}, x]*Zeta[3])/3 - 84*x*HPL[{2}, x]*Zeta[3] + 28*x^2*HPL[{2}, x]*Zeta[3] - (14*x^3*HPL[{2}, x]*Zeta[3])/3 + 28*x*HPL[{3}, x]*Zeta[3] + 3*x*HPL[{0, 0}, x]*Zeta[3] - x^2*HPL[{0, 0}, x]*Zeta[3] + (x^3*HPL[{0, 0}, x]*Zeta[3])/6 - (2*HPL[{1, 0}, x]*Zeta[3])/3 - x*HPL[{1, 0}, x]*Zeta[3] + 2*x^2*HPL[{1, 0}, x]*Zeta[3] - (x^3*HPL[{1, 0}, x]*Zeta[3])/3 - (112*HPL[{1, 1}, x]*Zeta[3])/3 - 56*x*HPL[{1, 1}, x]*Zeta[3] + 112*x^2*HPL[{1, 1}, x]*Zeta[3] - (56*x^3*HPL[{1, 1}, x]*Zeta[3])/3 - 2*x*HPL[{2, 0}, x]*Zeta[3] - 112*x*HPL[{2, 1}, x]*Zeta[3] - x*HPL[{0, 0, 0}, x]*Zeta[3] + 24*x*Zeta[3]^2 - 9*Zeta[5] - (171*x*Zeta[5])/2 + 51*x^2*Zeta[5] - (17*x^3*Zeta[5])/2 - 3*x*HPL[{0}, x]*Zeta[5])
  + ALARM*ep^6 );
U1a[x -> 1,ep] = (1-x)^(4 - 4*ep) * (
  + 1/24 + (11*ep)/36 + ep^2*(85/54 - Pi^2/18) + ep^3*(575/81 - (11*Pi^2)/27 - (5*Zeta[3])/6) + ep^4*(7322/243 - (170*Pi^2)/81 + Pi^4/135 - (55*Zeta[3])/9) + ep^5*(90524/729 - (2300*Pi^2)/243 + (22*Pi^4)/405 - (850*Zeta[3])/27 + (10*Pi^2*Zeta[3])/9 - (17*Zeta[5])/2)
  + ALARM*ep^6 );

U2[x,ep] = (
  + (1 - x)/ep
  + 4 - 4*x + 2*HPL[{1}, x] - 2*x*HPL[{1}, x]
  + ep*(12 - Pi^2 - 12*x + Pi^2*x + 8*HPL[{1}, x] - 8*x*HPL[{1}, x] + 4*HPL[{1, 1}, x] - 4*x*HPL[{1, 1}, x])
  + ep^2*(32 - 4*Pi^2 - 32*x + 4*Pi^2*x + 24*HPL[{1}, x] - 2*Pi^2*HPL[{1}, x] - 24*x*HPL[{1}, x] + 2*Pi^2*x*HPL[{1}, x] + 16*HPL[{1, 1}, x] - 16*x*HPL[{1, 1}, x] + 8*HPL[{1, 1, 1}, x] - 8*x*HPL[{1, 1, 1}, x] - 4*Zeta[3] + 4*x*Zeta[3])
  + ep^3*(80 - 12*Pi^2 + Pi^4/3 - 80*x + 12*Pi^2*x - (Pi^4*x)/3 + 64*HPL[{1}, x] - 8*Pi^2*HPL[{1}, x] - 64*x*HPL[{1}, x] + 8*Pi^2*x*HPL[{1}, x] + 48*HPL[{1, 1}, x] - 4*Pi^2*HPL[{1, 1}, x] - 48*x*HPL[{1, 1}, x] + 4*Pi^2*x*HPL[{1, 1}, x] + 32*HPL[{1, 1, 1}, x] - 32*x*HPL[{1, 1, 1}, x] + 16*HPL[{1, 1, 1, 1}, x] - 16*x*HPL[{1, 1, 1, 1}, x] - 16*Zeta[3] + 16*x*Zeta[3] - 8*HPL[{1}, x]*Zeta[3] + 8*x*HPL[{1}, x]*Zeta[3])
  + ep^4*(192 - 32*Pi^2 + (4*Pi^4)/3 - 192*x + 32*Pi^2*x - (4*Pi^4*x)/3 + 160*HPL[{1}, x] - 24*Pi^2*HPL[{1}, x] + (2*Pi^4*HPL[{1}, x])/3 - 160*x*HPL[{1}, x] + 24*Pi^2*x*HPL[{1}, x] - (2*Pi^4*x*HPL[{1}, x])/3 + 128*HPL[{1, 1}, x] - 16*Pi^2*HPL[{1, 1}, x] - 128*x*HPL[{1, 1}, x] + 16*Pi^2*x*HPL[{1, 1}, x] + 96*HPL[{1, 1, 1}, x] - 8*Pi^2*HPL[{1, 1, 1}, x] - 96*x*HPL[{1, 1, 1}, x] + 8*Pi^2*x*HPL[{1, 1, 1}, x] + 64*HPL[{1, 1, 1, 1}, x] - 64*x*HPL[{1, 1, 1, 1}, x] + 32*HPL[{1, 1, 1, 1, 1}, x] - 32*x*HPL[{1, 1, 1, 1, 1}, x] - 48*Zeta[3] + 4*Pi^2*Zeta[3] + 48*x*Zeta[3] - 4*Pi^2*x*Zeta[3] - 32*HPL[{1}, x]*Zeta[3] + 32*x*HPL[{1}, x]*Zeta[3] - 16*HPL[{1, 1}, x]*Zeta[3] + 16*x*HPL[{1, 1}, x]*Zeta[3] - 12*Zeta[5] + 12*x*Zeta[5])
  + ALARM*ep^5 );
U2[x -> 1,ep] = (1-x)^(1 - 2*ep) * (
  + ep^(-1) + 4 + ep*(12 - Pi^2) + ep^2*(32 - 4*Pi^2 - 4*Zeta[3]) + ep^3*(80 - 12*Pi^2 + Pi^4/3 - 16*Zeta[3]) + ep^4*(192 - 12*Zeta[5] - 48*Zeta[3] - 32*Pi^2 + 4*Pi^2*Zeta[3] + 4/3*Pi^4)
  + ALARM*ep^5 );

U3[x,ep] = (
  + (1 - x)/ep
  + 4 - 4*x - HPL[{0}, x] + x*HPL[{0}, x] + 2*HPL[{1}, x] - 2*x*HPL[{1}, x]
  + ep*(12 - Pi^2 - 12*x + Pi^2*x - 4*HPL[{0}, x] + 4*x*HPL[{0}, x] + 8*HPL[{1}, x] - 8*x*HPL[{1}, x] - 2*HPL[{2}, x] + 2*x*HPL[{2}, x] + HPL[{0, 0}, x] - x*HPL[{0, 0}, x] - 2*HPL[{1, 0}, x] + 2*x*HPL[{1, 0}, x] + 4*HPL[{1, 1}, x] - 4*x*HPL[{1, 1}, x])
  + ep^2*(32 - 4*Pi^2 - 32*x + 4*Pi^2*x - 12*HPL[{0}, x] + Pi^2*HPL[{0}, x] + 12*x*HPL[{0}, x] - Pi^2*x*HPL[{0}, x] + 24*HPL[{1}, x] - 2*Pi^2*HPL[{1}, x] - 24*x*HPL[{1}, x] + 2*Pi^2*x*HPL[{1}, x] - 8*HPL[{2}, x] + 8*x*HPL[{2}, x] + 2*HPL[{3}, x] - 2*x*HPL[{3}, x] + 4*HPL[{0, 0}, x] - 4*x*HPL[{0, 0}, x] - 8*HPL[{1, 0}, x] + 8*x*HPL[{1, 0}, x] + 16*HPL[{1, 1}, x] - 16*x*HPL[{1, 1}, x] - 4*HPL[{1, 2}, x] + 4*x*HPL[{1, 2}, x] + 2*HPL[{2, 0}, x] - 2*x*HPL[{2, 0}, x] - 4*HPL[{2, 1}, x] + 4*x*HPL[{2, 1}, x] - HPL[{0, 0, 0}, x] + x*HPL[{0, 0, 0}, x] + 2*HPL[{1, 0, 0}, x] - 2*x*HPL[{1, 0, 0}, x] - 4*HPL[{1, 1, 0}, x] + 4*x*HPL[{1, 1, 0}, x] + 8*HPL[{1, 1, 1}, x] - 8*x*HPL[{1, 1, 1}, x] - 4*Zeta[3] + 4*x*Zeta[3])
  + ep^3*(80 - 12*Pi^2 + Pi^4/3 - 80*x + 12*Pi^2*x - (Pi^4*x)/3 - 32*HPL[{0}, x] + 4*Pi^2*HPL[{0}, x] + 32*x*HPL[{0}, x] - 4*Pi^2*x*HPL[{0}, x] + 64*HPL[{1}, x] - 8*Pi^2*HPL[{1}, x] - 64*x*HPL[{1}, x] + 8*Pi^2*x*HPL[{1}, x] - 24*HPL[{2}, x] + 2*Pi^2*HPL[{2}, x] + 24*x*HPL[{2}, x] - 2*Pi^2*x*HPL[{2}, x] + 8*HPL[{3}, x] - 8*x*HPL[{3}, x] - 2*HPL[{4}, x] + 2*x*HPL[{4}, x] + 12*HPL[{0, 0}, x] - Pi^2*HPL[{0, 0}, x] - 12*x*HPL[{0, 0}, x] + Pi^2*x*HPL[{0, 0}, x] - 24*HPL[{1, 0}, x] + 2*Pi^2*HPL[{1, 0}, x] + 24*x*HPL[{1, 0}, x] - 2*Pi^2*x*HPL[{1, 0}, x] + 48*HPL[{1, 1}, x] - 4*Pi^2*HPL[{1, 1}, x] - 48*x*HPL[{1, 1}, x] + 4*Pi^2*x*HPL[{1, 1}, x] - 16*HPL[{1, 2}, x] + 16*x*HPL[{1, 2}, x] + 4*HPL[{1, 3}, x] - 4*x*HPL[{1, 3}, x] + 8*HPL[{2, 0}, x] - 8*x*HPL[{2, 0}, x] - 16*HPL[{2, 1}, x] + 16*x*HPL[{2, 1}, x] + 4*HPL[{2, 2}, x] - 4*x*HPL[{2, 2}, x] - 2*HPL[{3, 0}, x] + 2*x*HPL[{3, 0}, x] + 4*HPL[{3, 1}, x] - 4*x*HPL[{3, 1}, x] - 4*HPL[{0, 0, 0}, x] + 4*x*HPL[{0, 0, 0}, x] + 8*HPL[{1, 0, 0}, x] - 8*x*HPL[{1, 0, 0}, x] - 16*HPL[{1, 1, 0}, x] + 16*x*HPL[{1, 1, 0}, x] + 32*HPL[{1, 1, 1}, x] - 32*x*HPL[{1, 1, 1}, x] - 8*HPL[{1, 1, 2}, x] + 8*x*HPL[{1, 1, 2}, x] + 4*HPL[{1, 2, 0}, x] - 4*x*HPL[{1, 2, 0}, x] - 8*HPL[{1, 2, 1}, x] + 8*x*HPL[{1, 2, 1}, x] - 2*HPL[{2, 0, 0}, x] + 2*x*HPL[{2, 0, 0}, x] + 4*HPL[{2, 1, 0}, x] - 4*x*HPL[{2, 1, 0}, x] - 8*HPL[{2, 1, 1}, x] + 8*x*HPL[{2, 1, 1}, x] + HPL[{0, 0, 0, 0}, x] - x*HPL[{0, 0, 0, 0}, x] - 2*HPL[{1, 0, 0, 0}, x] + 2*x*HPL[{1, 0, 0, 0}, x] + 4*HPL[{1, 1, 0, 0}, x] - 4*x*HPL[{1, 1, 0, 0}, x] - 8*HPL[{1, 1, 1, 0}, x] + 8*x*HPL[{1, 1, 1, 0}, x] + 16*HPL[{1, 1, 1, 1}, x] - 16*x*HPL[{1, 1, 1, 1}, x] - 16*Zeta[3] + 16*x*Zeta[3] + 4*HPL[{0}, x]*Zeta[3] - 4*x*HPL[{0}, x]*Zeta[3] - 8*HPL[{1}, x]*Zeta[3] + 8*x*HPL[{1}, x]*Zeta[3])
  + ep^4*(192 - 32*Pi^2 + (4*Pi^4)/3 - 192*x + 32*Pi^2*x - (4*Pi^4*x)/3 - 80*HPL[{0}, x] + 12*Pi^2*HPL[{0}, x] - (Pi^4*HPL[{0}, x])/3 + 80*x*HPL[{0}, x] - 12*Pi^2*x*HPL[{0}, x] + (Pi^4*x*HPL[{0}, x])/3 + 160*HPL[{1}, x] - 24*Pi^2*HPL[{1}, x] + (2*Pi^4*HPL[{1}, x])/3 - 160*x*HPL[{1}, x] + 24*Pi^2*x*HPL[{1}, x] - (2*Pi^4*x*HPL[{1}, x])/3 - 64*HPL[{2}, x] + 8*Pi^2*HPL[{2}, x] + 64*x*HPL[{2}, x] - 8*Pi^2*x*HPL[{2}, x] + 24*HPL[{3}, x] - 2*Pi^2*HPL[{3}, x] - 24*x*HPL[{3}, x] + 2*Pi^2*x*HPL[{3}, x] - 8*HPL[{4}, x] + 8*x*HPL[{4}, x] + 2*HPL[{5}, x] - 2*x*HPL[{5}, x] + 32*HPL[{0, 0}, x] - 4*Pi^2*HPL[{0, 0}, x] - 32*x*HPL[{0, 0}, x] + 4*Pi^2*x*HPL[{0, 0}, x] - 64*HPL[{1, 0}, x] + 8*Pi^2*HPL[{1, 0}, x] + 64*x*HPL[{1, 0}, x] - 8*Pi^2*x*HPL[{1, 0}, x] + 128*HPL[{1, 1}, x] - 16*Pi^2*HPL[{1, 1}, x] - 128*x*HPL[{1, 1}, x] + 16*Pi^2*x*HPL[{1, 1}, x] - 48*HPL[{1, 2}, x] + 4*Pi^2*HPL[{1, 2}, x] + 48*x*HPL[{1, 2}, x] - 4*Pi^2*x*HPL[{1, 2}, x] + 16*HPL[{1, 3}, x] - 16*x*HPL[{1, 3}, x] - 4*HPL[{1, 4}, x] + 4*x*HPL[{1, 4}, x] + 24*HPL[{2, 0}, x] - 2*Pi^2*HPL[{2, 0}, x] - 24*x*HPL[{2, 0}, x] + 2*Pi^2*x*HPL[{2, 0}, x] - 48*HPL[{2, 1}, x] + 4*Pi^2*HPL[{2, 1}, x] + 48*x*HPL[{2, 1}, x] - 4*Pi^2*x*HPL[{2, 1}, x] + 16*HPL[{2, 2}, x] - 16*x*HPL[{2, 2}, x] - 4*HPL[{2, 3}, x] + 4*x*HPL[{2, 3}, x] - 8*HPL[{3, 0}, x] + 8*x*HPL[{3, 0}, x] + 16*HPL[{3, 1}, x] - 16*x*HPL[{3, 1}, x] - 4*HPL[{3, 2}, x] + 4*x*HPL[{3, 2}, x] + 2*HPL[{4, 0}, x] - 2*x*HPL[{4, 0}, x] - 4*HPL[{4, 1}, x] + 4*x*HPL[{4, 1}, x] - 12*HPL[{0, 0, 0}, x] + Pi^2*HPL[{0, 0, 0}, x] + 12*x*HPL[{0, 0, 0}, x] - Pi^2*x*HPL[{0, 0, 0}, x] + 24*HPL[{1, 0, 0}, x] - 2*Pi^2*HPL[{1, 0, 0}, x] - 24*x*HPL[{1, 0, 0}, x] + 2*Pi^2*x*HPL[{1, 0, 0}, x] - 48*HPL[{1, 1, 0}, x] + 4*Pi^2*HPL[{1, 1, 0}, x] + 48*x*HPL[{1, 1, 0}, x] - 4*Pi^2*x*HPL[{1, 1, 0}, x] + 96*HPL[{1, 1, 1}, x] - 8*Pi^2*HPL[{1, 1, 1}, x] - 96*x*HPL[{1, 1, 1}, x] + 8*Pi^2*x*HPL[{1, 1, 1}, x] - 32*HPL[{1, 1, 2}, x] + 32*x*HPL[{1, 1, 2}, x] + 8*HPL[{1, 1, 3}, x] - 8*x*HPL[{1, 1, 3}, x] + 16*HPL[{1, 2, 0}, x] - 16*x*HPL[{1, 2, 0}, x] - 32*HPL[{1, 2, 1}, x] + 32*x*HPL[{1, 2, 1}, x] + 8*HPL[{1, 2, 2}, x] - 8*x*HPL[{1, 2, 2}, x] - 4*HPL[{1, 3, 0}, x] + 4*x*HPL[{1, 3, 0}, x] + 8*HPL[{1, 3, 1}, x] - 8*x*HPL[{1, 3, 1}, x] - 8*HPL[{2, 0, 0}, x] + 8*x*HPL[{2, 0, 0}, x] + 16*HPL[{2, 1, 0}, x] - 16*x*HPL[{2, 1, 0}, x] - 32*HPL[{2, 1, 1}, x] + 32*x*HPL[{2, 1, 1}, x] + 8*HPL[{2, 1, 2}, x] - 8*x*HPL[{2, 1, 2}, x] - 4*HPL[{2, 2, 0}, x] + 4*x*HPL[{2, 2, 0}, x] + 8*HPL[{2, 2, 1}, x] - 8*x*HPL[{2, 2, 1}, x] + 2*HPL[{3, 0, 0}, x] - 2*x*HPL[{3, 0, 0}, x] - 4*HPL[{3, 1, 0}, x] + 4*x*HPL[{3, 1, 0}, x] + 8*HPL[{3, 1, 1}, x] - 8*x*HPL[{3, 1, 1}, x] + 4*HPL[{0, 0, 0, 0}, x] - 4*x*HPL[{0, 0, 0, 0}, x] - 8*HPL[{1, 0, 0, 0}, x] + 8*x*HPL[{1, 0, 0, 0}, x] + 16*HPL[{1, 1, 0, 0}, x] - 16*x*HPL[{1, 1, 0, 0}, x] - 32*HPL[{1, 1, 1, 0}, x] + 32*x*HPL[{1, 1, 1, 0}, x] + 64*HPL[{1, 1, 1, 1}, x] - 64*x*HPL[{1, 1, 1, 1}, x] - 16*HPL[{1, 1, 1, 2}, x] + 16*x*HPL[{1, 1, 1, 2}, x] + 8*HPL[{1, 1, 2, 0}, x] - 8*x*HPL[{1, 1, 2, 0}, x] - 16*HPL[{1, 1, 2, 1}, x] + 16*x*HPL[{1, 1, 2, 1}, x] - 4*HPL[{1, 2, 0, 0}, x] + 4*x*HPL[{1, 2, 0, 0}, x] + 8*HPL[{1, 2, 1, 0}, x] - 8*x*HPL[{1, 2, 1, 0}, x] - 16*HPL[{1, 2, 1, 1}, x] + 16*x*HPL[{1, 2, 1, 1}, x] + 2*HPL[{2, 0, 0, 0}, x] - 2*x*HPL[{2, 0, 0, 0}, x] - 4*HPL[{2, 1, 0, 0}, x] + 4*x*HPL[{2, 1, 0, 0}, x] + 8*HPL[{2, 1, 1, 0}, x] - 8*x*HPL[{2, 1, 1, 0}, x] - 16*HPL[{2, 1, 1, 1}, x] + 16*x*HPL[{2, 1, 1, 1}, x] - HPL[{0, 0, 0, 0, 0}, x] + x*HPL[{0, 0, 0, 0, 0}, x] + 2*HPL[{1, 0, 0, 0, 0}, x] - 2*x*HPL[{1, 0, 0, 0, 0}, x] - 4*HPL[{1, 1, 0, 0, 0}, x] + 4*x*HPL[{1, 1, 0, 0, 0}, x] + 8*HPL[{1, 1, 1, 0, 0}, x] - 8*x*HPL[{1, 1, 1, 0, 0}, x] - 16*HPL[{1, 1, 1, 1, 0}, x] + 16*x*HPL[{1, 1, 1, 1, 0}, x] + 32*HPL[{1, 1, 1, 1, 1}, x] - 32*x*HPL[{1, 1, 1, 1, 1}, x] - 48*Zeta[3] + 4*Pi^2*Zeta[3] + 48*x*Zeta[3] - 4*Pi^2*x*Zeta[3] + 16*HPL[{0}, x]*Zeta[3] - 16*x*HPL[{0}, x]*Zeta[3] - 32*HPL[{1}, x]*Zeta[3] + 32*x*HPL[{1}, x]*Zeta[3] + 8*HPL[{2}, x]*Zeta[3] - 8*x*HPL[{2}, x]*Zeta[3] - 4*HPL[{0, 0}, x]*Zeta[3] + 4*x*HPL[{0, 0}, x]*Zeta[3] + 8*HPL[{1, 0}, x]*Zeta[3] - 8*x*HPL[{1, 0}, x]*Zeta[3] - 16*HPL[{1, 1}, x]*Zeta[3] + 16*x*HPL[{1, 1}, x]*Zeta[3] - 12*Zeta[5] + 12*x*Zeta[5])
  + ALARM*ep^5 );
U3[x -> 1,ep] = (1-x)^(1 - 2*ep) * (
  + ep^(-1) + 4 + ep*(12 - Pi^2) + ep^2*(32 - 4*Pi^2 - 4*Zeta[3]) + ep^3*(80 - 12*Pi^2 + Pi^4/3 - 16*Zeta[3]) + ep^4 * (192 - 12*Zeta[5] - 48*Zeta[3] - 32*Pi^2 + 4*Pi^2*Zeta[3] + 4/3*Pi^4)
  + ALARM*ep^5 );

U4[x,ep] = (
  + (1 - x)/ep
  + 5 - 5*x + 3*HPL[{1}, x] - 3*x*HPL[{1}, x]
  + ep*(19 - (2*Pi^2)/3 - 19*x + (2*Pi^2*x)/3 + 15*HPL[{1}, x] - 15*x*HPL[{1}, x] + 9*HPL[{1, 1}, x] - 9*x*HPL[{1, 1}, x])
  + ep^2*(65 - (10*Pi^2)/3 - 65*x + (10*Pi^2*x)/3 + 57*HPL[{1}, x] - 2*Pi^2*HPL[{1}, x] - 57*x*HPL[{1}, x] + 2*Pi^2*x*HPL[{1}, x] + 45*HPL[{1, 1}, x] - 45*x*HPL[{1, 1}, x] + 27*HPL[{1, 1, 1}, x] - 27*x*HPL[{1, 1, 1}, x] - 8*Zeta[3] + 8*x*Zeta[3])
  + ep^3*(211 - (38*Pi^2)/3 - 211*x + (38*Pi^2*x)/3 + 195*HPL[{1}, x] - 10*Pi^2*HPL[{1}, x] - 195*x*HPL[{1}, x] + 10*Pi^2*x*HPL[{1}, x] + 171*HPL[{1, 1}, x] - 6*Pi^2*HPL[{1, 1}, x] - 171*x*HPL[{1, 1}, x] + 6*Pi^2*x*HPL[{1, 1}, x] + 135*HPL[{1, 1, 1}, x] - 135*x*HPL[{1, 1, 1}, x] + 81*HPL[{1, 1, 1, 1}, x] - 81*x*HPL[{1, 1, 1, 1}, x] - 40*Zeta[3] + 40*x*Zeta[3] - 24*HPL[{1}, x]*Zeta[3] + 24*x*HPL[{1}, x]*Zeta[3])
  + ep^4*(665 - (130*Pi^2)/3 - 665*x + (130*Pi^2*x)/3 + 633*HPL[{1}, x] - 38*Pi^2*HPL[{1}, x] - 633*x*HPL[{1}, x] + 38*Pi^2*x*HPL[{1}, x] + 585*HPL[{1, 1}, x] - 30*Pi^2*HPL[{1, 1}, x] - 585*x*HPL[{1, 1}, x] + 30*Pi^2*x*HPL[{1, 1}, x] + 513*HPL[{1, 1, 1}, x] - 18*Pi^2*HPL[{1, 1, 1}, x] - 513*x*HPL[{1, 1, 1}, x] + 18*Pi^2*x*HPL[{1, 1, 1}, x] + 405*HPL[{1, 1, 1, 1}, x] - 405*x*HPL[{1, 1, 1, 1}, x] + 243*HPL[{1, 1, 1, 1, 1}, x] - 243*x*HPL[{1, 1, 1, 1, 1}, x] - 152*Zeta[3] + (16*Pi^2*Zeta[3])/3 + 152*x*Zeta[3] - (16*Pi^2*x*Zeta[3])/3 - 120*HPL[{1}, x]*Zeta[3] + 120*x*HPL[{1}, x]*Zeta[3] - 72*HPL[{1, 1}, x]*Zeta[3] + 72*x*HPL[{1, 1}, x]*Zeta[3] - 48*Zeta[5] + 48*x*Zeta[5])
  + ALARM*ep^5 );
U4[x -> 1,ep] = (1-x)^(1 - 3*ep) * (
  + ep^(-1) + 5 + ep*(19 - (2*Pi^2)/3) + ep^2*(65 - (10*Pi^2)/3 - 8*Zeta[3]) + ep^3*(211 - (38*Pi^2)/3 - 40*Zeta[3]) + ep^4 * (665 - 48*Zeta[5] - 152*Zeta[3] - 130/3*Pi^2 + 16/3*Pi^2*Zeta[3])
  + ALARM*ep^5 );
  
U5[x,ep] = (
  + (Pi^2/6 - HPL[{2}, x] - HPL[{0, 0}, x])/ep
  - (2*Pi^2*HPL[{0}, x])/3 - 2*HPL[{2, 0}, x] - 5*HPL[{2, 1}, x] + 3*HPL[{0, 0, 0}, x] + Zeta[3]
  + ep*((17*Pi^4)/60 + 2*HPL[{4}, x] + (7*Pi^2*HPL[{0, 0}, x])/3 - 4*HPL[{2, 2}, x] + 6*HPL[{3, 0}, x] + 6*HPL[{3, 1}, x] + 2*HPL[{2, 0, 0}, x] - 4*HPL[{2, 1, 0}, x] - 19*HPL[{2, 1, 1}, x] - 7*HPL[{0, 0, 0, 0}, x] + 2*HPL[{0}, x]*Zeta[3])
  + ep^2*((-7*Pi^4*HPL[{0}, x])/30 + 2*Pi^2*HPL[{3}, x] - 6*HPL[{5}, x] + 2*Pi^2*HPL[{2, 0}, x] + 2*Pi^2*HPL[{2, 1}, x] + 4*HPL[{2, 3}, x] + 12*HPL[{3, 2}, x] - 14*HPL[{4, 0}, x] - 8*HPL[{4, 1}, x] - (17*Pi^2*HPL[{0, 0, 0}, x])/3 - 8*HPL[{2, 1, 2}, x] + 4*HPL[{2, 2, 0}, x] - 8*HPL[{2, 2, 1}, x] - 6*HPL[{3, 0, 0}, x] + 12*HPL[{3, 1, 0}, x] + 30*HPL[{3, 1, 1}, x] - 2*HPL[{2, 0, 0, 0}, x] + 4*HPL[{2, 1, 0, 0}, x] - 8*HPL[{2, 1, 1, 0}, x] - 65*HPL[{2, 1, 1, 1}, x] + 15*HPL[{0, 0, 0, 0, 0}, x] - (Pi^2*Zeta[3])/3 + 16*HPL[{2}, x]*Zeta[3] + 15*Zeta[5])
  + ep^3*((-47*Pi^6)/3780 + (2*Pi^4*HPL[{2}, x])/3 - 6*Pi^2*HPL[{4}, x] + 14*HPL[{6}, x] + (2*Pi^4*HPL[{0, 0}, x])/15 + 4*Pi^2*HPL[{2, 2}, x] - 4*HPL[{2, 4}, x] - 6*Pi^2*HPL[{3, 0}, x] - 12*HPL[{3, 3}, x] - 28*HPL[{4, 2}, x] + 30*HPL[{5, 0}, x] + 12*HPL[{5, 1}, x] - 2*Pi^2*HPL[{2, 0, 0}, x] + 4*Pi^2*HPL[{2, 1, 0}, x] + 10*Pi^2*HPL[{2, 1, 1}, x] + 8*HPL[{2, 1, 3}, x] + 8*HPL[{2, 2, 2}, x] - 4*HPL[{2, 3, 0}, x] + 8*HPL[{2, 3, 1}, x] + 24*HPL[{3, 1, 2}, x] - 12*HPL[{3, 2, 0}, x] + 24*HPL[{3, 2, 1}, x] + 14*HPL[{4, 0, 0}, x] - 28*HPL[{4, 1, 0}, x] - 52*HPL[{4, 1, 1}, x] + (37*Pi^2*HPL[{0, 0, 0, 0}, x])/3 - 16*HPL[{2, 1, 1, 2}, x] + 8*HPL[{2, 1, 2, 0}, x] - 16*HPL[{2, 1, 2, 1}, x] - 4*HPL[{2, 2, 0, 0}, x] + 8*HPL[{2, 2, 1, 0}, x] - 16*HPL[{2, 2, 1, 1}, x] + 6*HPL[{3, 0, 0, 0}, x] - 12*HPL[{3, 1, 0, 0}, x] + 24*HPL[{3, 1, 1, 0}, x] + 114*HPL[{3, 1, 1, 1}, x] + 2*HPL[{2, 0, 0, 0, 0}, x] - 4*HPL[{2, 1, 0, 0, 0}, x] + 8*HPL[{2, 1, 1, 0, 0}, x] - 16*HPL[{2, 1, 1, 1, 0}, x] - 211*HPL[{2, 1, 1, 1, 1}, x] - 31*HPL[{0, 0, 0, 0, 0, 0}, x] - (2*Pi^2*HPL[{0}, x]*Zeta[3])/3 - 24*HPL[{3}, x]*Zeta[3] + 8*HPL[{2, 0}, x]*Zeta[3] + 56*HPL[{2, 1}, x]*Zeta[3] - 4*HPL[{0, 0, 0}, x]*Zeta[3] - 9*Zeta[3]^2 + 6*HPL[{0}, x]*Zeta[5])
  + ALARM*ep^4 );
U5[x -> 1,ep] = (
  + (1-x)^(1 - 3*ep) * (
    + ep^(-2) + 3/ep + 9 - (2*Pi^2)/3 + ep*(27 - 2*Pi^2 - 8*Zeta[3]) + ep^2*(81 - 6*Pi^2 - 24*Zeta[3])
    + ALARM*ep^3 )
  + (1-x)^(1 - 2*ep) * (
    - ep^(-2) - 2/ep - 4 + Pi^2 + 2*ep*(-4 + Pi^2 + 2*Zeta[3]) + ep^2*(4*Pi^2 - Pi^4/3 + 8*(-2 + Zeta[3]))
    + ALARM*ep^3 ));

U6[x,ep] = (
  + HPL[{0, 0}, x]/ep
  - (2*Pi^2*HPL[{0}, x])/3 + 4*HPL[{3}, x] - 3*HPL[{0, 0, 0}, x]
  + ep*((7*Pi^4)/45 - 12*HPL[{4}, x] + Pi^2*HPL[{0, 0}, x] + 2*HPL[{3, 0}, x] + 16*HPL[{3, 1}, x] + 7*HPL[{0, 0, 0, 0}, x]) - 4*Zeta[3]
  + ep^2*((-4*Pi^4*HPL[{0}, x])/45 - (16*Pi^2*HPL[{3}, x])/3 + 28*HPL[{5}, x] + 8*HPL[{3, 2}, x] - 6*HPL[{4, 0}, x] - 48*HPL[{4, 1}, x] - (5*Pi^2*HPL[{0, 0, 0}, x])/3 - 2*HPL[{3, 0, 0}, x] + 8*HPL[{3, 1, 0}, x] + 64*HPL[{3, 1, 1}, x] - 15*HPL[{0, 0, 0, 0, 0}, x] + 4*Pi^2*Zeta[3] - 2*HPL[{0, 0}, x]*Zeta[3] - 24*Zeta[5])
  + ep^3*(-Pi^6/945 + 16*Pi^2*HPL[{4}, x] - 60*HPL[{6}, x] + (17*Pi^4*HPL[{0, 0}, x])/90 - (2*Pi^2*HPL[{3, 0}, x])/3 - (64*Pi^2*HPL[{3, 1}, x])/3 - 8*HPL[{3, 3}, x] - 24*HPL[{4, 2}, x] + 14*HPL[{5, 0}, x] + 112*HPL[{5, 1}, x] + 32*HPL[{3, 1, 2}, x] + 4*HPL[{3, 2, 0}, x] + 32*HPL[{3, 2, 1}, x] + 6*HPL[{4, 0, 0}, x] - 24*HPL[{4, 1, 0}, x] - 192*HPL[{4, 1, 1}, x] + 3*Pi^2*HPL[{0, 0, 0, 0}, x] + 2*HPL[{3, 0, 0, 0}, x] - 8*HPL[{3, 1, 0, 0}, x] + 32*HPL[{3, 1, 1, 0}, x] + 256*HPL[{3, 1, 1, 1}, x] + 31*HPL[{0, 0, 0, 0, 0, 0}, x] + (4*Pi^2*HPL[{0}, x]*Zeta[3])/3 - 56*HPL[{3}, x]*Zeta[3] + 6*HPL[{0, 0, 0}, x]*Zeta[3] + 24*Zeta[3]^2)
  + ALARM*ep^4 );
U6[x -> 1,ep] = (1-x)^(2 - 4*ep) * (
  + 1/(2*ep) + 3 + ep*(14 - (2*Pi^2)/3) + ep^2*(60 - 4*Pi^2 - 10*Zeta[3]) + ep^3*(248 - (56*Pi^2)/3 + (4*Pi^4)/45 - 60*Zeta[3])
  + ALARM*ep^4 );

U7[x,ep] = (
  + (Pi^2*HPL[{0}, x])/6 + HPL[{2, 0}, x] + HPL[{0, 0, 0}, x] + 2*Zeta[3]
  + ep*((-23*Pi^4)/180 - (Pi^2*HPL[{2}, x])/3 + 4*HPL[{4}, x] - (7*Pi^2*HPL[{0, 0}, x])/6 + 4*HPL[{2, 2}, x] - 3*HPL[{3, 0}, x] + HPL[{2, 0, 0}, x] + 2*HPL[{2, 1, 0}, x] - 4*HPL[{0, 0, 0, 0}, x] - 5*HPL[{0}, x]*Zeta[3])
  + ep^2*((7*Pi^4*HPL[{0}, x])/18 + Pi^2*HPL[{3}, x] - 16*HPL[{5}, x] - 2*Pi^2*HPL[{2, 0}, x] - (2*Pi^2*HPL[{2, 1}, x])/3 + 4*HPL[{2, 3}, x] - 12*HPL[{3, 2}, x] + 9*HPL[{4, 0}, x] + 16*HPL[{4, 1}, x] + (17*Pi^2*HPL[{0, 0, 0}, x])/6 + 8*HPL[{2, 1, 2}, x] + 16*HPL[{2, 2, 1}, x] - 3*HPL[{3, 0, 0}, x] - 6*HPL[{3, 1, 0}, x] - 3*HPL[{2, 0, 0, 0}, x] + 2*HPL[{2, 1, 0, 0}, x] + 4*HPL[{2, 1, 1, 0}, x] + 11*HPL[{0, 0, 0, 0, 0}, x] + Pi^2*Zeta[3] - 10*HPL[{2}, x]*Zeta[3] + 3*HPL[{0, 0}, x]*Zeta[3] + 4*Zeta[5])
  + ep^3*((-19*Pi^6)/270 + (22*Pi^4*HPL[{2}, x])/45 - (23*Pi^2*HPL[{4}, x])/3 + 44*HPL[{6}, x] - (31*Pi^4*HPL[{0, 0}, x])/45 - (14*Pi^2*HPL[{2, 2}, x])/3 - 12*HPL[{2, 4}, x] + 6*Pi^2*HPL[{3, 0}, x] + 2*Pi^2*HPL[{3, 1}, x] - 12*HPL[{3, 3}, x] + 36*HPL[{4, 2}, x] - 23*HPL[{5, 0}, x] - 64*HPL[{5, 1}, x] + (4*Pi^2*HPL[{2, 0, 0}, x])/3 - 4*Pi^2*HPL[{2, 1, 0}, x] - (4*Pi^2*HPL[{2, 1, 1}, x])/3 + 8*HPL[{2, 1, 3}, x] + 4*HPL[{2, 3, 0}, x] + 16*HPL[{2, 3, 1}, x] - 24*HPL[{3, 1, 2}, x] - 48*HPL[{3, 2, 1}, x] + 5*HPL[{4, 0, 0}, x] + 22*HPL[{4, 1, 0}, x] + 64*HPL[{4, 1, 1}, x] - (35*Pi^2*HPL[{0, 0, 0, 0}, x])/6 + 16*HPL[{2, 1, 1, 2}, x] + 32*HPL[{2, 1, 2, 1}, x] - 4*HPL[{2, 2, 0, 0}, x] + 4*HPL[{2, 2, 1, 0}, x] + 64*HPL[{2, 2, 1, 1}, x] + 9*HPL[{3, 0, 0, 0}, x] - 6*HPL[{3, 1, 0, 0}, x] - 12*HPL[{3, 1, 1, 0}, x] + 5*HPL[{2, 0, 0, 0, 0}, x] - 6*HPL[{2, 1, 0, 0, 0}, x] + 4*HPL[{2, 1, 1, 0, 0}, x] + 8*HPL[{2, 1, 1, 1, 0}, x] - 26*HPL[{0, 0, 0, 0, 0, 0}, x] + (10*Pi^2*HPL[{0}, x]*Zeta[3])/3 + 30*HPL[{3}, x]*Zeta[3] - 16*HPL[{2, 0}, x]*Zeta[3] - 20*HPL[{2, 1}, x]*Zeta[3] - HPL[{0, 0, 0}, x]*Zeta[3] + 11*Zeta[3]^2 - 25*HPL[{0}, x]*Zeta[5])
  + ALARM*ep^4 );
U7[x -> 1,ep] = (1-x)^(2 - 4*ep) * (
  - 1/2 - 4*ep + ep^2*(-22 + (2*Pi^2)/3) + ep^3*(-104 + (16*Pi^2)/3 + 10*Zeta[3])
  + ALARM*ep^4 );

U7a[x,ep] = (
  + 1 + Pi^2/6 - x - (Pi^2*x)/6 + HPL[{0}, x] + (Pi^2*HPL[{0}, x])/6 + HPL[{0, 0}, x] - x*HPL[{0, 0}, x] + HPL[{1, 0}, x] - x*HPL[{1, 0}, x] + HPL[{2, 0}, x] + HPL[{0, 0, 0}, x] + 2*Zeta[3]
  + ep*(7 - Pi^2/3 - (23*Pi^4)/180 - 7*x - (Pi^2*x)/3 + 4*HPL[{0}, x] - (5*Pi^2*HPL[{0}, x])/6 - x*HPL[{0}, x] + (5*Pi^2*x*HPL[{0}, x])/6 + 4*HPL[{1}, x] - (Pi^2*HPL[{1}, x])/3 - 4*x*HPL[{1}, x] + (Pi^2*x*HPL[{1}, x])/3 + 4*HPL[{2}, x] - (Pi^2*HPL[{2}, x])/3 + 4*HPL[{3}, x] - 4*x*HPL[{3}, x] + 4*HPL[{4}, x] + HPL[{0, 0}, x] - (7*Pi^2*HPL[{0, 0}, x])/6 - 2*x*HPL[{0, 0}, x] + 2*HPL[{1, 0}, x] - 2*x*HPL[{1, 0}, x] + 4*HPL[{1, 2}, x] - 4*x*HPL[{1, 2}, x] - HPL[{2, 0}, x] + x*HPL[{2, 0}, x] + 4*HPL[{2, 2}, x] - 3*HPL[{3, 0}, x] - 2*HPL[{0, 0, 0}, x] + 2*x*HPL[{0, 0, 0}, x] + HPL[{1, 0, 0}, x] - x*HPL[{1, 0, 0}, x] + 2*HPL[{1, 1, 0}, x] - 2*x*HPL[{1, 1, 0}, x] + HPL[{2, 0, 0}, x] + 2*HPL[{2, 1, 0}, x] - 4*HPL[{0, 0, 0, 0}, x] - Zeta[3] + x*Zeta[3] - 5*HPL[{0}, x]*Zeta[3])
  + ep^2*(33 - 3*Pi^2 + (2*Pi^4)/15 - 33*x + Pi^2*x - (2*Pi^4*x)/15 + 12*HPL[{0}, x] - 2*Pi^2*HPL[{0}, x] + (7*Pi^4*HPL[{0}, x])/18 - 7*x*HPL[{0}, x] + (5*Pi^2*x*HPL[{0}, x])/3 + 28*HPL[{1}, x] - (2*Pi^2*HPL[{1}, x])/3 - 28*x*HPL[{1}, x] + (2*Pi^2*x*HPL[{1}, x])/3 + 16*HPL[{2}, x] + (Pi^2*HPL[{2}, x])/3 - 4*x*HPL[{2}, x] - (Pi^2*x*HPL[{2}, x])/3 + 4*HPL[{3}, x] + Pi^2*HPL[{3}, x] - 8*x*HPL[{3}, x] - 8*HPL[{4}, x] + 8*x*HPL[{4}, x] - 16*HPL[{5}, x] + (Pi^2*HPL[{0, 0}, x])/2 - 3*x*HPL[{0, 0}, x] - (Pi^2*x*HPL[{0, 0}, x])/2 + 6*HPL[{1, 0}, x] - 2*Pi^2*HPL[{1, 0}, x] - 6*x*HPL[{1, 0}, x] + 2*Pi^2*x*HPL[{1, 0}, x] + 16*HPL[{1, 1}, x] - (2*Pi^2*HPL[{1, 1}, x])/3 - 16*x*HPL[{1, 1}, x] + (2*Pi^2*x*HPL[{1, 1}, x])/3 + 8*HPL[{1, 2}, x] - 8*x*HPL[{1, 2}, x] + 4*HPL[{1, 3}, x] - 4*x*HPL[{1, 3}, x] - 2*Pi^2*HPL[{2, 0}, x] + 2*x*HPL[{2, 0}, x] + 16*HPL[{2, 1}, x] - (2*Pi^2*HPL[{2, 1}, x])/3 - 4*HPL[{2, 2}, x] + 4*x*HPL[{2, 2}, x] + 4*HPL[{2, 3}, x] + 3*HPL[{3, 0}, x] - 3*x*HPL[{3, 0}, x] + 16*HPL[{3, 1}, x] - 16*x*HPL[{3, 1}, x] - 12*HPL[{3, 2}, x] + 9*HPL[{4, 0}, x] + 16*HPL[{4, 1}, x] - 3*HPL[{0, 0, 0}, x] + (17*Pi^2*HPL[{0, 0, 0}, x])/6 + 4*x*HPL[{0, 0, 0}, x] + 2*HPL[{1, 0, 0}, x] - 2*x*HPL[{1, 0, 0}, x] + 4*HPL[{1, 1, 0}, x] - 4*x*HPL[{1, 1, 0}, x] + 8*HPL[{1, 1, 2}, x] - 8*x*HPL[{1, 1, 2}, x] + 16*HPL[{1, 2, 1}, x] - 16*x*HPL[{1, 2, 1}, x] - HPL[{2, 0, 0}, x] + x*HPL[{2, 0, 0}, x] - 2*HPL[{2, 1, 0}, x] + 2*x*HPL[{2, 1, 0}, x] + 8*HPL[{2, 1, 2}, x] + 16*HPL[{2, 2, 1}, x] - 3*HPL[{3, 0, 0}, x] - 6*HPL[{3, 1, 0}, x] + 3*HPL[{0, 0, 0, 0}, x] - 3*x*HPL[{0, 0, 0, 0}, x] - 3*HPL[{1, 0, 0, 0}, x] + 3*x*HPL[{1, 0, 0, 0}, x] + 2*HPL[{1, 1, 0, 0}, x] - 2*x*HPL[{1, 1, 0, 0}, x] + 4*HPL[{1, 1, 1, 0}, x] - 4*x*HPL[{1, 1, 1, 0}, x] - 3*HPL[{2, 0, 0, 0}, x] + 2*HPL[{2, 1, 0, 0}, x] + 4*HPL[{2, 1, 1, 0}, x] + 11*HPL[{0, 0, 0, 0, 0}, x] - 10*Zeta[3] + Pi^2*Zeta[3] + 2*x*Zeta[3] - 7*HPL[{0}, x]*Zeta[3] + 7*x*HPL[{0}, x]*Zeta[3] - 10*HPL[{1}, x]*Zeta[3] + 10*x*HPL[{1}, x]*Zeta[3] - 10*HPL[{2}, x]*Zeta[3] + 3*HPL[{0, 0}, x]*Zeta[3] + 4*Zeta[5])
  + ep^3*(131 - (41*Pi^2)/3 + (22*Pi^4)/45 - (19*Pi^6)/270 - 131*x + (31*Pi^2*x)/3 - (4*Pi^4*x)/15 + 32*HPL[{0}, x] - (14*Pi^2*HPL[{0}, x])/3 + (4*Pi^4*HPL[{0}, x])/45 - 33*x*HPL[{0}, x] + (11*Pi^2*x*HPL[{0}, x])/3 - (4*Pi^4*x*HPL[{0}, x])/45 + 132*HPL[{1}, x] - (20*Pi^2*HPL[{1}, x])/3 + (22*Pi^4*HPL[{1}, x])/45 - 132*x*HPL[{1}, x] + (20*Pi^2*x*HPL[{1}, x])/3 - (22*Pi^4*x*HPL[{1}, x])/45 + 48*HPL[{2}, x] - (14*Pi^2*HPL[{2}, x])/3 + (22*Pi^4*HPL[{2}, x])/45 - 28*x*HPL[{2}, x] - (2*Pi^2*x*HPL[{2}, x])/3 - (17*Pi^2*HPL[{3}, x])/3 - 12*x*HPL[{3}, x] + (17*Pi^2*x*HPL[{3}, x])/3 - 12*HPL[{4}, x] - (23*Pi^2*HPL[{4}, x])/3 + 16*x*HPL[{4}, x] + 12*HPL[{5}, x] - 12*x*HPL[{5}, x] + 44*HPL[{6}, x] - 4*HPL[{0, 0}, x] + (4*Pi^2*HPL[{0, 0}, x])/3 - (31*Pi^4*HPL[{0, 0}, x])/45 - x*HPL[{0, 0}, x] - Pi^2*x*HPL[{0, 0}, x] + 22*HPL[{1, 0}, x] - 4*Pi^2*HPL[{1, 0}, x] - 22*x*HPL[{1, 0}, x] + 4*Pi^2*x*HPL[{1, 0}, x] + 112*HPL[{1, 1}, x] - (4*Pi^2*HPL[{1, 1}, x])/3 - 112*x*HPL[{1, 1}, x] + (4*Pi^2*x*HPL[{1, 1}, x])/3 + 24*HPL[{1, 2}, x] - (14*Pi^2*HPL[{1, 2}, x])/3 - 24*x*HPL[{1, 2}, x] + (14*Pi^2*x*HPL[{1, 2}, x])/3 + 8*HPL[{1, 3}, x] - 8*x*HPL[{1, 3}, x] - 12*HPL[{1, 4}, x] + 12*x*HPL[{1, 4}, x] + 4*HPL[{2, 0}, x] + 2*Pi^2*HPL[{2, 0}, x] + 2*x*HPL[{2, 0}, x] - 2*Pi^2*x*HPL[{2, 0}, x] + 64*HPL[{2, 1}, x] + (2*Pi^2*HPL[{2, 1}, x])/3 - 16*x*HPL[{2, 1}, x] - (2*Pi^2*x*HPL[{2, 1}, x])/3 - (14*Pi^2*HPL[{2, 2}, x])/3 + 8*x*HPL[{2, 2}, x] - 4*HPL[{2, 3}, x] + 4*x*HPL[{2, 3}, x] - 12*HPL[{2, 4}, x] + 4*HPL[{3, 0}, x] + 6*Pi^2*HPL[{3, 0}, x] - 6*x*HPL[{3, 0}, x] + 16*HPL[{3, 1}, x] + 2*Pi^2*HPL[{3, 1}, x] - 32*x*HPL[{3, 1}, x] + 12*HPL[{3, 2}, x] - 12*x*HPL[{3, 2}, x] - 12*HPL[{3, 3}, x] - 5*HPL[{4, 0}, x] + 5*x*HPL[{4, 0}, x] - 32*HPL[{4, 1}, x] + 32*x*HPL[{4, 1}, x] + 36*HPL[{4, 2}, x] - 23*HPL[{5, 0}, x] - 64*HPL[{5, 1}, x] - 4*HPL[{0, 0, 0}, x] - (Pi^2*HPL[{0, 0, 0}, x])/6 + 7*x*HPL[{0, 0, 0}, x] + (Pi^2*x*HPL[{0, 0, 0}, x])/6 + 2*HPL[{1, 0, 0}, x] + (4*Pi^2*HPL[{1, 0, 0}, x])/3 - 2*x*HPL[{1, 0, 0}, x] - (4*Pi^2*x*HPL[{1, 0, 0}, x])/3 + 16*HPL[{1, 1, 0}, x] - 4*Pi^2*HPL[{1, 1, 0}, x] - 16*x*HPL[{1, 1, 0}, x] + 4*Pi^2*x*HPL[{1, 1, 0}, x] + 64*HPL[{1, 1, 1}, x] - (4*Pi^2*HPL[{1, 1, 1}, x])/3 - 64*x*HPL[{1, 1, 1}, x] + (4*Pi^2*x*HPL[{1, 1, 1}, x])/3 + 16*HPL[{1, 1, 2}, x] - 16*x*HPL[{1, 1, 2}, x] + 8*HPL[{1, 1, 3}, x] - 8*x*HPL[{1, 1, 3}, x] + 32*HPL[{1, 2, 1}, x] - 32*x*HPL[{1, 2, 1}, x] + 4*HPL[{1, 3, 0}, x] - 4*x*HPL[{1, 3, 0}, x] + 16*HPL[{1, 3, 1}, x] - 16*x*HPL[{1, 3, 1}, x] - 4*HPL[{2, 0, 0}, x] + (4*Pi^2*HPL[{2, 0, 0}, x])/3 + 2*x*HPL[{2, 0, 0}, x] + 4*HPL[{2, 1, 0}, x] - 4*Pi^2*HPL[{2, 1, 0}, x] + 4*x*HPL[{2, 1, 0}, x] + 64*HPL[{2, 1, 1}, x] - (4*Pi^2*HPL[{2, 1, 1}, x])/3 - 8*HPL[{2, 1, 2}, x] + 8*x*HPL[{2, 1, 2}, x] + 8*HPL[{2, 1, 3}, x] - 16*HPL[{2, 2, 1}, x] + 16*x*HPL[{2, 2, 1}, x] + 4*HPL[{2, 3, 0}, x] + 16*HPL[{2, 3, 1}, x] - HPL[{3, 0, 0}, x] + x*HPL[{3, 0, 0}, x] + 10*HPL[{3, 1, 0}, x] - 10*x*HPL[{3, 1, 0}, x] + 64*HPL[{3, 1, 1}, x] - 64*x*HPL[{3, 1, 1}, x] - 24*HPL[{3, 1, 2}, x] - 48*HPL[{3, 2, 1}, x] + 5*HPL[{4, 0, 0}, x] + 22*HPL[{4, 1, 0}, x] + 64*HPL[{4, 1, 1}, x] + 5*HPL[{0, 0, 0, 0}, x] - (35*Pi^2*HPL[{0, 0, 0, 0}, x])/6 - 6*x*HPL[{0, 0, 0, 0}, x] - 6*HPL[{1, 0, 0, 0}, x] + 6*x*HPL[{1, 0, 0, 0}, x] + 4*HPL[{1, 1, 0, 0}, x] - 4*x*HPL[{1, 1, 0, 0}, x] + 8*HPL[{1, 1, 1, 0}, x] - 8*x*HPL[{1, 1, 1, 0}, x] + 16*HPL[{1, 1, 1, 2}, x] - 16*x*HPL[{1, 1, 1, 2}, x] + 32*HPL[{1, 1, 2, 1}, x] - 32*x*HPL[{1, 1, 2, 1}, x] - 4*HPL[{1, 2, 0, 0}, x] + 4*x*HPL[{1, 2, 0, 0}, x] + 4*HPL[{1, 2, 1, 0}, x] - 4*x*HPL[{1, 2, 1, 0}, x] + 64*HPL[{1, 2, 1, 1}, x] - 64*x*HPL[{1, 2, 1, 1}, x] + 3*HPL[{2, 0, 0, 0}, x] - 3*x*HPL[{2, 0, 0, 0}, x] - 2*HPL[{2, 1, 0, 0}, x] + 2*x*HPL[{2, 1, 0, 0}, x] - 4*HPL[{2, 1, 1, 0}, x] + 4*x*HPL[{2, 1, 1, 0}, x] + 16*HPL[{2, 1, 1, 2}, x] + 32*HPL[{2, 1, 2, 1}, x] - 4*HPL[{2, 2, 0, 0}, x] + 4*HPL[{2, 2, 1, 0}, x] + 64*HPL[{2, 2, 1, 1}, x] + 9*HPL[{3, 0, 0, 0}, x] - 6*HPL[{3, 1, 0, 0}, x] - 12*HPL[{3, 1, 1, 0}, x] - 4*HPL[{0, 0, 0, 0, 0}, x] + 4*x*HPL[{0, 0, 0, 0, 0}, x] + 5*HPL[{1, 0, 0, 0, 0}, x] - 5*x*HPL[{1, 0, 0, 0, 0}, x] - 6*HPL[{1, 1, 0, 0, 0}, x] + 6*x*HPL[{1, 1, 0, 0, 0}, x] + 4*HPL[{1, 1, 1, 0, 0}, x] - 4*x*HPL[{1, 1, 1, 0, 0}, x] + 8*HPL[{1, 1, 1, 1, 0}, x] - 8*x*HPL[{1, 1, 1, 1, 0}, x] + 5*HPL[{2, 0, 0, 0, 0}, x] - 6*HPL[{2, 1, 0, 0, 0}, x] + 4*HPL[{2, 1, 1, 0, 0}, x] + 8*HPL[{2, 1, 1, 1, 0}, x] - 26*HPL[{0, 0, 0, 0, 0, 0}, x] - 46*Zeta[3] + (16*Pi^2*Zeta[3])/3 + 22*x*Zeta[3] - (16*Pi^2*x*Zeta[3])/3 - 16*HPL[{0}, x]*Zeta[3] + (10*Pi^2*HPL[{0}, x]*Zeta[3])/3 + 14*x*HPL[{0}, x]*Zeta[3] - 20*HPL[{1}, x]*Zeta[3] + 20*x*HPL[{1}, x]*Zeta[3] + 10*HPL[{2}, x]*Zeta[3] - 10*x*HPL[{2}, x]*Zeta[3] + 30*HPL[{3}, x]*Zeta[3] + 5*HPL[{0, 0}, x]*Zeta[3] - 5*x*HPL[{0, 0}, x]*Zeta[3] - 16*HPL[{1, 0}, x]*Zeta[3] + 16*x*HPL[{1, 0}, x]*Zeta[3] - 20*HPL[{1, 1}, x]*Zeta[3] + 20*x*HPL[{1, 1}, x]*Zeta[3] - 16*HPL[{2, 0}, x]*Zeta[3] - 20*HPL[{2, 1}, x]*Zeta[3] - HPL[{0, 0, 0}, x]*Zeta[3] + 11*Zeta[3]^2 - 17*Zeta[5] + 17*x*Zeta[5] - 25*HPL[{0}, x]*Zeta[5])
  + ALARM*ep^4 );
U7a[x -> 1,ep] = (1-x)^(3 - 4*ep) * (
  - 1/6 - (25*ep)/18 + ep^2*(-212/27 + (2*Pi^2)/9) + ep^3*(-3062/81 + (50*Pi^2)/27 + (10*Zeta[3])/3)
  + ALARM*ep^4 );

U8[x,ep] = (
  - HPL[{0}, x]/ep^2
  + ((2*Pi^2)/3 - 4*HPL[{2}, x] + 2*HPL[{0, 0}, x])/ep
  - (Pi^2*HPL[{0}, x])/3 + 8*HPL[{3}, x] - 2*HPL[{2, 0}, x] - 16*HPL[{2, 1}, x] - 4*HPL[{0, 0, 0}, x] + 4*Zeta[3]
  + ep*(-Pi^4/15 + (16*Pi^2*HPL[{2}, x])/3 - 16*HPL[{4}, x] + (2*Pi^2*HPL[{0, 0}, x])/3 - 8*HPL[{2, 2}, x] + 4*HPL[{3, 0}, x] + 32*HPL[{3, 1}, x] + 2*HPL[{2, 0, 0}, x] - 8*HPL[{2, 1, 0}, x] - 64*HPL[{2, 1, 1}, x] + 8*HPL[{0, 0, 0, 0}, x] + 2*HPL[{0}, x]*Zeta[3])
  + ep^2*(-(Pi^4*HPL[{0}, x])/10 - (32*Pi^2*HPL[{3}, x])/3 + 32*HPL[{5}, x] + (2*Pi^2*HPL[{2, 0}, x])/3 + (64*Pi^2*HPL[{2, 1}, x])/3 + 8*HPL[{2, 3}, x] + 16*HPL[{3, 2}, x] - 8*HPL[{4, 0}, x] - 64*HPL[{4, 1}, x] - (4*Pi^2*HPL[{0, 0, 0}, x])/3 - 32*HPL[{2, 1, 2}, x] - 4*HPL[{2, 2, 0}, x] - 32*HPL[{2, 2, 1}, x] - 4*HPL[{3, 0, 0}, x] + 16*HPL[{3, 1, 0}, x] + 128*HPL[{3, 1, 1}, x] - 2*HPL[{2, 0, 0, 0}, x] + 8*HPL[{2, 1, 0, 0}, x] - 32*HPL[{2, 1, 1, 0}, x] - 256*HPL[{2, 1, 1, 1}, x] - 16*HPL[{0, 0, 0, 0, 0}, x] - (16*Pi^2*Zeta[3])/3 + 56*HPL[{2}, x]*Zeta[3] - 4*HPL[{0, 0}, x]*Zeta[3] + 24*Zeta[5])
  + ALARM*ep^3 );
U8[x -> 1,ep] = (1-x)^(1 - 4*ep) * (
  + ep^(-2) + 4/ep + 16 - (4*Pi^2)/3 + ep*(64 - (16*Pi^2)/3 - 20*Zeta[3]) + ep^2*(256 - (64*Pi^2)/3 + (8*Pi^4)/45 - 80*Zeta[3])
  + ALARM*ep^3 );

U9[x,ep] = (
  - Pi^2/3 - 4*HPL[{-1, 0}, x] + 2*HPL[{0, 0}, x]
  + ep*((-2*Pi^2)/3 + 3*Pi^2*HPL[{-1}, x] - Pi^2*HPL[{0}, x] + 8*HPL[{3}, x] + 4*HPL[{-2, 0}, x] - 8*HPL[{-1, 0}, x] - 16*HPL[{-1, 2}, x] + 4*HPL[{0, 0}, x] + 4*HPL[{-1, -1, 0}, x] + 2*HPL[{-1, 0, 0}, x] - 4*HPL[{0, 0, 0}, x] - 14*Zeta[3])
  + ep^2*((-4*Pi^2)/3 + (8*Pi^4)/45 - 3*Pi^2*HPL[{-2}, x] + 6*Pi^2*HPL[{-1}, x] - 2*Pi^2*HPL[{0}, x] + 16*HPL[{3}, x] - 16*HPL[{4}, x] - 4*HPL[{-3, 0}, x] + 8*HPL[{-2, 0}, x] + 16*HPL[{-2, 2}, x] - 3*Pi^2*HPL[{-1, -1}, x] - 16*HPL[{-1, 0}, x] + (7*Pi^2*HPL[{-1, 0}, x])/3 - 32*HPL[{-1, 2}, x] + 8*HPL[{-1, 3}, x] + 8*HPL[{0, 0}, x] + (Pi^2*HPL[{0, 0}, x])/3 + 4*HPL[{3, 0}, x] + 32*HPL[{3, 1}, x] - 4*HPL[{-2, -1, 0}, x] - 2*HPL[{-2, 0, 0}, x] - 4*HPL[{-1, -2, 0}, x] + 8*HPL[{-1, -1, 0}, x] + 16*HPL[{-1, -1, 2}, x] + 4*HPL[{-1, 0, 0}, x] - 8*HPL[{-1, 2, 0}, x] - 64*HPL[{-1, 2, 1}, x] - 8*HPL[{0, 0, 0}, x] - 4*HPL[{-1, -1, -1, 0}, x] - 2*HPL[{-1, -1, 0, 0}, x] + 6*HPL[{0, 0, 0, 0}, x] - 28*Zeta[3] + 46*HPL[{-1}, x]*Zeta[3] - 2*HPL[{0}, x]*Zeta[3])
  + ep^3*((-8*Pi^2)/3 + (16*Pi^4)/45 + 3*Pi^2*HPL[{-3}, x] - 6*Pi^2*HPL[{-2}, x] + 12*Pi^2*HPL[{-1}, x] - (16*Pi^4*HPL[{-1}, x])/15 - 4*Pi^2*HPL[{0}, x] + (4*Pi^4*HPL[{0}, x])/15 + 32*HPL[{3}, x] - (32*Pi^2*HPL[{3}, x])/3 - 32*HPL[{4}, x] + 24*HPL[{5}, x] + 4*HPL[{-4, 0}, x] - 8*HPL[{-3, 0}, x] - 16*HPL[{-3, 2}, x] + 3*Pi^2*HPL[{-2, -1}, x] + 16*HPL[{-2, 0}, x] - (7*Pi^2*HPL[{-2, 0}, x])/3 + 32*HPL[{-2, 2}, x] - 8*HPL[{-2, 3}, x] + 3*Pi^2*HPL[{-1, -2}, x] - 6*Pi^2*HPL[{-1, -1}, x] - 32*HPL[{-1, 0}, x] + (14*Pi^2*HPL[{-1, 0}, x])/3 - 64*HPL[{-1, 2}, x] + (64*Pi^2*HPL[{-1, 2}, x])/3 + 16*HPL[{-1, 3}, x] + 16*HPL[{0, 0}, x] + (2*Pi^2*HPL[{0, 0}, x])/3 + 8*HPL[{3, 0}, x] + 64*HPL[{3, 1}, x] + 16*HPL[{3, 2}, x] - 8*HPL[{4, 0}, x] - 64*HPL[{4, 1}, x] + 4*HPL[{-3, -1, 0}, x] + 2*HPL[{-3, 0, 0}, x] + 4*HPL[{-2, -2, 0}, x] - 8*HPL[{-2, -1, 0}, x] - 16*HPL[{-2, -1, 2}, x] - 4*HPL[{-2, 0, 0}, x] + 8*HPL[{-2, 2, 0}, x] + 64*HPL[{-2, 2, 1}, x] + 4*HPL[{-1, -3, 0}, x] - 8*HPL[{-1, -2, 0}, x] - 16*HPL[{-1, -2, 2}, x] + 3*Pi^2*HPL[{-1, -1, -1}, x] + 16*HPL[{-1, -1, 0}, x] - (7*Pi^2*HPL[{-1, -1, 0}, x])/3 + 32*HPL[{-1, -1, 2}, x] - 8*HPL[{-1, -1, 3}, x] + 8*HPL[{-1, 0, 0}, x] - (5*Pi^2*HPL[{-1, 0, 0}, x])/3 - 16*HPL[{-1, 2, 0}, x] - 128*HPL[{-1, 2, 1}, x] - 32*HPL[{-1, 2, 2}, x] + 4*HPL[{-1, 3, 0}, x] + 32*HPL[{-1, 3, 1}, x] - 16*HPL[{0, 0, 0}, x] + (Pi^2*HPL[{0, 0, 0}, x])/3 - 4*HPL[{3, 0, 0}, x] + 16*HPL[{3, 1, 0}, x] + 128*HPL[{3, 1, 1}, x] + 4*HPL[{-2, -1, -1, 0}, x] + 2*HPL[{-2, -1, 0, 0}, x] + 4*HPL[{-1, -2, -1, 0}, x] + 2*HPL[{-1, -2, 0, 0}, x] + 4*HPL[{-1, -1, -2, 0}, x] - 8*HPL[{-1, -1, -1, 0}, x] - 16*HPL[{-1, -1, -1, 2}, x] - 4*HPL[{-1, -1, 0, 0}, x] + 8*HPL[{-1, -1, 2, 0}, x] + 64*HPL[{-1, -1, 2, 1}, x] + 8*HPL[{-1, 2, 0, 0}, x] - 32*HPL[{-1, 2, 1, 0}, x] - 256*HPL[{-1, 2, 1, 1}, x] + 12*HPL[{0, 0, 0, 0}, x] + 4*HPL[{-1, -1, -1, -1, 0}, x] + 2*HPL[{-1, -1, -1, 0, 0}, x] - 2*HPL[{-1, 0, 0, 0, 0}, x] - 8*HPL[{0, 0, 0, 0, 0}, x] - 56*Zeta[3] + 16*Pi^2*Zeta[3] - 46*HPL[{-2}, x]*Zeta[3] + 92*HPL[{-1}, x]*Zeta[3] - 4*HPL[{0}, x]*Zeta[3] - 46*HPL[{-1, -1}, x]*Zeta[3] + 10*HPL[{-1, 0}, x]*Zeta[3] - 2*HPL[{0, 0}, x]*Zeta[3] - 94*Zeta[5])
  + ALARM*ep^4 );
U9[x -> 1,ep] = (1-x)^(3 - 4*ep) * (
  + 1/3 + (22*ep)/9 + ep^2*(340/27 - (4*Pi^2)/9) + ep^3*(4600/81 - (88*Pi^2)/27 - (20*Zeta[3])/3)
  + ALARM*ep^4 );

U10[x,ep] = (
  + Pi^2/6 + HPL[{0, 0}, x] + HPL[{1, 0}, x]
  + ep*(Pi^2/3 - (5*Pi^2*HPL[{0}, x])/6 - (Pi^2*HPL[{1}, x])/3 + 4*HPL[{3}, x] + 2*HPL[{0, 0}, x] + 2*HPL[{1, 0}, x] + 4*HPL[{1, 2}, x] - HPL[{2, 0}, x] - 2*HPL[{0, 0, 0}, x] + HPL[{1, 0, 0}, x] + 2*HPL[{1, 1, 0}, x] - Zeta[3])
  + ep^2*((2*Pi^2)/3 + (2*Pi^4)/15 - (5*Pi^2*HPL[{0}, x])/3 - (2*Pi^2*HPL[{1}, x])/3 + (Pi^2*HPL[{2}, x])/3 + 8*HPL[{3}, x] - 8*HPL[{4}, x] + 4*HPL[{0, 0}, x] + (Pi^2*HPL[{0, 0}, x])/2 + 4*HPL[{1, 0}, x] - 2*Pi^2*HPL[{1, 0}, x] - (2*Pi^2*HPL[{1, 1}, x])/3 + 8*HPL[{1, 2}, x] + 4*HPL[{1, 3}, x] - 2*HPL[{2, 0}, x] - 4*HPL[{2, 2}, x] + 3*HPL[{3, 0}, x] + 16*HPL[{3, 1}, x] - 4*HPL[{0, 0, 0}, x] + 2*HPL[{1, 0, 0}, x] + 4*HPL[{1, 1, 0}, x] + 8*HPL[{1, 1, 2}, x] + 16*HPL[{1, 2, 1}, x] - HPL[{2, 0, 0}, x] - 2*HPL[{2, 1, 0}, x] + 3*HPL[{0, 0, 0, 0}, x] - 3*HPL[{1, 0, 0, 0}, x] + 2*HPL[{1, 1, 0, 0}, x] + 4*HPL[{1, 1, 1, 0}, x] - 2*Zeta[3] - 7*HPL[{0}, x]*Zeta[3] - 10*HPL[{1}, x]*Zeta[3])
  + ep^3*((4*Pi^2)/3 + (4*Pi^4)/15 - (10*Pi^2*HPL[{0}, x])/3 + (4*Pi^4*HPL[{0}, x])/45 - (4*Pi^2*HPL[{1}, x])/3 + (22*Pi^4*HPL[{1}, x])/45 + (2*Pi^2*HPL[{2}, x])/3 + 16*HPL[{3}, x] - (17*Pi^2*HPL[{3}, x])/3 - 16*HPL[{4}, x] + 12*HPL[{5}, x] + 8*HPL[{0, 0}, x] + Pi^2*HPL[{0, 0}, x] + 8*HPL[{1, 0}, x] - 4*Pi^2*HPL[{1, 0}, x] - (4*Pi^2*HPL[{1, 1}, x])/3 + 16*HPL[{1, 2}, x] - (14*Pi^2*HPL[{1, 2}, x])/3 + 8*HPL[{1, 3}, x] - 12*HPL[{1, 4}, x] - 4*HPL[{2, 0}, x] + 2*Pi^2*HPL[{2, 0}, x] + (2*Pi^2*HPL[{2, 1}, x])/3 - 8*HPL[{2, 2}, x] - 4*HPL[{2, 3}, x] + 6*HPL[{3, 0}, x] + 32*HPL[{3, 1}, x] + 12*HPL[{3, 2}, x] - 5*HPL[{4, 0}, x] - 32*HPL[{4, 1}, x] - 8*HPL[{0, 0, 0}, x] - (Pi^2*HPL[{0, 0, 0}, x])/6 + 4*HPL[{1, 0, 0}, x] + (4*Pi^2*HPL[{1, 0, 0}, x])/3 + 8*HPL[{1, 1, 0}, x] - 4*Pi^2*HPL[{1, 1, 0}, x] - (4*Pi^2*HPL[{1, 1, 1}, x])/3 + 16*HPL[{1, 1, 2}, x] + 8*HPL[{1, 1, 3}, x] + 32*HPL[{1, 2, 1}, x] + 4*HPL[{1, 3, 0}, x] + 16*HPL[{1, 3, 1}, x] - 2*HPL[{2, 0, 0}, x] - 4*HPL[{2, 1, 0}, x] - 8*HPL[{2, 1, 2}, x] - 16*HPL[{2, 2, 1}, x] - HPL[{3, 0, 0}, x] + 10*HPL[{3, 1, 0}, x] + 64*HPL[{3, 1, 1}, x] + 6*HPL[{0, 0, 0, 0}, x] - 6*HPL[{1, 0, 0, 0}, x] + 4*HPL[{1, 1, 0, 0}, x] + 8*HPL[{1, 1, 1, 0}, x] + 16*HPL[{1, 1, 1, 2}, x] + 32*HPL[{1, 1, 2, 1}, x] - 4*HPL[{1, 2, 0, 0}, x] + 4*HPL[{1, 2, 1, 0}, x] + 64*HPL[{1, 2, 1, 1}, x] + 3*HPL[{2, 0, 0, 0}, x] - 2*HPL[{2, 1, 0, 0}, x] - 4*HPL[{2, 1, 1, 0}, x] - 4*HPL[{0, 0, 0, 0, 0}, x] + 5*HPL[{1, 0, 0, 0, 0}, x] - 6*HPL[{1, 1, 0, 0, 0}, x] + 4*HPL[{1, 1, 1, 0, 0}, x] + 8*HPL[{1, 1, 1, 1, 0}, x] - 4*Zeta[3] + (16*Pi^2*Zeta[3])/3 - 14*HPL[{0}, x]*Zeta[3] - 20*HPL[{1}, x]*Zeta[3] + 10*HPL[{2}, x]*Zeta[3] + 5*HPL[{0, 0}, x]*Zeta[3] - 16*HPL[{1, 0}, x]*Zeta[3] - 20*HPL[{1, 1}, x]*Zeta[3] - 17*Zeta[5])
  + ep^4*((8*Pi^2)/3 + (8*Pi^4)/15 - (83*Pi^6)/3780 - (20*Pi^2*HPL[{0}, x])/3 + (8*Pi^4*HPL[{0}, x])/45 - (8*Pi^2*HPL[{1}, x])/3 + (44*Pi^4*HPL[{1}, x])/45 + (4*Pi^2*HPL[{2}, x])/3 - (22*Pi^4*HPL[{2}, x])/45 + 32*HPL[{3}, x] - (34*Pi^2*HPL[{3}, x])/3 - 32*HPL[{4}, x] + 11*Pi^2*HPL[{4}, x] + 24*HPL[{5}, x] - 16*HPL[{6}, x] + 16*HPL[{0, 0}, x] + 2*Pi^2*HPL[{0, 0}, x] - (7*Pi^4*HPL[{0, 0}, x])/90 + 16*HPL[{1, 0}, x] - 8*Pi^2*HPL[{1, 0}, x] + (17*Pi^4*HPL[{1, 0}, x])/90 - (8*Pi^2*HPL[{1, 1}, x])/3 + (44*Pi^4*HPL[{1, 1}, x])/45 + 32*HPL[{1, 2}, x] - (28*Pi^2*HPL[{1, 2}, x])/3 + 16*HPL[{1, 3}, x] - 6*Pi^2*HPL[{1, 3}, x] - 24*HPL[{1, 4}, x] + 20*HPL[{1, 5}, x] - 8*HPL[{2, 0}, x] + 4*Pi^2*HPL[{2, 0}, x] + (4*Pi^2*HPL[{2, 1}, x])/3 - 16*HPL[{2, 2}, x] + (14*Pi^2*HPL[{2, 2}, x])/3 - 8*HPL[{2, 3}, x] + 12*HPL[{2, 4}, x] + 12*HPL[{3, 0}, x] - (8*Pi^2*HPL[{3, 0}, x])/3 + 64*HPL[{3, 1}, x] - 22*Pi^2*HPL[{3, 1}, x] + 24*HPL[{3, 2}, x] - 4*HPL[{3, 3}, x] - 10*HPL[{4, 0}, x] - 64*HPL[{4, 1}, x] - 20*HPL[{4, 2}, x] + 7*HPL[{5, 0}, x] + 48*HPL[{5, 1}, x] - 16*HPL[{0, 0, 0}, x] - (Pi^2*HPL[{0, 0, 0}, x])/3 + 8*HPL[{1, 0, 0}, x] + (8*Pi^2*HPL[{1, 0, 0}, x])/3 + 16*HPL[{1, 1, 0}, x] - 8*Pi^2*HPL[{1, 1, 0}, x] - (8*Pi^2*HPL[{1, 1, 1}, x])/3 + 32*HPL[{1, 1, 2}, x] - (28*Pi^2*HPL[{1, 1, 2}, x])/3 + 16*HPL[{1, 1, 3}, x] - 24*HPL[{1, 1, 4}, x] + (10*Pi^2*HPL[{1, 2, 0}, x])/3 + 64*HPL[{1, 2, 1}, x] - 20*Pi^2*HPL[{1, 2, 1}, x] - 16*HPL[{1, 2, 3}, x] + 8*HPL[{1, 3, 0}, x] + 32*HPL[{1, 3, 1}, x] + 16*HPL[{1, 3, 2}, x] - 8*HPL[{1, 4, 0}, x] - 48*HPL[{1, 4, 1}, x] - 4*HPL[{2, 0, 0}, x] - (4*Pi^2*HPL[{2, 0, 0}, x])/3 - 8*HPL[{2, 1, 0}, x] + 4*Pi^2*HPL[{2, 1, 0}, x] + (4*Pi^2*HPL[{2, 1, 1}, x])/3 - 16*HPL[{2, 1, 2}, x] - 8*HPL[{2, 1, 3}, x] - 32*HPL[{2, 2, 1}, x] - 4*HPL[{2, 3, 0}, x] - 16*HPL[{2, 3, 1}, x] - 2*HPL[{3, 0, 0}, x] + 20*HPL[{3, 1, 0}, x] + 128*HPL[{3, 1, 1}, x] + 40*HPL[{3, 1, 2}, x] + 4*HPL[{3, 2, 0}, x] + 48*HPL[{3, 2, 1}, x] + 3*HPL[{4, 0, 0}, x] - 18*HPL[{4, 1, 0}, x] - 128*HPL[{4, 1, 1}, x] + 12*HPL[{0, 0, 0, 0}, x] - (Pi^2*HPL[{0, 0, 0, 0}, x])/6 - 12*HPL[{1, 0, 0, 0}, x] - (2*Pi^2*HPL[{1, 0, 0, 0}, x])/3 + 8*HPL[{1, 1, 0, 0}, x] + (8*Pi^2*HPL[{1, 1, 0, 0}, x])/3 + 16*HPL[{1, 1, 1, 0}, x] - 8*Pi^2*HPL[{1, 1, 1, 0}, x] - (8*Pi^2*HPL[{1, 1, 1, 1}, x])/3 + 32*HPL[{1, 1, 1, 2}, x] + 16*HPL[{1, 1, 1, 3}, x] + 64*HPL[{1, 1, 2, 1}, x] + 8*HPL[{1, 1, 3, 0}, x] + 32*HPL[{1, 1, 3, 1}, x] - 8*HPL[{1, 2, 0, 0}, x] + 8*HPL[{1, 2, 1, 0}, x] + 128*HPL[{1, 2, 1, 1}, x] + 16*HPL[{1, 2, 1, 2}, x] + 4*HPL[{1, 2, 2, 0}, x] + 12*HPL[{1, 3, 1, 0}, x] + 64*HPL[{1, 3, 1, 1}, x] + 6*HPL[{2, 0, 0, 0}, x] - 4*HPL[{2, 1, 0, 0}, x] - 8*HPL[{2, 1, 1, 0}, x] - 16*HPL[{2, 1, 1, 2}, x] - 32*HPL[{2, 1, 2, 1}, x] + 4*HPL[{2, 2, 0, 0}, x] - 4*HPL[{2, 2, 1, 0}, x] - 64*HPL[{2, 2, 1, 1}, x] - HPL[{3, 0, 0, 0}, x] - 6*HPL[{3, 1, 0, 0}, x] + 36*HPL[{3, 1, 1, 0}, x] + 256*HPL[{3, 1, 1, 1}, x] - 8*HPL[{0, 0, 0, 0, 0}, x] + 10*HPL[{1, 0, 0, 0, 0}, x] - 12*HPL[{1, 1, 0, 0, 0}, x] + 8*HPL[{1, 1, 1, 0, 0}, x] + 16*HPL[{1, 1, 1, 1, 0}, x] + 32*HPL[{1, 1, 1, 1, 2}, x] + 64*HPL[{1, 1, 1, 2, 1}, x] - 8*HPL[{1, 1, 2, 0, 0}, x] + 8*HPL[{1, 1, 2, 1, 0}, x] + 128*HPL[{1, 1, 2, 1, 1}, x] + 8*HPL[{1, 2, 0, 0, 0}, x] - 12*HPL[{1, 2, 1, 0, 0}, x] + 24*HPL[{1, 2, 1, 1, 0}, x] + 256*HPL[{1, 2, 1, 1, 1}, x] - 5*HPL[{2, 0, 0, 0, 0}, x] + 6*HPL[{2, 1, 0, 0, 0}, x] - 4*HPL[{2, 1, 1, 0, 0}, x] - 8*HPL[{2, 1, 1, 1, 0}, x] + 5*HPL[{0, 0, 0, 0, 0, 0}, x] - 7*HPL[{1, 0, 0, 0, 0, 0}, x] + 10*HPL[{1, 1, 0, 0, 0, 0}, x] - 12*HPL[{1, 1, 1, 0, 0, 0}, x] + 8*HPL[{1, 1, 1, 1, 0, 0}, x] + 16*HPL[{1, 1, 1, 1, 1, 0}, x] - 8*Zeta[3] + (32*Pi^2*Zeta[3])/3 - 28*HPL[{0}, x]*Zeta[3] + 4*Pi^2*HPL[{0}, x]*Zeta[3] - 40*HPL[{1}, x]*Zeta[3] + 20*Pi^2*HPL[{1}, x]*Zeta[3] + 20*HPL[{2}, x]*Zeta[3] - 66*HPL[{3}, x]*Zeta[3] + 10*HPL[{0, 0}, x]*Zeta[3] - 32*HPL[{1, 0}, x]*Zeta[3] - 40*HPL[{1, 1}, x]*Zeta[3] - 36*HPL[{1, 2}, x]*Zeta[3] + 16*HPL[{2, 0}, x]*Zeta[3] + 20*HPL[{2, 1}, x]*Zeta[3] - 3*HPL[{0, 0, 0}, x]*Zeta[3] + 12*HPL[{1, 0, 0}, x]*Zeta[3] - 32*HPL[{1, 1, 0}, x]*Zeta[3] - 40*HPL[{1, 1, 1}, x]*Zeta[3] + 34*Zeta[3]^2 - 34*Zeta[5] - 31*HPL[{0}, x]*Zeta[5] - 82*HPL[{1}, x]*Zeta[5])
  + ALARM*ep^5 );
U10[x -> 1,ep] = (1-x)^(1 - 4*ep) * (
  + 1 + 8*ep + ep^2*(44 - (4*Pi^2)/3) + ep^3*(208 - (32*Pi^2)/3 - 20*Zeta[3]) + ep^4*(912 - (176*Pi^2)/3 + (8*Pi^4)/45 - 160*Zeta[3])
  + ALARM*ep^5 );

U11[x,ep] = (
  + HPL[{0}, x]/(ep^2*(1 + x))
  + (-(Pi^2/(1 + x)) + (4*HPL[{2}, x])/(1 + x) - (4*HPL[{-1, 0}, x])/(1 + x) + (2*HPL[{0, 0}, x])/(1 + x))/ep
  - (8*Pi^2*HPL[{0}, x])/(3*(1 + x)) + (8*HPL[{3}, x])/(1 + x) - (16*HPL[{-1, 2}, x])/(1 + x) + (4*Pi^2*HPL[{-1}, x])/(1 + x) + (16*HPL[{2, 1}, x])/(1 + x) + (16*HPL[{-1, -1, 0}, x])/(1 + x) - (8*HPL[{-1, 0, 0}, x])/(1 + x) - (6*HPL[{0, 0, 0}, x])/(1 + x) - (30*Zeta[3])/(1 + x)
  + ep*(Pi^4/(1 + x) - (14*Pi^2*HPL[{2}, x])/(3*(1 + x)) - (24*HPL[{4}, x])/(1 + x) - (16*Pi^2*HPL[{-1, -1}, x])/(1 + x) + (32*Pi^2*HPL[{-1, 0}, x])/(3*(1 + x)) - (32*HPL[{-1, 3}, x])/(1 + x) + (7*Pi^2*HPL[{0, 0}, x])/(3*(1 + x)) + (6*HPL[{3, 0}, x])/(1 + x) + (32*HPL[{3, 1}, x])/(1 + x) + (64*HPL[{-1, -1, 2}, x])/(1 + x) - (64*HPL[{-1, 2, 1}, x])/(1 + x) - (4*HPL[{2, 0, 0}, x])/(1 + x) + (4*HPL[{2, 1, 0}, x])/(1 + x) + (64*HPL[{2, 1, 1}, x])/(1 + x) - (64*HPL[{-1, -1, -1, 0}, x])/(1 + x) + (32*HPL[{-1, -1, 0, 0}, x])/(1 + x) + (24*HPL[{-1, 0, 0, 0}, x])/(1 + x) + (12*HPL[{0, 0, 0, 0}, x])/(1 + x) + (120*HPL[{-1}, x]*Zeta[3])/(1 + x) - (16*HPL[{0}, x]*Zeta[3])/(1 + x))
  + ep^2*((-4*Pi^4*HPL[{-1}, x])/(1 + x) + (Pi^4*HPL[{0}, x])/(10*(1 + x)) - (34*Pi^2*HPL[{3}, x])/(3*(1 + x)) + (48*HPL[{5}, x])/(1 + x) + (56*Pi^2*HPL[{-1, 2}, x])/(3*(1 + x)) + (96*HPL[{-1, 4}, x])/(1 + x) + (10*Pi^2*HPL[{2, 0}, x])/(3*(1 + x)) - (20*Pi^2*HPL[{2, 1}, x])/(1 + x) - (16*HPL[{2, 3}, x])/(1 + x) + (24*HPL[{3, 2}, x])/(1 + x) - (14*HPL[{4, 0}, x])/(1 + x) - (96*HPL[{4, 1}, x])/(1 + x) + (64*Pi^2*HPL[{-1, -1, -1}, x])/(1 + x) - (128*Pi^2*HPL[{-1, -1, 0}, x])/(3*(1 + x)) + (128*HPL[{-1, -1, 3}, x])/(1 + x) - (28*Pi^2*HPL[{-1, 0, 0}, x])/(3*(1 + x)) - (24*HPL[{-1, 3, 0}, x])/(1 + x) - (128*HPL[{-1, 3, 1}, x])/(1 + x) - (7*Pi^2*HPL[{0, 0, 0}, x])/(3*(1 + x)) + (16*HPL[{2, 1, 2}, x])/(1 + x) + (4*HPL[{2, 2, 0}, x])/(1 + x) - (2*HPL[{3, 0, 0}, x])/(1 + x) + (20*HPL[{3, 1, 0}, x])/(1 + x) + (128*HPL[{3, 1, 1}, x])/(1 + x) - (256*HPL[{-1, -1, -1, 2}, x])/(1 + x) + (256*HPL[{-1, -1, 2, 1}, x])/(1 + x) + (16*HPL[{-1, 2, 0, 0}, x])/(1 + x) - (16*HPL[{-1, 2, 1, 0}, x])/(1 + x) - (256*HPL[{-1, 2, 1, 1}, x])/(1 + x) + (8*HPL[{2, 0, 0, 0}, x])/(1 + x) - (12*HPL[{2, 1, 0, 0}, x])/(1 + x) + (24*HPL[{2, 1, 1, 0}, x])/(1 + x) + (256*HPL[{2, 1, 1, 1}, x])/(1 + x) + (256*HPL[{-1, -1, -1, -1, 0}, x])/(1 + x) - (128*HPL[{-1, -1, -1, 0, 0}, x])/(1 + x) - (96*HPL[{-1, -1, 0, 0, 0}, x])/(1 + x) - (48*HPL[{-1, 0, 0, 0, 0}, x])/(1 + x) - (22*HPL[{0, 0, 0, 0, 0}, x])/(1 + x) + (52*Pi^2*Zeta[3])/(1 + x) - (36*HPL[{2}, x]*Zeta[3])/(1 + x) - (480*HPL[{-1, -1}, x]*Zeta[3])/(1 + x) + (64*HPL[{-1, 0}, x]*Zeta[3])/(1 + x) + (10*HPL[{0, 0}, x]*Zeta[3])/(1 + x) - (154*Zeta[5])/(1 + x))
  + ALARM*ep^3 );
U11[x -> 1,ep] = (1-x)^(1 - 4*ep) * (
  - 1/(2*ep^2) - 2/ep - 8 + (2*Pi^2)/3 + ep*(-32 + (8*Pi^2)/3 + 10*Zeta[3]) + ep^2*(-128 + (32*Pi^2)/3 - (4*Pi^4)/45 + 40*Zeta[3])
  + ALARM*ep^3 );

U12d[x,ep] = (
  + 2/(ep^3*(-1 + x))
  + ((2*HPL[{0}, x])/(-1 + x) + (8*HPL[{1}, x])/(-1 + x))/ep^2
  + ((-7*Pi^2)/(3*(-1 + x)) + (6*HPL[{2}, x])/(-1 + x) + (8*HPL[{1, 0}, x])/(-1 + x) + (32*HPL[{1, 1}, x])/(-1 + x))/ep
  - (2*Pi^2*HPL[{0}, x])/(3*(-1 + x)) - (28*Pi^2*HPL[{1}, x])/(3*(-1 + x)) - (4*HPL[{3}, x])/(-1 + x) + (24*HPL[{1, 2}, x])/(-1 + x) + (18*HPL[{2, 1}, x])/(-1 + x) - (4*HPL[{0, 0, 0}, x])/(-1 + x) + (32*HPL[{1, 1, 0}, x])/(-1 + x) + (128*HPL[{1, 1, 1}, x])/(-1 + x) - (34*Zeta[3])/(-1 + x)
  + ep*((-7*Pi^4)/(90*(-1 + x)) - (4*Pi^2*HPL[{2}, x])/(-1 + x) - (8*Pi^2*HPL[{0, 0}, x])/(3*(-1 + x)) - (8*Pi^2*HPL[{1, 0}, x])/(3*(-1 + x)) - (112*Pi^2*HPL[{1, 1}, x])/(3*(-1 + x)) - (16*HPL[{1, 3}, x])/(-1 + x) - (8*HPL[{3, 0}, x])/(-1 + x) - (20*HPL[{3, 1}, x])/(-1 + x) + (96*HPL[{1, 1, 2}, x])/(-1 + x) + (72*HPL[{1, 2, 1}, x])/(-1 + x) + (54*HPL[{2, 1, 1}, x])/(-1 + x) + (12*HPL[{0, 0, 0, 0}, x])/(-1 + x) - (16*HPL[{1, 0, 0, 0}, x])/(-1 + x) + (128*HPL[{1, 1, 1, 0}, x])/(-1 + x) + (512*HPL[{1, 1, 1, 1}, x])/(-1 + x) - (12*HPL[{0}, x]*Zeta[3])/(-1 + x) - (136*HPL[{1}, x]*Zeta[3])/(-1 + x))
  + ep^2*((17*Pi^4*HPL[{0}, x])/(15*(-1 + x)) - (14*Pi^4*HPL[{1}, x])/(45*(-1 + x)) + (8*HPL[{5}, x])/(-1 + x) - (16*Pi^2*HPL[{1, 2}, x])/(-1 + x) - (12*Pi^2*HPL[{2, 1}, x])/(-1 + x) - (16*HPL[{3, 2}, x])/(-1 + x) + (24*HPL[{4, 0}, x])/(-1 + x) + (24*HPL[{4, 1}, x])/(-1 + x) + (28*Pi^2*HPL[{0, 0, 0}, x])/(3*(-1 + x)) - (32*Pi^2*HPL[{1, 0, 0}, x])/(3*(-1 + x)) - (32*Pi^2*HPL[{1, 1, 0}, x])/(3*(-1 + x)) - (448*Pi^2*HPL[{1, 1, 1}, x])/(3*(-1 + x)) - (64*HPL[{1, 1, 3}, x])/(-1 + x) - (32*HPL[{1, 3, 0}, x])/(-1 + x) - (80*HPL[{1, 3, 1}, x])/(-1 + x) + (8*HPL[{3, 0, 0}, x])/(-1 + x) - (16*HPL[{3, 1, 0}, x])/(-1 + x) - (76*HPL[{3, 1, 1}, x])/(-1 + x) + (384*HPL[{1, 1, 1, 2}, x])/(-1 + x) + (288*HPL[{1, 1, 2, 1}, x])/(-1 + x) + (216*HPL[{1, 2, 1, 1}, x])/(-1 + x) + (162*HPL[{2, 1, 1, 1}, x])/(-1 + x) - (28*HPL[{0, 0, 0, 0, 0}, x])/(-1 + x) + (48*HPL[{1, 0, 0, 0, 0}, x])/(-1 + x) - (64*HPL[{1, 1, 0, 0, 0}, x])/(-1 + x) + (512*HPL[{1, 1, 1, 1, 0}, x])/(-1 + x) + (2048*HPL[{1, 1, 1, 1, 1}, x])/(-1 + x) + (130*Pi^2*Zeta[3])/(3*(-1 + x)) - (48*HPL[{2}, x]*Zeta[3])/(-1 + x) + (8*HPL[{0, 0}, x]*Zeta[3])/(-1 + x) - (48*HPL[{1, 0}, x]*Zeta[3])/(-1 + x) - (544*HPL[{1, 1}, x]*Zeta[3])/(-1 + x) - (198*Zeta[5])/(-1 + x))
  + ALARM*ep^3 );
U12d[x -> 1,ep] = (1-x)^(-1 - 4*ep) * (
  - 2/ep^3 + (8*Pi^2)/(3*ep) - (13*ep*Pi^4)/45 + 36*Zeta[3] + ep^2*(-48*Pi^2*Zeta[3] + 396*Zeta[5])
  + ALARM*ep^3 );

U12t[x,ep] = (
  - 3/(ep^3*(-1 + x))
  + (-(HPL[{0}, x]/(-1 + x)) - (12*HPL[{1}, x])/(-1 + x))/ep^2
  + ((13*Pi^2)/(3*(-1 + x)) - (4*HPL[{2}, x])/(-1 + x) + (2*HPL[{0, 0}, x])/(-1 + x) - (4*HPL[{1, 0}, x])/(-1 + x) - (48*HPL[{1, 1}, x])/(-1 + x))/ep
  - (Pi^2*HPL[{0}, x])/(3*(-1 + x)) + (52*Pi^2*HPL[{1}, x])/(3*(-1 + x)) + (8*HPL[{3}, x])/(-1 + x) - (16*HPL[{1, 2}, x])/(-1 + x) - (2*HPL[{2, 0}, x])/(-1 + x) - (16*HPL[{2, 1}, x])/(-1 + x) - (4*HPL[{0, 0, 0}, x])/(-1 + x) + (8*HPL[{1, 0, 0}, x])/(-1 + x) - (16*HPL[{1, 1, 0}, x])/(-1 + x) - (192*HPL[{1, 1, 1}, x])/(-1 + x) + (50*Zeta[3])/(-1 + x)
 + ep*((-101*Pi^4)/(90*(-1 + x)) + (16*Pi^2*HPL[{2}, x])/(3*(-1 + x)) - (16*HPL[{4}, x])/(-1 + x) + (2*Pi^2*HPL[{0, 0}, x])/(3*(-1 + x)) - (4*Pi^2*HPL[{1, 0}, x])/(3*(-1 + x)) + (208*Pi^2*HPL[{1, 1}, x])/(3*(-1 + x)) + (32*HPL[{1, 3}, x])/(-1 + x) - (8*HPL[{2, 2}, x])/(-1 + x) + (4*HPL[{3, 0}, x])/(-1 + x) + (32*HPL[{3, 1}, x])/(-1 + x) - (64*HPL[{1, 1, 2}, x])/(-1 + x) - (8*HPL[{1, 2, 0}, x])/(-1 + x) - (64*HPL[{1, 2, 1}, x])/(-1 + x) + (2*HPL[{2, 0, 0}, x])/(-1 + x) - (8*HPL[{2, 1, 0}, x])/(-1 + x) - (64*HPL[{2, 1, 1}, x])/(-1 + x) + (8*HPL[{0, 0, 0, 0}, x])/(-1 + x) - (16*HPL[{1, 0, 0, 0}, x])/(-1 + x) + (32*HPL[{1, 1, 0, 0}, x])/(-1 + x) - (64*HPL[{1, 1, 1, 0}, x])/(-1 + x) - (768*HPL[{1, 1, 1, 1}, x])/(-1 + x) + (2*HPL[{0}, x]*Zeta[3])/(-1 + x) + (200*HPL[{1}, x]*Zeta[3])/(-1 + x))
 + ep^2*(-(Pi^4*HPL[{0}, x])/(10*(-1 + x)) - (202*Pi^4*HPL[{1}, x])/(45*(-1 + x)) - (32*Pi^2*HPL[{3}, x])/(3*(-1 + x)) + (32*HPL[{5}, x])/(-1 + x) + (64*Pi^2*HPL[{1, 2}, x])/(3*(-1 + x)) - (64*HPL[{1, 4}, x])/(-1 + x) + (2*Pi^2*HPL[{2, 0}, x])/(3*(-1 + x)) + (64*Pi^2*HPL[{2, 1}, x])/(3*(-1 + x)) + (8*HPL[{2, 3}, x])/(-1 + x) + (16*HPL[{3, 2}, x])/(-1 + x) - (8*HPL[{4, 0}, x])/(-1 + x) - (64*HPL[{4, 1}, x])/(-1 + x) - (4*Pi^2*HPL[{0, 0, 0}, x])/(3*(-1 + x)) + (8*Pi^2*HPL[{1, 0, 0}, x])/(3*(-1 + x)) - (16*Pi^2*HPL[{1, 1, 0}, x])/(3*(-1 + x)) + (832*Pi^2*HPL[{1, 1, 1}, x])/(3*(-1 + x)) + (128*HPL[{1, 1, 3}, x])/(-1 + x) - (32*HPL[{1, 2, 2}, x])/(-1 + x) + (16*HPL[{1, 3, 0}, x])/(-1 + x) + (128*HPL[{1, 3, 1}, x])/(-1 + x) - (32*HPL[{2, 1, 2}, x])/(-1 + x) - (4*HPL[{2, 2, 0}, x])/(-1 + x) - (32*HPL[{2, 2, 1}, x])/(-1 + x) - (4*HPL[{3, 0, 0}, x])/(-1 + x) + (16*HPL[{3, 1, 0}, x])/(-1 + x) + (128*HPL[{3, 1, 1}, x])/(-1 + x) - (256*HPL[{1, 1, 1, 2}, x])/(-1 + x) - (32*HPL[{1, 1, 2, 0}, x])/(-1 + x) - (256*HPL[{1, 1, 2, 1}, x])/(-1 + x) + (8*HPL[{1, 2, 0, 0}, x])/(-1 + x) - (32*HPL[{1, 2, 1, 0}, x])/(-1 + x) - (256*HPL[{1, 2, 1, 1}, x])/(-1 + x) - (2*HPL[{2, 0, 0, 0}, x])/(-1 + x) + (8*HPL[{2, 1, 0, 0}, x])/(-1 + x) - (32*HPL[{2, 1, 1, 0}, x])/(-1 + x) - (256*HPL[{2, 1, 1, 1}, x])/(-1 + x) - (16*HPL[{0, 0, 0, 0, 0}, x])/(-1 + x) + (32*HPL[{1, 0, 0, 0, 0}, x])/(-1 + x) - (64*HPL[{1, 1, 0, 0, 0}, x])/(-1 + x) + (128*HPL[{1, 1, 1, 0, 0}, x])/(-1 + x) - (256*HPL[{1, 1, 1, 1, 0}, x])/(-1 + x) - (3072*HPL[{1, 1, 1, 1, 1}, x])/(-1 + x) - (190*Pi^2*Zeta[3])/(3*(-1 + x)) + (56*HPL[{2}, x]*Zeta[3])/(-1 + x) - (4*HPL[{0, 0}, x]*Zeta[3])/(-1 + x) + (8*HPL[{1, 0}, x]*Zeta[3])/(-1 + x) + (800*HPL[{1, 1}, x]*Zeta[3])/(-1 + x) + (234*Zeta[5])/(-1 + x))
 + ALARM*ep^3 );
U12t[x -> 1,ep] = (1-x)^(-1 - 4*ep) * (
 + 3/ep^3 - (13*Pi^2)/(3*ep) + (19*ep*Pi^4)/30 - 70*Zeta[3] + ep^2*((298*Pi^2*Zeta[3])/3 - 686*Zeta[5])
 + ALARM*ep^3 );

U13[x,ep] = (
  + 4/(ep^3*(-1 + x))
  + ((2*HPL[{0}, x])/(-1 + x) + (16*HPL[{1}, x])/(-1 + x))/ep^2
  + ((-16*Pi^2)/(3*(-1 + x)) + (8*HPL[{2}, x])/(-1 + x) - (4*HPL[{0, 0}, x])/(-1 + x) + (8*HPL[{1, 0}, x])/(-1 + x) + (64*HPL[{1, 1}, x])/(-1 + x))/ep
  + (2*Pi^2*HPL[{0}, x])/(3*(-1 + x)) - (64*Pi^2*HPL[{1}, x])/(3*(-1 + x)) - (16*HPL[{3}, x])/(-1 + x) + (32*HPL[{1, 2}, x])/(-1 + x) + (4*HPL[{2, 0}, x])/(-1 + x) + (32*HPL[{2, 1}, x])/(-1 + x) + (8*HPL[{0, 0, 0}, x])/(-1 + x) - (12*HPL[{1, 0, 0}, x])/(-1 + x) + (32*HPL[{1, 1, 0}, x])/(-1 + x) + (256*HPL[{1, 1, 1}, x])/(-1 + x) - (44*Zeta[3])/(-1 + x)
  + ep*((7*Pi^4)/(5*(-1 + x)) - (32*Pi^2*HPL[{2}, x])/(3*(-1 + x)) + (32*HPL[{4}, x])/(-1 + x) - (4*Pi^2*HPL[{0, 0}, x])/(3*(-1 + x)) - (256*Pi^2*HPL[{1, 1}, x])/(3*(-1 + x)) - (48*HPL[{1, 3}, x])/(-1 + x) + (16*HPL[{2, 2}, x])/(-1 + x) - (8*HPL[{3, 0}, x])/(-1 + x) - (64*HPL[{3, 1}, x])/(-1 + x) + (128*HPL[{1, 1, 2}, x])/(-1 + x) + (16*HPL[{1, 2, 0}, x])/(-1 + x) + (128*HPL[{1, 2, 1}, x])/(-1 + x) - (4*HPL[{2, 0, 0}, x])/(-1 + x) + (16*HPL[{2, 1, 0}, x])/(-1 + x) + (128*HPL[{2, 1, 1}, x])/(-1 + x) - (16*HPL[{0, 0, 0, 0}, x])/(-1 + x) + (20*HPL[{1, 0, 0, 0}, x])/(-1 + x) - (40*HPL[{1, 1, 0, 0}, x])/(-1 + x) + (128*HPL[{1, 1, 1, 0}, x])/(-1 + x) + (1024*HPL[{1, 1, 1, 1}, x])/(-1 + x) - (4*HPL[{0}, x]*Zeta[3])/(-1 + x) - (200*HPL[{1}, x]*Zeta[3])/(-1 + x))
  + ep^2*((Pi^4*HPL[{0}, x])/(5*(-1 + x)) + (50*Pi^4*HPL[{1}, x])/(9*(-1 + x)) + (64*Pi^2*HPL[{3}, x])/(3*(-1 + x)) - (64*HPL[{5}, x])/(-1 + x) - (128*Pi^2*HPL[{1, 2}, x])/(3*(-1 + x)) + (80*HPL[{1, 4}, x])/(-1 + x) - (4*Pi^2*HPL[{2, 0}, x])/(3*(-1 + x)) - (128*Pi^2*HPL[{2, 1}, x])/(3*(-1 + x)) - (16*HPL[{2, 3}, x])/(-1 + x) - (32*HPL[{3, 2}, x])/(-1 + x) + (16*HPL[{4, 0}, x])/(-1 + x) + (128*HPL[{4, 1}, x])/(-1 + x) + (8*Pi^2*HPL[{0, 0, 0}, x])/(3*(-1 + x)) - (4*Pi^2*HPL[{1, 0, 0}, x])/(3*(-1 + x)) - (16*Pi^2*HPL[{1, 1, 0}, x])/(3*(-1 + x)) - (1024*Pi^2*HPL[{1, 1, 1}, x])/(3*(-1 + x)) - (160*HPL[{1, 1, 3}, x])/(-1 + x) + (64*HPL[{1, 2, 2}, x])/(-1 + x) - (24*HPL[{1, 3, 0}, x])/(-1 + x) - (192*HPL[{1, 3, 1}, x])/(-1 + x) + (64*HPL[{2, 1, 2}, x])/(-1 + x) + (8*HPL[{2, 2, 0}, x])/(-1 + x) + (64*HPL[{2, 2, 1}, x])/(-1 + x) + (8*HPL[{3, 0, 0}, x])/(-1 + x) - (32*HPL[{3, 1, 0}, x])/(-1 + x) - (256*HPL[{3, 1, 1}, x])/(-1 + x) + (512*HPL[{1, 1, 1, 2}, x])/(-1 + x) + (64*HPL[{1, 1, 2, 0}, x])/(-1 + x) + (512*HPL[{1, 1, 2, 1}, x])/(-1 + x) - (16*HPL[{1, 2, 0, 0}, x])/(-1 + x) + (64*HPL[{1, 2, 1, 0}, x])/(-1 + x) + (512*HPL[{1, 2, 1, 1}, x])/(-1 + x) + (4*HPL[{2, 0, 0, 0}, x])/(-1 + x) - (16*HPL[{2, 1, 0, 0}, x])/(-1 + x) + (64*HPL[{2, 1, 1, 0}, x])/(-1 + x) + (512*HPL[{2, 1, 1, 1}, x])/(-1 + x) + (32*HPL[{0, 0, 0, 0, 0}, x])/(-1 + x) - (36*HPL[{1, 0, 0, 0, 0}, x])/(-1 + x) + (56*HPL[{1, 1, 0, 0, 0}, x])/(-1 + x) - (144*HPL[{1, 1, 1, 0, 0}, x])/(-1 + x) + (512*HPL[{1, 1, 1, 1, 0}, x])/(-1 + x) + (4096*HPL[{1, 1, 1, 1, 1}, x])/(-1 + x) + (140*Pi^2*Zeta[3])/(3*(-1 + x)) - (112*HPL[{2}, x]*Zeta[3])/(-1 + x) + (8*HPL[{0, 0}, x]*Zeta[3])/(-1 + x) - (16*HPL[{1, 0}, x]*Zeta[3])/(-1 + x) - (848*HPL[{1, 1}, x]*Zeta[3])/(-1 + x) - (120*Zeta[5])/(-1 + x))
  + ALARM*ep^3 );
U13[x -> 1,ep] = (1-x)^(-1 - 4*ep) * (
  - 4/ep^3 + (16*Pi^2)/(3*ep) - (32*ep*Pi^4)/45 + 80*Zeta[3] + ep^2*((-320*Pi^2*Zeta[3])/3 + 816*Zeta[5])
  + ALARM*ep^3 );

U14[x,ep] = (
  + ((-1 + x)^(-1) - x^(-1))/ep^3
  + (-(HPL[{0}, x]/(-1 + x)) + HPL[{0}, x]/x + (4*HPL[{1}, x])/(-1 + x) - (4*HPL[{1}, x])/x)/ep^2
  + ((-2*Pi^2)/(3*(-1 + x)) + (2*Pi^2)/(3*x) - (4*HPL[{2}, x])/(-1 + x) + (4*HPL[{2}, x])/x + (16*HPL[{1, 1}, x])/(-1 + x) - (16*HPL[{1, 1}, x])/x)/ep
  + (5*Pi^2*HPL[{0}, x])/(3*(-1 + x)) - (5*Pi^2*HPL[{0}, x])/(3*x) - (13*Pi^2*HPL[{1}, x])/(3*(-1 + x)) + (13*Pi^2*HPL[{1}, x])/(3*x) + (2*HPL[{2, 0}, x])/(-1 + x) - (2*HPL[{2, 0}, x])/x - (16*HPL[{2, 1}, x])/(-1 + x) + (16*HPL[{2, 1}, x])/x + (4*HPL[{0, 0, 0}, x])/(-1 + x) - (4*HPL[{0, 0, 0}, x])/x + (6*HPL[{1, 1, 0}, x])/(-1 + x) - (6*HPL[{1, 1, 0}, x])/x + (64*HPL[{1, 1, 1}, x])/(-1 + x) - (64*HPL[{1, 1, 1}, x])/x - (6*Zeta[3])/(-1 + x) + (6*Zeta[3])/x
  + ep*(-Pi^4/(2*(-1 + x)) + Pi^4/(2*x) + (10*Pi^2*HPL[{2}, x])/(3*(-1 + x)) - (10*Pi^2*HPL[{2}, x])/(3*x) + (16*HPL[{4}, x])/(-1 + x) - (16*HPL[{4}, x])/x - (4*Pi^2*HPL[{0, 0}, x])/(-1 + x) + (4*Pi^2*HPL[{0, 0}, x])/x + (Pi^2*HPL[{1, 0}, x])/(-1 + x) - (Pi^2*HPL[{1, 0}, x])/x - (20*Pi^2*HPL[{1, 1}, x])/(-1 + x) + (20*Pi^2*HPL[{1, 1}, x])/x + (8*HPL[{2, 2}, x])/(-1 + x) - (8*HPL[{2, 2}, x])/x - (8*HPL[{3, 0}, x])/(-1 + x) + (8*HPL[{3, 0}, x])/x + (24*HPL[{1, 1, 2}, x])/(-1 + x) - (24*HPL[{1, 1, 2}, x])/x + (6*HPL[{1, 2, 0}, x])/(-1 + x) - (6*HPL[{1, 2, 0}, x])/x - (2*HPL[{2, 0, 0}, x])/(-1 + x) + (2*HPL[{2, 0, 0}, x])/x - (4*HPL[{2, 1, 0}, x])/(-1 + x) + (4*HPL[{2, 1, 0}, x])/x - (64*HPL[{2, 1, 1}, x])/(-1 + x) + (64*HPL[{2, 1, 1}, x])/x - (16*HPL[{0, 0, 0, 0}, x])/(-1 + x) + (16*HPL[{0, 0, 0, 0}, x])/x + (6*HPL[{1, 0, 0, 0}, x])/(-1 + x) - (6*HPL[{1, 0, 0, 0}, x])/x - (2*HPL[{1, 1, 0, 0}, x])/(-1 + x) + (2*HPL[{1, 1, 0, 0}, x])/x + (32*HPL[{1, 1, 1, 0}, x])/(-1 + x) - (32*HPL[{1, 1, 1, 0}, x])/x + (256*HPL[{1, 1, 1, 1}, x])/(-1 + x) - (256*HPL[{1, 1, 1, 1}, x])/x + (10*HPL[{0}, x]*Zeta[3])/(-1 + x) - (10*HPL[{0}, x]*Zeta[3])/x - (50*HPL[{1}, x]*Zeta[3])/(-1 + x) + (50*HPL[{1}, x]*Zeta[3])/x)
  + ep^2*((11*Pi^4*HPL[{0}, x])/(10*(-1 + x)) - (11*Pi^4*HPL[{0}, x])/(10*x) + (8*Pi^4*HPL[{1}, x])/(45*(-1 + x)) - (8*Pi^4*HPL[{1}, x])/(45*x) + (4*Pi^2*HPL[{3}, x])/(-1 + x) - (4*Pi^2*HPL[{3}, x])/x - (64*HPL[{5}, x])/(-1 + x) + (64*HPL[{5}, x])/x - (10*Pi^2*HPL[{1, 2}, x])/(3*(-1 + x)) + (10*Pi^2*HPL[{1, 2}, x])/(3*x) + (24*HPL[{1, 4}, x])/(-1 + x) - (24*HPL[{1, 4}, x])/x - (8*Pi^2*HPL[{2, 0}, x])/(3*(-1 + x)) + (8*Pi^2*HPL[{2, 0}, x])/(3*x) + (56*Pi^2*HPL[{2, 1}, x])/(3*(-1 + x)) - (56*Pi^2*HPL[{2, 1}, x])/(3*x) - (8*HPL[{2, 3}, x])/(-1 + x) + (8*HPL[{2, 3}, x])/x - (32*HPL[{3, 2}, x])/(-1 + x) + (32*HPL[{3, 2}, x])/x + (24*HPL[{4, 0}, x])/(-1 + x) - (24*HPL[{4, 0}, x])/x + (64*HPL[{4, 1}, x])/(-1 + x) - (64*HPL[{4, 1}, x])/x + (28*Pi^2*HPL[{0, 0, 0}, x])/(3*(-1 + x)) - (28*Pi^2*HPL[{0, 0, 0}, x])/(3*x) - (19*Pi^2*HPL[{1, 0, 0}, x])/(3*(-1 + x)) + (19*Pi^2*HPL[{1, 0, 0}, x])/(3*x) - (10*Pi^2*HPL[{1, 1, 0}, x])/(3*(-1 + x)) + (10*Pi^2*HPL[{1, 1, 0}, x])/(3*x) - (84*Pi^2*HPL[{1, 1, 1}, x])/(-1 + x) + (84*Pi^2*HPL[{1, 1, 1}, x])/x - (8*HPL[{1, 1, 3}, x])/(-1 + x) + (8*HPL[{1, 1, 3}, x])/x + (24*HPL[{1, 2, 2}, x])/(-1 + x) - (24*HPL[{1, 2, 2}, x])/x - (14*HPL[{1, 3, 0}, x])/(-1 + x) + (14*HPL[{1, 3, 0}, x])/x - (16*HPL[{2, 1, 2}, x])/(-1 + x) + (16*HPL[{2, 1, 2}, x])/x - (8*HPL[{2, 2, 0}, x])/(-1 + x) + (8*HPL[{2, 2, 0}, x])/x + (32*HPL[{2, 2, 1}, x])/(-1 + x) - (32*HPL[{2, 2, 1}, x])/x + (8*HPL[{3, 0, 0}, x])/(-1 + x) - (8*HPL[{3, 0, 0}, x])/x - (8*HPL[{3, 1, 0}, x])/(-1 + x) + (8*HPL[{3, 1, 0}, x])/x + (128*HPL[{1, 1, 1, 2}, x])/(-1 + x) - (128*HPL[{1, 1, 1, 2}, x])/x + (20*HPL[{1, 1, 2, 0}, x])/(-1 + x) - (20*HPL[{1, 1, 2, 0}, x])/x + (96*HPL[{1, 1, 2, 1}, x])/(-1 + x) - (96*HPL[{1, 1, 2, 1}, x])/x - (10*HPL[{1, 2, 0, 0}, x])/(-1 + x) + (10*HPL[{1, 2, 0, 0}, x])/x + (4*HPL[{1, 2, 1, 0}, x])/(-1 + x) - (4*HPL[{1, 2, 1, 0}, x])/x - (10*HPL[{2, 0, 0, 0}, x])/(-1 + x) + (10*HPL[{2, 0, 0, 0}, x])/x - (4*HPL[{2, 1, 0, 0}, x])/(-1 + x) + (4*HPL[{2, 1, 0, 0}, x])/x - (32*HPL[{2, 1, 1, 0}, x])/(-1 + x) + (32*HPL[{2, 1, 1, 0}, x])/x - (256*HPL[{2, 1, 1, 1}, x])/(-1 + x) + (256*HPL[{2, 1, 1, 1}, x])/x + (48*HPL[{0, 0, 0, 0, 0}, x])/(-1 + x) - (48*HPL[{0, 0, 0, 0, 0}, x])/x - (28*HPL[{1, 0, 0, 0, 0}, x])/(-1 + x) + (28*HPL[{1, 0, 0, 0, 0}, x])/x + (10*HPL[{1, 1, 0, 0, 0}, x])/(-1 + x) - (10*HPL[{1, 1, 0, 0, 0}, x])/x - (16*HPL[{1, 1, 1, 0, 0}, x])/(-1 + x) + (16*HPL[{1, 1, 1, 0, 0}, x])/x + (136*HPL[{1, 1, 1, 1, 0}, x])/(-1 + x) - (136*HPL[{1, 1, 1, 1, 0}, x])/x + (1024*HPL[{1, 1, 1, 1, 1}, x])/(-1 + x) - (1024*HPL[{1, 1, 1, 1, 1}, x])/x + (44*HPL[{2}, x]*Zeta[3])/(-1 + x) - (44*HPL[{2}, x]*Zeta[3])/x - (16*HPL[{0, 0}, x]*Zeta[3])/(-1 + x) + (16*HPL[{0, 0}, x]*Zeta[3])/x + (2*HPL[{1, 0}, x]*Zeta[3])/(-1 + x) - (2*HPL[{1, 0}, x]*Zeta[3])/x - (232*HPL[{1, 1}, x]*Zeta[3])/(-1 + x) + (232*HPL[{1, 1}, x]*Zeta[3])/x - (22*Zeta[5])/(-1 + x) + (22*Zeta[5])/x)
  + ALARM*ep^3 );
U14[x -> 1,ep] = (1-x)^(-1 - 4*ep) * (
  - ep^(-3) + (4*Pi^2)/(3*ep) - (8*ep*Pi^4)/45 + 20*Zeta[3] + ep^2*((-80*Pi^2*Zeta[3])/3 + 204*Zeta[5])
  + ALARM*ep^3 );

U15[x,ep] = (
  + 1/(ep^3*(-1 + x))
  + ((2*HPL[{0}, x])/(-1 + x) + (4*HPL[{1}, x])/(-1 + x))/ep^2
  + ((-2*Pi^2)/(-1 + x) + (8*HPL[{2}, x])/(-1 + x) + (4*HPL[{1, 0}, x])/(-1 + x) + (16*HPL[{1, 1}, x])/(-1 + x))/ep
  - (7*Pi^2*HPL[{0}, x])/(3*(-1 + x)) - (19*Pi^2*HPL[{1}, x])/(3*(-1 + x)) + (16*HPL[{1, 2}, x])/(-1 + x) + (2*HPL[{2, 0}, x])/(-1 + x) + (32*HPL[{2, 1}, x])/(-1 + x) - (2*HPL[{0, 0, 0}, x])/(-1 + x) - (2*HPL[{1, 0, 0}, x])/(-1 + x) + (10*HPL[{1, 1, 0}, x])/(-1 + x) + (64*HPL[{1, 1, 1}, x])/(-1 + x) - (24*Zeta[3])/(-1 + x)
  + ep*((23*Pi^4)/(18*(-1 + x)) - (10*Pi^2*HPL[{2}, x])/(-1 + x) - (8*HPL[{4}, x])/(-1 + x) + (5*Pi^2*HPL[{0, 0}, x])/(3*(-1 + x)) - (3*Pi^2*HPL[{1, 0}, x])/(-1 + x) - (68*Pi^2*HPL[{1, 1}, x])/(3*(-1 + x)) - (8*HPL[{1, 3}, x])/(-1 + x) + (8*HPL[{2, 2}, x])/(-1 + x) + (2*HPL[{3, 0}, x])/(-1 + x) + (40*HPL[{1, 1, 2}, x])/(-1 + x) + (6*HPL[{1, 2, 0}, x])/(-1 + x) + (64*HPL[{1, 2, 1}, x])/(-1 + x) - (6*HPL[{2, 0, 0}, x])/(-1 + x) + (12*HPL[{2, 1, 0}, x])/(-1 + x) + (128*HPL[{2, 1, 1}, x])/(-1 + x) + (4*HPL[{0, 0, 0, 0}, x])/(-1 + x) - (10*HPL[{1, 1, 0, 0}, x])/(-1 + x) + (32*HPL[{1, 1, 1, 0}, x])/(-1 + x) + (256*HPL[{1, 1, 1, 1}, x])/(-1 + x) - (18*HPL[{0}, x]*Zeta[3])/(-1 + x) - (66*HPL[{1}, x]*Zeta[3])/(-1 + x))
  + ep^2*((Pi^4*HPL[{0}, x])/(5*(-1 + x)) + (124*Pi^4*HPL[{1}, x])/(45*(-1 + x)) - (2*Pi^2*HPL[{3}, x])/(3*(-1 + x)) + (16*HPL[{5}, x])/(-1 + x) - (62*Pi^2*HPL[{1, 2}, x])/(3*(-1 + x)) + (8*Pi^2*HPL[{2, 0}, x])/(3*(-1 + x)) - (124*Pi^2*HPL[{2, 1}, x])/(3*(-1 + x)) - (24*HPL[{2, 3}, x])/(-1 + x) + (8*HPL[{3, 2}, x])/(-1 + x) - (6*HPL[{4, 0}, x])/(-1 + x) - (32*HPL[{4, 1}, x])/(-1 + x) - (Pi^2*HPL[{0, 0, 0}, x])/(-1 + x) + (7*Pi^2*HPL[{1, 0, 0}, x])/(3*(-1 + x)) - (10*Pi^2*HPL[{1, 1, 0}, x])/(3*(-1 + x)) - (260*Pi^2*HPL[{1, 1, 1}, x])/(3*(-1 + x)) - (40*HPL[{1, 1, 3}, x])/(-1 + x) + (24*HPL[{1, 2, 2}, x])/(-1 + x) - (2*HPL[{1, 3, 0}, x])/(-1 + x) - (32*HPL[{1, 3, 1}, x])/(-1 + x) + (48*HPL[{2, 1, 2}, x])/(-1 + x) + (8*HPL[{2, 2, 0}, x])/(-1 + x) + (32*HPL[{2, 2, 1}, x])/(-1 + x) + (2*HPL[{3, 0, 0}, x])/(-1 + x) + (4*HPL[{3, 1, 0}, x])/(-1 + x) + (128*HPL[{1, 1, 1, 2}, x])/(-1 + x) + (20*HPL[{1, 1, 2, 0}, x])/(-1 + x) + (160*HPL[{1, 1, 2, 1}, x])/(-1 + x) - (10*HPL[{1, 2, 0, 0}, x])/(-1 + x) + (28*HPL[{1, 2, 1, 0}, x])/(-1 + x) + (256*HPL[{1, 2, 1, 1}, x])/(-1 + x) + (10*HPL[{2, 0, 0, 0}, x])/(-1 + x) - (20*HPL[{2, 1, 0, 0}, x])/(-1 + x) + (56*HPL[{2, 1, 1, 0}, x])/(-1 + x) + (512*HPL[{2, 1, 1, 1}, x])/(-1 + x) - (6*HPL[{0, 0, 0, 0, 0}, x])/(-1 + x) + (2*HPL[{1, 0, 0, 0, 0}, x])/(-1 + x) + (10*HPL[{1, 1, 0, 0, 0}, x])/(-1 + x) - (40*HPL[{1, 1, 1, 0, 0}, x])/(-1 + x) + (120*HPL[{1, 1, 1, 1, 0}, x])/(-1 + x) + (1024*HPL[{1, 1, 1, 1, 1}, x])/(-1 + x) + (106*Pi^2*Zeta[3])/(3*(-1 + x)) - (92*HPL[{2}, x]*Zeta[3])/(-1 + x) + (14*HPL[{0, 0}, x]*Zeta[3])/(-1 + x) - (22*HPL[{1, 0}, x]*Zeta[3])/(-1 + x) - (224*HPL[{1, 1}, x]*Zeta[3])/(-1 + x) - (40*Zeta[5])/(-1 + x))
  + ALARM*ep^3 );
U15[x -> 1,ep] = (1-x)^(-1 - 4*ep) * (
  - ep^(-3) + (4*Pi^2)/(3*ep) - (8*ep*Pi^4)/45 + 20*Zeta[3] + ep^2*((-80*Pi^2*Zeta[3])/3 + 204*Zeta[5])
  + ALARM*ep^3 );

U16[x,ep] = (
  - (2*HPL[{0}, x])/(ep^2*x)
  + ((2*Pi^2)/(3*x) - (8*HPL[{2}, x])/x - (8*HPL[{-1, 0}, x])/x + (8*HPL[{0, 0}, x])/x)/ep
  + (6*Pi^2*HPL[{-1}, x])/x - (2*Pi^2*HPL[{0}, x])/x + (32*HPL[{3}, x])/x + (16*HPL[{-2, 0}, x])/x - (32*HPL[{-1, 2}, x])/x - (4*HPL[{2, 0}, x])/x - (32*HPL[{2, 1}, x])/x + (8*HPL[{-1, -1, 0}, x])/x + (4*HPL[{-1, 0, 0}, x])/x - (24*HPL[{0, 0, 0}, x])/x - (8*Zeta[3])/x
  + ep*((2*Pi^4)/(5*x) - (12*Pi^2*HPL[{-2}, x])/x + (32*Pi^2*HPL[{2}, x])/(3*x) - (96*HPL[{4}, x])/x - (32*HPL[{-3, 0}, x])/x + (64*HPL[{-2, 2}, x])/x - (6*Pi^2*HPL[{-1, -1}, x])/x + (14*Pi^2*HPL[{-1, 0}, x])/(3*x) + (16*HPL[{-1, 3}, x])/x + (16*Pi^2*HPL[{0, 0}, x])/(3*x) - (16*HPL[{2, 2}, x])/x + (16*HPL[{3, 0}, x])/x + (128*HPL[{3, 1}, x])/x - (16*HPL[{-2, -1, 0}, x])/x - (8*HPL[{-2, 0, 0}, x])/x - (8*HPL[{-1, -2, 0}, x])/x + (32*HPL[{-1, -1, 2}, x])/x - (16*HPL[{-1, 2, 0}, x])/x - (128*HPL[{-1, 2, 1}, x])/x + (4*HPL[{2, 0, 0}, x])/x - (16*HPL[{2, 1, 0}, x])/x - (128*HPL[{2, 1, 1}, x])/x - (8*HPL[{-1, -1, -1, 0}, x])/x - (4*HPL[{-1, -1, 0, 0}, x])/x + (64*HPL[{0, 0, 0, 0}, x])/x + (92*HPL[{-1}, x]*Zeta[3])/x + (20*HPL[{0}, x]*Zeta[3])/x)
  + ep^2*((24*Pi^2*HPL[{-3}, x])/x - (32*Pi^4*HPL[{-1}, x])/(15*x) - (Pi^4*HPL[{0}, x])/x - (128*Pi^2*HPL[{3}, x])/(3*x) + (256*HPL[{5}, x])/x + (64*HPL[{-4, 0}, x])/x - (128*HPL[{-3, 2}, x])/x + (12*Pi^2*HPL[{-2, -1}, x])/x - (28*Pi^2*HPL[{-2, 0}, x])/(3*x) - (32*HPL[{-2, 3}, x])/x + (6*Pi^2*HPL[{-1, -2}, x])/x + (128*Pi^2*HPL[{-1, 2}, x])/(3*x) + (4*Pi^2*HPL[{2, 0}, x])/(3*x) + (128*Pi^2*HPL[{2, 1}, x])/(3*x) + (16*HPL[{2, 3}, x])/x + (64*HPL[{3, 2}, x])/x - (48*HPL[{4, 0}, x])/x - (384*HPL[{4, 1}, x])/x + (32*HPL[{-3, -1, 0}, x])/x + (16*HPL[{-3, 0, 0}, x])/x + (16*HPL[{-2, -2, 0}, x])/x - (64*HPL[{-2, -1, 2}, x])/x + (32*HPL[{-2, 2, 0}, x])/x + (256*HPL[{-2, 2, 1}, x])/x + (8*HPL[{-1, -3, 0}, x])/x - (32*HPL[{-1, -2, 2}, x])/x + (6*Pi^2*HPL[{-1, -1, -1}, x])/x - (14*Pi^2*HPL[{-1, -1, 0}, x])/(3*x) - (16*HPL[{-1, -1, 3}, x])/x - (10*Pi^2*HPL[{-1, 0, 0}, x])/(3*x) - (64*HPL[{-1, 2, 2}, x])/x + (8*HPL[{-1, 3, 0}, x])/x + (64*HPL[{-1, 3, 1}, x])/x - (40*Pi^2*HPL[{0, 0, 0}, x])/(3*x) - (64*HPL[{2, 1, 2}, x])/x - (8*HPL[{2, 2, 0}, x])/x - (64*HPL[{2, 2, 1}, x])/x - (16*HPL[{3, 0, 0}, x])/x + (64*HPL[{3, 1, 0}, x])/x + (512*HPL[{3, 1, 1}, x])/x + (16*HPL[{-2, -1, -1, 0}, x])/x + (8*HPL[{-2, -1, 0, 0}, x])/x + (8*HPL[{-1, -2, -1, 0}, x])/x + (4*HPL[{-1, -2, 0, 0}, x])/x + (8*HPL[{-1, -1, -2, 0}, x])/x - (32*HPL[{-1, -1, -1, 2}, x])/x + (16*HPL[{-1, -1, 2, 0}, x])/x + (128*HPL[{-1, -1, 2, 1}, x])/x + (16*HPL[{-1, 2, 0, 0}, x])/x - (64*HPL[{-1, 2, 1, 0}, x])/x - (512*HPL[{-1, 2, 1, 1}, x])/x - (4*HPL[{2, 0, 0, 0}, x])/x + (16*HPL[{2, 1, 0, 0}, x])/x - (64*HPL[{2, 1, 1, 0}, x])/x - (512*HPL[{2, 1, 1, 1}, x])/x + (8*HPL[{-1, -1, -1, -1, 0}, x])/x + (4*HPL[{-1, -1, -1, 0, 0}, x])/x - (4*HPL[{-1, 0, 0, 0, 0}, x])/x - (160*HPL[{0, 0, 0, 0, 0}, x])/x - (4*Pi^2*Zeta[3])/x - (184*HPL[{-2}, x]*Zeta[3])/x + (112*HPL[{2}, x]*Zeta[3])/x - (92*HPL[{-1, -1}, x]*Zeta[3])/x + (20*HPL[{-1, 0}, x]*Zeta[3])/x - (48*HPL[{0, 0}, x]*Zeta[3])/x - (32*Zeta[5])/x)
  + ALARM*ep^3 );
U16[x -> 1,ep] = (1-x)^(1 - 4*ep) * (
  + 32 + 2/ep^2 + 8/ep - (8*Pi^2)/3 + ep^2*(512 - (128*Pi^2)/3 + (16*Pi^4)/45 - 160*Zeta[3]) + ep*(128 - (32*Pi^2)/3 - 40*Zeta[3])
  + ALARM*ep^3 );
  
U17d[x,ep] = (
  - 3/(ep^3*(-1 + x))
  + ((-3*HPL[{0}, x])/(-1 + x) - (9*HPL[{1}, x])/(-1 + x))/ep^2
  + ((13*Pi^2)/(6*(-1 + x)) - (7*HPL[{2}, x])/(-1 + x) + (2*HPL[{0, 0}, x])/(-1 + x) - (6*HPL[{1, 0}, x])/(-1 + x) - (27*HPL[{1, 1}, x])/(-1 + x))/ep
  + (7*Pi^2*HPL[{0}, x])/(3*(-1 + x)) + (19*Pi^2*HPL[{1}, x])/(3*(-1 + x)) + (6*HPL[{3}, x])/(-1 + x) - (14*HPL[{1, 2}, x])/(-1 + x) + (4*HPL[{2, 0}, x])/(-1 + x) - (17*HPL[{2, 1}, x])/(-1 + x) + (4*HPL[{1, 0, 0}, x])/(-1 + x) - (12*HPL[{1, 1, 0}, x])/(-1 + x) - (81*HPL[{1, 1, 1}, x])/(-1 + x) + (23*Zeta[3])/(-1 + x)
  + ep*(-Pi^4/(180*(-1 + x)) + (6*Pi^2*HPL[{2}, x])/(-1 + x) - (4*HPL[{4}, x])/(-1 + x) - (2*Pi^2*HPL[{0, 0}, x])/(3*(-1 + x)) + (14*Pi^2*HPL[{1, 0}, x])/(3*(-1 + x)) + (56*Pi^2*HPL[{1, 1}, x])/(3*(-1 + x)) + (12*HPL[{1, 3}, x])/(-1 + x) + (8*HPL[{2, 2}, x])/(-1 + x) + (18*HPL[{3, 1}, x])/(-1 + x) - (28*HPL[{1, 1, 2}, x])/(-1 + x) + (8*HPL[{1, 2, 0}, x])/(-1 + x) - (34*HPL[{1, 2, 1}, x])/(-1 + x) - (4*HPL[{2, 0, 0}, x])/(-1 + x) + (8*HPL[{2, 1, 0}, x])/(-1 + x) - (43*HPL[{2, 1, 1}, x])/(-1 + x) - (4*HPL[{0, 0, 0, 0}, x])/(-1 + x) + (8*HPL[{1, 1, 0, 0}, x])/(-1 + x) - (24*HPL[{1, 1, 1, 0}, x])/(-1 + x) - (243*HPL[{1, 1, 1, 1}, x])/(-1 + x) + (14*HPL[{0}, x]*Zeta[3])/(-1 + x) + (70*HPL[{1}, x]*Zeta[3])/(-1 + x))
  + ep^2*((-37*Pi^4*HPL[{0}, x])/(30*(-1 + x)) - (Pi^4*HPL[{1}, x])/(90*(-1 + x)) - (4*Pi^2*HPL[{3}, x])/(-1 + x) + (12*Pi^2*HPL[{1, 2}, x])/(-1 + x) - (8*HPL[{1, 4}, x])/(-1 + x) - (4*Pi^2*HPL[{2, 0}, x])/(-1 + x) + (14*Pi^2*HPL[{2, 1}, x])/(-1 + x) - (8*HPL[{2, 3}, x])/(-1 + x) - (8*HPL[{4, 0}, x])/(-1 + x) - (20*HPL[{4, 1}, x])/(-1 + x) - (8*Pi^2*HPL[{0, 0, 0}, x])/(3*(-1 + x)) - (4*Pi^2*HPL[{1, 0, 0}, x])/(3*(-1 + x)) + (28*Pi^2*HPL[{1, 1, 0}, x])/(3*(-1 + x)) + (166*Pi^2*HPL[{1, 1, 1}, x])/(3*(-1 + x)) + (24*HPL[{1, 1, 3}, x])/(-1 + x) + (16*HPL[{1, 2, 2}, x])/(-1 + x) + (36*HPL[{1, 3, 1}, x])/(-1 + x) + (16*HPL[{2, 1, 2}, x])/(-1 + x) - (8*HPL[{2, 2, 0}, x])/(-1 + x) + (16*HPL[{2, 2, 1}, x])/(-1 + x) + (54*HPL[{3, 1, 1}, x])/(-1 + x) - (56*HPL[{1, 1, 1, 2}, x])/(-1 + x) + (16*HPL[{1, 1, 2, 0}, x])/(-1 + x) - (68*HPL[{1, 1, 2, 1}, x])/(-1 + x) - (8*HPL[{1, 2, 0, 0}, x])/(-1 + x) + (16*HPL[{1, 2, 1, 0}, x])/(-1 + x) - (86*HPL[{1, 2, 1, 1}, x])/(-1 + x) + (4*HPL[{2, 0, 0, 0}, x])/(-1 + x) - (8*HPL[{2, 1, 0, 0}, x])/(-1 + x) + (16*HPL[{2, 1, 1, 0}, x])/(-1 + x) - (113*HPL[{2, 1, 1, 1}, x])/(-1 + x) + (12*HPL[{0, 0, 0, 0, 0}, x])/(-1 + x) - (8*HPL[{1, 0, 0, 0, 0}, x])/(-1 + x) + (16*HPL[{1, 1, 1, 0, 0}, x])/(-1 + x) - (48*HPL[{1, 1, 1, 1, 0}, x])/(-1 + x) - (729*HPL[{1, 1, 1, 1, 1}, x])/(-1 + x) - (55*Pi^2*Zeta[3])/(3*(-1 + x)) + (40*HPL[{2}, x]*Zeta[3])/(-1 + x) - (12*HPL[{0, 0}, x]*Zeta[3])/(-1 + x) + (28*HPL[{1, 0}, x]*Zeta[3])/(-1 + x) + (212*HPL[{1, 1}, x]*Zeta[3])/(-1 + x) + (95*Zeta[5])/(-1 + x))
  + ALARM*ep^3 );
U17d[x -> 1,ep] = (1-x)^(-1 - 3*ep) * (
  + 3/ep^3 - (2*Pi^2)/ep - 24*Zeta[3] + ep^2*(16*Pi^2*Zeta[3] - 144*Zeta[5])
  + ALARM*ep^3 );

U17t[x,ep] = (
  + 1/(ep^3*(-1 + x))
  + ((2*HPL[{0}, x])/(-1 + x) + (4*HPL[{1}, x])/(-1 + x))/ep^2
  + ((-2*Pi^2)/(-1 + x) + (8*HPL[{2}, x])/(-1 + x) + (4*HPL[{1, 0}, x])/(-1 + x) + (16*HPL[{1, 1}, x])/(-1 + x))/ep
  - (7*Pi^2*HPL[{0}, x])/(3*(-1 + x)) - (19*Pi^2*HPL[{1}, x])/(3*(-1 + x)) + (16*HPL[{1, 2}, x])/(-1 + x) + (2*HPL[{2, 0}, x])/(-1 + x) + (32*HPL[{2, 1}, x])/(-1 + x) - (2*HPL[{0, 0, 0}, x])/(-1 + x) - (2*HPL[{1, 0, 0}, x])/(-1 + x) + (10*HPL[{1, 1, 0}, x])/(-1 + x) + (64*HPL[{1, 1, 1}, x])/(-1 + x) - (24*Zeta[3])/(-1 + x)
  + ep*((23*Pi^4)/(18*(-1 + x)) - (10*Pi^2*HPL[{2}, x])/(-1 + x) - (8*HPL[{4}, x])/(-1 + x) + (5*Pi^2*HPL[{0, 0}, x])/(3*(-1 + x)) - (3*Pi^2*HPL[{1, 0}, x])/(-1 + x) - (68*Pi^2*HPL[{1, 1}, x])/(3*(-1 + x)) - (8*HPL[{1, 3}, x])/(-1 + x) + (8*HPL[{2, 2}, x])/(-1 + x) + (2*HPL[{3, 0}, x])/(-1 + x) + (40*HPL[{1, 1, 2}, x])/(-1 + x) + (6*HPL[{1, 2, 0}, x])/(-1 + x) + (64*HPL[{1, 2, 1}, x])/(-1 + x) - (6*HPL[{2, 0, 0}, x])/(-1 + x) + (12*HPL[{2, 1, 0}, x])/(-1 + x) + (128*HPL[{2, 1, 1}, x])/(-1 + x) + (4*HPL[{0, 0, 0, 0}, x])/(-1 + x) - (10*HPL[{1, 1, 0, 0}, x])/(-1 + x) + (32*HPL[{1, 1, 1, 0}, x])/(-1 + x) + (256*HPL[{1, 1, 1, 1}, x])/(-1 + x) - (18*HPL[{0}, x]*Zeta[3])/(-1 + x) - (66*HPL[{1}, x]*Zeta[3])/(-1 + x))
  + ep^2*((Pi^4*HPL[{0}, x])/(5*(-1 + x)) + (124*Pi^4*HPL[{1}, x])/(45*(-1 + x)) - (2*Pi^2*HPL[{3}, x])/(3*(-1 + x)) + (16*HPL[{5}, x])/(-1 + x) - (62*Pi^2*HPL[{1, 2}, x])/(3*(-1 + x)) + (8*Pi^2*HPL[{2, 0}, x])/(3*(-1 + x)) - (124*Pi^2*HPL[{2, 1}, x])/(3*(-1 + x)) - (24*HPL[{2, 3}, x])/(-1 + x) + (8*HPL[{3, 2}, x])/(-1 + x) - (6*HPL[{4, 0}, x])/(-1 + x) - (32*HPL[{4, 1}, x])/(-1 + x) - (Pi^2*HPL[{0, 0, 0}, x])/(-1 + x) + (7*Pi^2*HPL[{1, 0, 0}, x])/(3*(-1 + x)) - (10*Pi^2*HPL[{1, 1, 0}, x])/(3*(-1 + x)) - (260*Pi^2*HPL[{1, 1, 1}, x])/(3*(-1 + x)) - (40*HPL[{1, 1, 3}, x])/(-1 + x) + (24*HPL[{1, 2, 2}, x])/(-1 + x) - (2*HPL[{1, 3, 0}, x])/(-1 + x) - (32*HPL[{1, 3, 1}, x])/(-1 + x) + (48*HPL[{2, 1, 2}, x])/(-1 + x) + (8*HPL[{2, 2, 0}, x])/(-1 + x) + (32*HPL[{2, 2, 1}, x])/(-1 + x) + (2*HPL[{3, 0, 0}, x])/(-1 + x) + (4*HPL[{3, 1, 0}, x])/(-1 + x) + (128*HPL[{1, 1, 1, 2}, x])/(-1 + x) + (20*HPL[{1, 1, 2, 0}, x])/(-1 + x) + (160*HPL[{1, 1, 2, 1}, x])/(-1 + x) - (10*HPL[{1, 2, 0, 0}, x])/(-1 + x) + (28*HPL[{1, 2, 1, 0}, x])/(-1 + x) + (256*HPL[{1, 2, 1, 1}, x])/(-1 + x) + (10*HPL[{2, 0, 0, 0}, x])/(-1 + x) - (20*HPL[{2, 1, 0, 0}, x])/(-1 + x) + (56*HPL[{2, 1, 1, 0}, x])/(-1 + x) + (512*HPL[{2, 1, 1, 1}, x])/(-1 + x) - (6*HPL[{0, 0, 0, 0, 0}, x])/(-1 + x) + (2*HPL[{1, 0, 0, 0, 0}, x])/(-1 + x) + (10*HPL[{1, 1, 0, 0, 0}, x])/(-1 + x) - (40*HPL[{1, 1, 1, 0, 0}, x])/(-1 + x) + (120*HPL[{1, 1, 1, 1, 0}, x])/(-1 + x) + (1024*HPL[{1, 1, 1, 1, 1}, x])/(-1 + x) + (106*Pi^2*Zeta[3])/(3*(-1 + x)) - (92*HPL[{2}, x]*Zeta[3])/(-1 + x) + (14*HPL[{0, 0}, x]*Zeta[3])/(-1 + x) - (22*HPL[{1, 0}, x]*Zeta[3])/(-1 + x) - (224*HPL[{1, 1}, x]*Zeta[3])/(-1 + x) - (40*Zeta[5])/(-1 + x))
  + ALARM*ep^3 );
U17t[x -> 1,ep] = (1-x)^(-1 - 4*ep) * (
  - ep^(-3) + (4*Pi^2)/(3*ep) - (8*ep*Pi^4)/45 + 20*Zeta[3] + ep^2*((-80*Pi^2*Zeta[3])/3 + 204*Zeta[5])
  + ALARM*ep^3 );

U18[x,ep] = (
  + 1/(2*ep^3*x)
  + (HPL[{0}, x]/(1 + x) + (2*HPL[{1}, x])/x)/ep^2
  + (-Pi^2/(2*x) - Pi^2/(1 + x) + (4*HPL[{2}, x])/(1 + x) + (2*HPL[{-1, 0}, x])/x - (4*HPL[{-1, 0}, x])/(1 + x) - (2*HPL[{0, 0}, x])/x + (2*HPL[{0, 0}, x])/(1 + x) + (8*HPL[{1, 1}, x])/x)/ep
  - (2*Pi^2*HPL[{-1}, x])/x + (4*Pi^2*HPL[{-1}, x])/(1 + x) + (4*Pi^2*HPL[{0}, x])/(3*x) - (8*Pi^2*HPL[{0}, x])/(3*(1 + x)) - (7*Pi^2*HPL[{1}, x])/(3*x) - (8*HPL[{3}, x])/x + (8*HPL[{3}, x])/(1 + x) - (4*HPL[{-2, 0}, x])/x + (8*HPL[{-1, 2}, x])/x - (16*HPL[{-1, 2}, x])/(1 + x) + (2*HPL[{2, 0}, x])/x + (16*HPL[{2, 1}, x])/(1 + x) - (8*HPL[{-1, -1, 0}, x])/x + (16*HPL[{-1, -1, 0}, x])/(1 + x) + (4*HPL[{-1, 0, 0}, x])/x - (8*HPL[{-1, 0, 0}, x])/(1 + x) + (8*HPL[{0, 0, 0}, x])/x - (6*HPL[{0, 0, 0}, x])/(1 + x) - (2*HPL[{1, 0, 0}, x])/x + (2*HPL[{1, 1, 0}, x])/x + (32*HPL[{1, 1, 1}, x])/x - Zeta[3]/x - (30*Zeta[3])/(1 + x)
  + ep*((-7*Pi^4)/(20*x) + Pi^4/(1 + x) + (4*Pi^2*HPL[{-2}, x])/x - (2*Pi^2*HPL[{2}, x])/(3*x) - (14*Pi^2*HPL[{2}, x])/(3*(1 + x)) + (32*HPL[{4}, x])/x - (24*HPL[{4}, x])/(1 + x) + (8*HPL[{-3, 0}, x])/x - (16*HPL[{-2, 2}, x])/x + (8*Pi^2*HPL[{-1, -1}, x])/x - (16*Pi^2*HPL[{-1, -1}, x])/(1 + x) - (16*Pi^2*HPL[{-1, 0}, x])/(3*x) + (32*Pi^2*HPL[{-1, 0}, x])/(3*(1 + x)) + (16*HPL[{-1, 3}, x])/x - (32*HPL[{-1, 3}, x])/(1 + x) - (10*Pi^2*HPL[{0, 0}, x])/(3*x) + (7*Pi^2*HPL[{0, 0}, x])/(3*(1 + x)) + (5*Pi^2*HPL[{1, 0}, x])/(3*x) - (10*Pi^2*HPL[{1, 1}, x])/x - (8*HPL[{1, 3}, x])/x + (8*HPL[{2, 2}, x])/x - (8*HPL[{3, 0}, x])/x + (6*HPL[{3, 0}, x])/(1 + x) - (32*HPL[{3, 1}, x])/x + (32*HPL[{3, 1}, x])/(1 + x) + (16*HPL[{-2, -1, 0}, x])/x - (8*HPL[{-2, 0, 0}, x])/x - (32*HPL[{-1, -1, 2}, x])/x + (64*HPL[{-1, -1, 2}, x])/(1 + x) + (32*HPL[{-1, 2, 1}, x])/x - (64*HPL[{-1, 2, 1}, x])/(1 + x) + (8*HPL[{1, 1, 2}, x])/x + (2*HPL[{1, 2, 0}, x])/x + (2*HPL[{2, 0, 0}, x])/x - (4*HPL[{2, 0, 0}, x])/(1 + x) + (4*HPL[{2, 1, 0}, x])/x + (4*HPL[{2, 1, 0}, x])/(1 + x) + (64*HPL[{2, 1, 1}, x])/(1 + x) + (32*HPL[{-1, -1, -1, 0}, x])/x - (64*HPL[{-1, -1, -1, 0}, x])/(1 + x) - (16*HPL[{-1, -1, 0, 0}, x])/x + (32*HPL[{-1, -1, 0, 0}, x])/(1 + x) - (12*HPL[{-1, 0, 0, 0}, x])/x + (24*HPL[{-1, 0, 0, 0}, x])/(1 + x) - (24*HPL[{0, 0, 0, 0}, x])/x + (12*HPL[{0, 0, 0, 0}, x])/(1 + x) + (4*HPL[{1, 0, 0, 0}, x])/x - (6*HPL[{1, 1, 0, 0}, x])/x + (12*HPL[{1, 1, 1, 0}, x])/x + (128*HPL[{1, 1, 1, 1}, x])/x - (60*HPL[{-1}, x]*Zeta[3])/x + (120*HPL[{-1}, x]*Zeta[3])/(1 + x) - (16*HPL[{0}, x]*Zeta[3])/(1 + x) - (18*HPL[{1}, x]*Zeta[3])/x)
  + ep^2*((-8*Pi^2*HPL[{-3}, x])/x + (2*Pi^4*HPL[{-1}, x])/x - (4*Pi^4*HPL[{-1}, x])/(1 + x) + (4*Pi^4*HPL[{0}, x])/(5*x) + (Pi^4*HPL[{0}, x])/(10*(1 + x)) + (Pi^4*HPL[{1}, x])/(5*x) + (12*Pi^2*HPL[{3}, x])/x - (34*Pi^2*HPL[{3}, x])/(3*(1 + x)) - (96*HPL[{5}, x])/x + (48*HPL[{5}, x])/(1 + x) - (16*HPL[{-4, 0}, x])/x + (32*HPL[{-3, 2}, x])/x - (16*Pi^2*HPL[{-2, -1}, x])/x + (32*Pi^2*HPL[{-2, 0}, x])/(3*x) - (32*HPL[{-2, 3}, x])/x - (28*Pi^2*HPL[{-1, 2}, x])/(3*x) + (56*Pi^2*HPL[{-1, 2}, x])/(3*(1 + x)) - (48*HPL[{-1, 4}, x])/x + (96*HPL[{-1, 4}, x])/(1 + x) - (2*Pi^2*HPL[{1, 2}, x])/(3*x) + (16*HPL[{1, 4}, x])/x - (4*Pi^2*HPL[{2, 0}, x])/x + (10*Pi^2*HPL[{2, 0}, x])/(3*(1 + x)) - (4*Pi^2*HPL[{2, 1}, x])/(3*x) - (20*Pi^2*HPL[{2, 1}, x])/(1 + x) + (8*HPL[{2, 3}, x])/x - (16*HPL[{2, 3}, x])/(1 + x) - (32*HPL[{3, 2}, x])/x + (24*HPL[{3, 2}, x])/(1 + x) + (24*HPL[{4, 0}, x])/x - (14*HPL[{4, 0}, x])/(1 + x) + (128*HPL[{4, 1}, x])/x - (96*HPL[{4, 1}, x])/(1 + x) - (32*HPL[{-3, -1, 0}, x])/x + (16*HPL[{-3, 0, 0}, x])/x + (64*HPL[{-2, -1, 2}, x])/x - (64*HPL[{-2, 2, 1}, x])/x - (32*Pi^2*HPL[{-1, -1, -1}, x])/x + (64*Pi^2*HPL[{-1, -1, -1}, x])/(1 + x) + (64*Pi^2*HPL[{-1, -1, 0}, x])/(3*x) - (128*Pi^2*HPL[{-1, -1, 0}, x])/(3*(1 + x)) - (64*HPL[{-1, -1, 3}, x])/x + (128*HPL[{-1, -1, 3}, x])/(1 + x) + (14*Pi^2*HPL[{-1, 0, 0}, x])/(3*x) - (28*Pi^2*HPL[{-1, 0, 0}, x])/(3*(1 + x)) + (12*HPL[{-1, 3, 0}, x])/x - (24*HPL[{-1, 3, 0}, x])/(1 + x) + (64*HPL[{-1, 3, 1}, x])/x - (128*HPL[{-1, 3, 1}, x])/(1 + x) + (8*Pi^2*HPL[{0, 0, 0}, x])/x - (7*Pi^2*HPL[{0, 0, 0}, x])/(3*(1 + x)) - (Pi^2*HPL[{1, 0, 0}, x])/x + (8*Pi^2*HPL[{1, 1, 0}, x])/(3*x) - (124*Pi^2*HPL[{1, 1, 1}, x])/(3*x) - (24*HPL[{1, 1, 3}, x])/x + (8*HPL[{1, 2, 2}, x])/x - (6*HPL[{1, 3, 0}, x])/x - (32*HPL[{1, 3, 1}, x])/x + (16*HPL[{2, 1, 2}, x])/x + (16*HPL[{2, 1, 2}, x])/(1 + x) + (4*HPL[{2, 2, 0}, x])/(1 + x) + (32*HPL[{2, 2, 1}, x])/x - (2*HPL[{3, 0, 0}, x])/(1 + x) - (24*HPL[{3, 1, 0}, x])/x + (20*HPL[{3, 1, 0}, x])/(1 + x) - (128*HPL[{3, 1, 1}, x])/x + (128*HPL[{3, 1, 1}, x])/(1 + x) - (64*HPL[{-2, -1, -1, 0}, x])/x + (32*HPL[{-2, -1, 0, 0}, x])/x + (24*HPL[{-2, 0, 0, 0}, x])/x + (128*HPL[{-1, -1, -1, 2}, x])/x - (256*HPL[{-1, -1, -1, 2}, x])/(1 + x) - (128*HPL[{-1, -1, 2, 1}, x])/x + (256*HPL[{-1, -1, 2, 1}, x])/(1 + x) - (8*HPL[{-1, 2, 0, 0}, x])/x + (16*HPL[{-1, 2, 0, 0}, x])/(1 + x) + (8*HPL[{-1, 2, 1, 0}, x])/x - (16*HPL[{-1, 2, 1, 0}, x])/(1 + x) + (128*HPL[{-1, 2, 1, 1}, x])/x - (256*HPL[{-1, 2, 1, 1}, x])/(1 + x) + (48*HPL[{1, 1, 1, 2}, x])/x + (8*HPL[{1, 1, 2, 0}, x])/x + (32*HPL[{1, 1, 2, 1}, x])/x + (2*HPL[{1, 2, 0, 0}, x])/x + (4*HPL[{1, 2, 1, 0}, x])/x - (6*HPL[{2, 0, 0, 0}, x])/x + (8*HPL[{2, 0, 0, 0}, x])/(1 + x) + (4*HPL[{2, 1, 0, 0}, x])/x - (12*HPL[{2, 1, 0, 0}, x])/(1 + x) + (8*HPL[{2, 1, 1, 0}, x])/x + (24*HPL[{2, 1, 1, 0}, x])/(1 + x) + (256*HPL[{2, 1, 1, 1}, x])/(1 + x) - (128*HPL[{-1, -1, -1, -1, 0}, x])/x + (256*HPL[{-1, -1, -1, -1, 0}, x])/(1 + x) + (64*HPL[{-1, -1, -1, 0, 0}, x])/x - (128*HPL[{-1, -1, -1, 0, 0}, x])/(1 + x) + (48*HPL[{-1, -1, 0, 0, 0}, x])/x - (96*HPL[{-1, -1, 0, 0, 0}, x])/(1 + x) + (24*HPL[{-1, 0, 0, 0, 0}, x])/x - (48*HPL[{-1, 0, 0, 0, 0}, x])/(1 + x) + (64*HPL[{0, 0, 0, 0, 0}, x])/x - (22*HPL[{0, 0, 0, 0, 0}, x])/(1 + x) - (6*HPL[{1, 0, 0, 0, 0}, x])/x + (10*HPL[{1, 1, 0, 0, 0}, x])/x - (20*HPL[{1, 1, 1, 0, 0}, x])/x + (56*HPL[{1, 1, 1, 1, 0}, x])/x + (512*HPL[{1, 1, 1, 1, 1}, x])/x + (Pi^2*Zeta[3])/x + (52*Pi^2*Zeta[3])/(1 + x) + (120*HPL[{-2}, x]*Zeta[3])/x - (20*HPL[{2}, x]*Zeta[3])/x - (36*HPL[{2}, x]*Zeta[3])/(1 + x) + (240*HPL[{-1, -1}, x]*Zeta[3])/x - (480*HPL[{-1, -1}, x]*Zeta[3])/(1 + x) - (32*HPL[{-1, 0}, x]*Zeta[3])/x + (64*HPL[{-1, 0}, x]*Zeta[3])/(1 + x) + (4*HPL[{0, 0}, x]*Zeta[3])/x + (10*HPL[{0, 0}, x]*Zeta[3])/(1 + x) + (14*HPL[{1, 0}, x]*Zeta[3])/x - (92*HPL[{1, 1}, x]*Zeta[3])/x - (3*Zeta[5])/x - (154*Zeta[5])/(1 + x))
  + ALARM*ep^3 );
U18[x -> 1,ep] = (1-x)^(-4*ep) * (
  + 1/(2*ep^3) - (2*Pi^2)/(3*ep) + (4*ep*Pi^4)/45 - 10*Zeta[3] + ep^2*((40*Pi^2*Zeta[3])/3 - 102*Zeta[5])
  + ALARM*ep^3 );

U19[x,ep] = (
  - (4*HPL[{0}, x])/(ep^2*(1 + x))
  + ((10*Pi^2)/(3*(1 + x)) - (16*HPL[{2}, x])/(1 + x) + (8*HPL[{-1, 0}, x])/(1 + x) - (4*HPL[{0, 0}, x])/(1 + x))/ep
  - (20*Pi^2*HPL[{-1}, x])/(3*(1 + x)) + (6*Pi^2*HPL[{0}, x])/(1 + x) - (16*HPL[{3}, x])/(1 + x) - (8*HPL[{-2, 0}, x])/(1 + x) + (32*HPL[{-1, 2}, x])/(1 + x) - (8*HPL[{2, 0}, x])/(1 + x) - (64*HPL[{2, 1}, x])/(1 + x) - (16*HPL[{-1, -1, 0}, x])/(1 + x) + (8*HPL[{-1, 0, 0}, x])/(1 + x) + (16*HPL[{0, 0, 0}, x])/(1 + x) + (68*Zeta[3])/(1 + x)
  + ep*((-38*Pi^4)/(15*(1 + x)) + (6*Pi^2*HPL[{-2}, x])/(1 + x) + (64*Pi^2*HPL[{2}, x])/(3*(1 + x)) + (64*HPL[{4}, x])/(1 + x) + (8*HPL[{-3, 0}, x])/(1 + x) - (32*HPL[{-2, 2}, x])/(1 + x) + (40*Pi^2*HPL[{-1, -1}, x])/(3*(1 + x)) - (12*Pi^2*HPL[{-1, 0}, x])/(1 + x) + (32*HPL[{-1, 3}, x])/(1 + x) - (6*Pi^2*HPL[{0, 0}, x])/(1 + x) - (32*HPL[{2, 2}, x])/(1 + x) - (8*HPL[{3, 0}, x])/(1 + x) - (64*HPL[{3, 1}, x])/(1 + x) + (8*HPL[{-2, -1, 0}, x])/(1 + x) + (4*HPL[{-2, 0, 0}, x])/(1 + x) + (16*HPL[{-1, -2, 0}, x])/(1 + x) - (64*HPL[{-1, -1, 2}, x])/(1 + x) + (16*HPL[{-1, 2, 0}, x])/(1 + x) + (128*HPL[{-1, 2, 1}, x])/(1 + x) + (8*HPL[{2, 0, 0}, x])/(1 + x) - (32*HPL[{2, 1, 0}, x])/(1 + x) - (256*HPL[{2, 1, 1}, x])/(1 + x) + (32*HPL[{-1, -1, -1, 0}, x])/(1 + x) - (16*HPL[{-1, -1, 0, 0}, x])/(1 + x) - (32*HPL[{-1, 0, 0, 0}, x])/(1 + x) - (36*HPL[{0, 0, 0, 0}, x])/(1 + x) - (136*HPL[{-1}, x]*Zeta[3])/(1 + x) + (12*HPL[{0}, x]*Zeta[3])/(1 + x))
  + ep^2*((-6*Pi^2*HPL[{-3}, x])/(1 + x) + (76*Pi^4*HPL[{-1}, x])/(15*(1 + x)) - (2*Pi^4*HPL[{0}, x])/(9*(1 + x)) + (64*Pi^2*HPL[{3}, x])/(3*(1 + x)) - (144*HPL[{5}, x])/(1 + x) - (8*HPL[{-4, 0}, x])/(1 + x) + (32*HPL[{-3, 2}, x])/(1 + x) - (6*Pi^2*HPL[{-2, -1}, x])/(1 + x) + (14*Pi^2*HPL[{-2, 0}, x])/(3*(1 + x)) + (16*HPL[{-2, 3}, x])/(1 + x) - (12*Pi^2*HPL[{-1, -2}, x])/(1 + x) - (128*Pi^2*HPL[{-1, 2}, x])/(3*(1 + x)) - (128*HPL[{-1, 4}, x])/(1 + x) + (8*Pi^2*HPL[{2, 0}, x])/(3*(1 + x)) + (256*Pi^2*HPL[{2, 1}, x])/(3*(1 + x)) + (32*HPL[{2, 3}, x])/(1 + x) - (32*HPL[{3, 2}, x])/(1 + x) + (32*HPL[{4, 0}, x])/(1 + x) + (256*HPL[{4, 1}, x])/(1 + x) - (8*HPL[{-3, -1, 0}, x])/(1 + x) - (4*HPL[{-3, 0, 0}, x])/(1 + x) - (8*HPL[{-2, -2, 0}, x])/(1 + x) + (32*HPL[{-2, -1, 2}, x])/(1 + x) - (16*HPL[{-2, 2, 0}, x])/(1 + x) - (128*HPL[{-2, 2, 1}, x])/(1 + x) - (16*HPL[{-1, -3, 0}, x])/(1 + x) + (64*HPL[{-1, -2, 2}, x])/(1 + x) - (80*Pi^2*HPL[{-1, -1, -1}, x])/(3*(1 + x)) + (24*Pi^2*HPL[{-1, -1, 0}, x])/(1 + x) - (64*HPL[{-1, -1, 3}, x])/(1 + x) + (12*Pi^2*HPL[{-1, 0, 0}, x])/(1 + x) + (64*HPL[{-1, 2, 2}, x])/(1 + x) + (16*HPL[{-1, 3, 0}, x])/(1 + x) + (128*HPL[{-1, 3, 1}, x])/(1 + x) + (22*Pi^2*HPL[{0, 0, 0}, x])/(3*(1 + x)) - (128*HPL[{2, 1, 2}, x])/(1 + x) - (16*HPL[{2, 2, 0}, x])/(1 + x) - (128*HPL[{2, 2, 1}, x])/(1 + x) + (8*HPL[{3, 0, 0}, x])/(1 + x) - (32*HPL[{3, 1, 0}, x])/(1 + x) - (256*HPL[{3, 1, 1}, x])/(1 + x) - (8*HPL[{-2, -1, -1, 0}, x])/(1 + x) - (4*HPL[{-2, -1, 0, 0}, x])/(1 + x) - (16*HPL[{-1, -2, -1, 0}, x])/(1 + x) - (8*HPL[{-1, -2, 0, 0}, x])/(1 + x) - (32*HPL[{-1, -1, -2, 0}, x])/(1 + x) + (128*HPL[{-1, -1, -1, 2}, x])/(1 + x) - (32*HPL[{-1, -1, 2, 0}, x])/(1 + x) - (256*HPL[{-1, -1, 2, 1}, x])/(1 + x) - (16*HPL[{-1, 2, 0, 0}, x])/(1 + x) + (64*HPL[{-1, 2, 1, 0}, x])/(1 + x) + (512*HPL[{-1, 2, 1, 1}, x])/(1 + x) - (8*HPL[{2, 0, 0, 0}, x])/(1 + x) + (32*HPL[{2, 1, 0, 0}, x])/(1 + x) - (128*HPL[{2, 1, 1, 0}, x])/(1 + x) - (1024*HPL[{2, 1, 1, 1}, x])/(1 + x) - (64*HPL[{-1, -1, -1, -1, 0}, x])/(1 + x) + (32*HPL[{-1, -1, -1, 0, 0}, x])/(1 + x) + (64*HPL[{-1, -1, 0, 0, 0}, x])/(1 + x) + (72*HPL[{-1, 0, 0, 0, 0}, x])/(1 + x) + (72*HPL[{0, 0, 0, 0, 0}, x])/(1 + x) - (232*Pi^2*Zeta[3])/(3*(1 + x)) + (92*HPL[{-2}, x]*Zeta[3])/(1 + x) + (224*HPL[{2}, x]*Zeta[3])/(1 + x) + (272*HPL[{-1, -1}, x]*Zeta[3])/(1 + x) - (24*HPL[{-1, 0}, x]*Zeta[3])/(1 + x) + (4*HPL[{0, 0}, x]*Zeta[3])/(1 + x) + (188*Zeta[5])/(1 + x))
  + ALARM*ep^3 );
U19[x -> 1,ep] = (1-x)^(1-4*ep) * (
  + 32 + 2/ep^2 + 8/ep - (8*Pi^2)/3 + ep^2*(512 - (128*Pi^2)/3 + (16*Pi^4)/45 - 160*Zeta[3]) + ep*(128 - (32*Pi^2)/3 - 40*Zeta[3])
  + ALARM*ep^3 );

