-1/12 + x/3 + (-1/6 + x/3 - (2*x^2)/3)*HPL[{0}, x] + 
 (-1/3 + (2*x)/3 - (2*x^2)/3)*HPL[{1}, x]
