function [coord]=rv(eo)
global mus
EC=eo(1);
A=eo(2);
INC=eo(3);
OM=eo(4);
W=eo(5);
TA=eo(6);
%Posición y velocidad en el sistema perofocal
rF=A*(1-EC^2)/(1+EC*cos(TA))*[cos(TA),sin(TA),0];
vF=sqrt(mus/A/(1-EC^2))*[-sin(TA),cos(TA)+EC,0];
%Matriz de cambio de ejes
mC=[cos(OM)*cos(W)-sin(OM)*sin(W)*cos(INC), sin(OM)*cos(W)+cos(OM)*sin(W)*cos(INC), sin(W)*sin(INC);-cos(OM)*sin(W)-sin(OM)*cos(W)*cos(INC), -sin(OM)*sin(W)+cos(OM)*cos(W)*cos(INC), cos(W)*sin(INC); sin(OM)*sin(INC), -cos(OM)*sin(INC), cos(INC)];
%Posición y velocidad en el sistema de ejes de la eclíptica (eje X a punto de Aries)
r=mC^-1*rF';
v=mC^-1*vF';
coord=[r;v];
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
