r=Function[{a,e,ϕ},a(1-e^2)/(1+e Cos[ϕ])] (* Standardowa trajektoria eliptyczna *)

heliocentric=Function[{a,e,Ω,ω,i,ϕ},{r[a,e,ϕ] (Cos[Ω]Cos[ω+ϕ]-Sin[Ω] Sin[ω+ϕ]Cos[i]),r[a,e,ϕ] (Sin[Ω]Cos[ω+ϕ]+Cos[Ω]Sin[ω+ϕ]Cos[i]),r[a,e,ϕ] Sin[ω+ϕ]Sin[i]}];

Manipulate[Show[ParametricPlot3D[
heliocentric[a,e,Ω,ω,i,ϕ],{ϕ,0,2Pi},PlotRange->6{{-1,1},{-1,1},{-1,1}}],Graphics3D[{Opacity[0.3],
InfinitePlane[{
heliocentric[a,e,Ω,ω,i,0],
heliocentric[a,e,Ω,ω,i,Pi/2],
heliocentric[a,e,Ω,ω,i,Pi]
}],InfinitePlane[{{0,0,0},{0,1,0},{1,0,0}}]}]],
{{a,2},0,12},
{{e,8/9},0,1},
{{i,Pi/6},-Pi,Pi},{{ω,1},-Pi,Pi},{{Ω,1},-Pi,Pi}]