Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- cols = RGBColor /@ {"#2537b5", "#ff3804", "#ffffff", "#000000"};
- a = 384748;
- rp = 362600;
- ra = 405400;
- e = (ra - rp)/(ra + rp);
- b = a Sqrt[1 - e^2];
- ea = 146600000;
- ee = .01671123;
- eb = ea Sqrt[1 - ee^2];
- er = 6371;
- mr = 1737;
- sr = 696342;
- In[15]:= Manipulate[
- Module[{earthOrbitScaleFactor, moonOrbitScaleFactor,
- earthRadiusScaleFactor, moonRadiusScaleFactor},
- earthOrbitScaleFactor = .05;
- moonOrbitScaleFactor = 40;
- earthRadiusScaleFactor = 40;
- moonRadiusScaleFactor = 40;
- Graphics[{cols[[2]], Disk[{0, 0}, sr], cols[[1]],
- Disk[earthOrbitScaleFactor {eb Cos[\[Theta]], ea Sin[\[Theta]]},
- earthRadiusScaleFactor*er], White,
- Disk[earthOrbitScaleFactor {eb Cos[\[Theta]] +
- moonOrbitScaleFactor *a Cos[13 \[Theta]],
- ea Sin[\[Theta]] + moonOrbitScaleFactor*b Sin[13 \[Theta]]},
- moonRadiusScaleFactor*mr]}, Background -> Black,
- ImageSize -> {540, 540},
- PlotRange -> 1.5 earthOrbitScaleFactor * ea]
- ], {\[Theta], 0., 2 \[Pi]}]
- Out[15]= Manipulate[Module[{earthOrbitScaleFactor$, \
- moonOrbitScaleFactor$,
- earthRadiusScaleFactor$, moonRadiusScaleFactor$},
- earthOrbitScaleFactor$ = 0.05; moonOrbitScaleFactor$ = 40;
- earthRadiusScaleFactor$ = 40; moonRadiusScaleFactor$ = 40;
- Graphics[{cols[[2]], Disk[{0, 0}, sr], cols[[1]],
- Disk[earthOrbitScaleFactor$*{eb*Cos[\[Theta]], \
- ea*Sin[\[Theta]]},
- earthRadiusScaleFactor$*er], White,
- Disk[earthOrbitScaleFactor$*{eb*Cos[\[Theta]] + \
- moonOrbitScaleFactor$*a*
- Cos[13*\[Theta]], ea*Sin[\[Theta]] + \
- moonOrbitScaleFactor$*b*Sin[13*\[Theta]]},
- moonRadiusScaleFactor$*mr]}, Background -> Black,
- ImageSize -> {540, 540}, PlotRange -> \
- 1.5*earthOrbitScaleFactor$*ea]],
- {{\[Theta], 6.283185307179586}, 0., 2*Pi}]
Advertisement
Add Comment
Please, Sign In to add comment