;Edward Burin des Roziers
;Asen5050
;9/24/03

function k2c, peiOmegawnu

p = peiOmegawnu[0]
e = peiOmegawnu[1]
i = peiOmegawnu[2]
Omega = peiOmegawnu[3]
w = peiOmegawnu[4]
nu = peiOmegawnu[5]

mu = 3.986d5

pospqw = dblarr(3)
velpqw = dblarr(3)
ijk_pqw = dblarr(3,3)

x = 1 + e * cos(nu)
pospqw = [(p*cos(nu))/x,(p*sin(nu))/x,0]
velpqw = [-1*sqrt(mu/p)*sin(nu),sqrt(mu/p)*(e+cos(nu)),0]

posijk = (rot3(-Omega) ## (rot1(-i) ## (rot3(-w) ## pospqw)))
velijk = (rot3(-Omega) ## (rot1(-i) ## (rot3(-w) ## velpqw)))

posvel = [posijk[0],posijk[1],posijk[2],velijk[0],velijk[1],velijk[2]]
return, posvel

end