d/dh[ h^2 / (r*(1+sqrt(1-e*h^2/r^2)))]