Planet

  0<a<pi     , a(step)=pi/10

  0<b<2*pi  , b(step)=pi/10

Fx(a,b)=R*cos(b)*sin(a)

Fy(a,b)=R*cos(a)

Fz(a,b)=R*sin(b)*sin(a)

         R=5

satelite of planet

  0<a<pi     , a(step)=pi/10

  0<b<2*pi , b(step)=pi/10

Fx(a,b)=R*cos(b)*sin(a)+4

Fy(a,b)=R*cos(a)+4

Fz(a,b)=R*sin(b)*sin(a)-6

         R=2

Planet

  0<a<pi     , a(step)=pi/10

  0<b<2*pi  , b(step)=pi/10

Fx(a,b)=R*cos(b)*sin(a)

Fy(a,b)=R*cos(a)

Fz(a,b)=R*sin(b)*sin(a)

         R=5

Rings of Planet

  0<a<2*pi  , a(step)=pi/10

  6<r<8       , b(step)=0.5

Fx(a,r)=r*sin(a)

Fy(a,r)=0

Fz(a,r)=r*cos(a)

Sun

  0<a<pi    , a(step)=pi/10

  0<b<2*pi , b(step)=pi/10

Fx(a,b)=R*cos(b)*sin(a)

Fy(a,b)=R*cos(a)

Fz(a,b)=R*sin(b)*sin(a)

         R=5

Earth

  0<a<pi     , a(step)=pi/10

  0<b<2*pi , b(step)=pi/10

  0<g<2*pi , g(step)=pi/6

Fx(a,b)=cos(b)*sin(a)+10*sin(g)-2

Fy(a,b)=cos(a)

Fz(a,b)=sin(b)*sin(a)+8*cos(g)

Moon

  0<a<pi     , a(step)=pi/10

  0<b<2*pi , b(step)=pi/10

  0<g<2*pi , g(step)=pi/6

Fx(a,b)=cos(b)*sin(a)+10*sin(g)-2

Fy(a,b)=cos(a)

Fz(a,b)=sin(b)*sin(a)+8*cos(g)+1

Black Hole 1 (defined with a square plane curve)

  -6<x<6    , x(step)=0.2

  -6<y<6    , y(step)=0.2

Fx(x,y)=x

Fy(x,y)=y

Fz(x,y)=-1/(x^2+y^2)

Black Hole 2 (defined with a cyclical plane curve)

   0<r<6.5      , r(step)=0.1

   0<a<2*pi    , a(step)=pi/40

Fx(x,y)= x=r*cos(a)

Fy(x,y)= y=r*sin(a)

Fz(x,y)=-1/(x^2+y^2)