solve (v/2-gm/rad=rad^2*v_phi^2/(2*(rad+h)^2)-gm/(rad+h)) for h