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