gradient of 4*K(4*rho/((1+rho)^2+z^2))/sqrt((rho+1)^2+z^2)