-6 - 33/(2*x) - (21*x)/2 + (33*x^2)/2 + nl*(2 - x^(-1) - x + x^2) + 
 (9/(-1 + x) - 9/x - 27*x + 9*x^2)*HPL[{0}, x] + 
 (36 - 18/x - 18*x + 18*x^2)*HPL[{1}, x] + 
 (121/16 - (11*nl)/12 + nl^2/36 - (3*Pi^2)/2)*PlusDistribution[-1, 1 - x] + 
 (33/2 - nl)*PlusDistribution[0, 1 - x] + 18*PlusDistribution[1, 1 - x] + 
 ep*(-6 - 33/(2*x) - (21*x)/2 + (33*x^2)/2 + nl*(2 - x^(-1) - x + x^2) + 
   (9/(-1 + x) - 9/x - 27*x + 9*x^2)*HPL[{0}, x] + 
   (36 - 18/x - 18*x + 18*x^2)*HPL[{1}, x] + 
   (121/16 - (11*nl)/12 + nl^2/36 - (3*Pi^2)/2)*PlusDistribution[-1, 1 - x] + 
   (33/2 - nl)*PlusDistribution[0, 1 - x] + 18*PlusDistribution[1, 1 - x]) + 
 ep^2*(-6 - 33/(2*x) - (21*x)/2 + (33*x^2)/2 + nl*(2 - x^(-1) - x + x^2) + 
   (9/(-1 + x) - 9/x - 27*x + 9*x^2)*HPL[{0}, x] + 
   (36 - 18/x - 18*x + 18*x^2)*HPL[{1}, x] + 
   (121/16 - (11*nl)/12 + nl^2/36 - (3*Pi^2)/2)*PlusDistribution[-1, 1 - x] + 
   (33/2 - nl)*PlusDistribution[0, 1 - x] + 18*PlusDistribution[1, 1 - x]) + 
 ep^3*(-6 - 33/(2*x) - (21*x)/2 + (33*x^2)/2 + nl*(2 - x^(-1) - x + x^2) + 
   (9/(-1 + x) - 9/x - 27*x + 9*x^2)*HPL[{0}, x] + 
   (36 - 18/x - 18*x + 18*x^2)*HPL[{1}, x] + 
   (121/16 - (11*nl)/12 + nl^2/36 - (3*Pi^2)/2)*PlusDistribution[-1, 1 - x] + 
   (33/2 - nl)*PlusDistribution[0, 1 - x] + 18*PlusDistribution[1, 1 - x])
