-16*ep**2/(3*d**2)