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