%Obtencion de punto de equilibrio en el eje X(Saddle equilibrium points) function feq=eqY(yi) global GEO wt mu [GU U]=Gpot(0,yi,0,GEO,mu); feq=wt^2*yi+GU(2); end