1-4*r**3*Qa*D1_d/d**5 + r**3*U_md*D1_d/d**3 - r**3*D1_U_md/(2*d**2) + r**3*D1_Qa/d**4 + 4*r*Qa*D1_d/d**3 + r*D1_U_md/2 - 2*r*D1_Qa/d**2 2
