function res = dHz(phi,R,x,z)
  res = R*(R-x.*cos(phi))/sqrt(R^2-2*x.*R.*cos(phi)+x.^2 +z.^2).^3;
endfunction
