function[X] = COE2RV(coe, mu);
r = (coe(1)*(1-coe(2)^2))/(1+coe(2)*cos(coe(6)));
h = sqrt(mu*coe(1)*(1-norm(coe(2))^2));
IP = [cos(coe(4))*cos(coe(5))-sin(coe(4))*cos(coe(3))*sin(coe(5)), -cos(coe(4))*sin(coe(5))-sin(coe(4))*cos(coe(3))*cos(coe(5)), sin(coe(4))*sin(coe(3)); sin(coe(4))*cos(coe(5))+cos(coe(4))*cos(coe(3))*sin(coe(5)), -sin(coe(4))*sin(coe(5))+cos(coe(4))*cos(coe(3))*cos(coe(5)), -cos(coe(4))*sin(coe(3)); sin(coe(3))*sin(coe(5)), sin(coe(3))*cos(coe(5)), cos(coe(3))];
R = IP*[r*cos(coe(6)); r*sin(coe(6)); 0;];
v = IP*[(-mu/h)*sin(coe(6)); (mu/h)*(cos(coe(6))+coe(2)); 0;];
[X] = [R; v];
end