function feq=HillRad(yi,C) global GEO wt mu U=pot(0,yi,0,GEO,mu); feq=1/2*wt^2*yi^2+U-C-1/2*(sqrt(mu/yi)-wt*yi)^2; end