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