function [eo]=elo(r,v,muc) %Función para obtener los elementos orbitales a partir de la posición y velocidad
AU=149.598E6;
eo=zeros(1,6);%Inicializamos el vector de los elementos orbitales
r;
v;
h=cross(r,v);

%EC
EC=cross(v,h)/muc-r/norm(r);
eo(1)=norm(EC);
%SMA(UA)
epsilon=norm(v)^2/2-muc/norm(r);
SMA=-muc/2/epsilon;
eo(2)=SMA/AU;
%TA(grados)
if r'*v>0
    TA=acos(r'*EC/norm(r)/norm(EC))*180/pi;
else
    TA=(2*pi-acos(r'*EC/norm(r)/norm(EC)))*180/pi;
end
eo(6)=TA;
%INC(grados)
INC=acos(h(3)/norm(h))*180/pi;
eo(3)=INC;
%%OM(grados)
K=[0 0 1]';
n=cross(K,h)/norm(cross(K,h));
if n(2)>0
OM=acos(n(1))*180/pi;
else 
    OM=(2*pi-acos(n(1)))*180/pi;
end
eo(4)=OM;
%%W(grados)
if EC(3)>0
    W=acos(EC'*n/norm(EC))*180/pi;
else
    W=(2*pi-acos(EC'*n/norm(EC)))*180/pi;
end
eo(5)=W;