d**4*np.log(d)/4 - d**4/16 - d**2*(d**2/8 - 1/8) + d**2/8 - 1/16