• API
• FAQ
• Tools
• Archive
SHARE
TWEET

# Untitled

a guest Jun 17th, 2019 52 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
1. f = 6.6743*10^(-11);
2.
3.     m1 = 7.34767309*10^(22);
4.     m2 = 5.97219*10^(24);
5.     e = 0.0549006; a = 384748;
6.     p = a*(1 - e^2);
7.     c = Sqrt[ f*(m1 + m2)*p];
8.     eq1 = {x''[t] == - ((
9.          f*(m1 + m2) *x[t])/(x[t]^2 + y[t]^2 + z[t]^2)^3),
10.             y''[t] == -((
11.          f*(m1 + m2) *y[t])/(x[t]^2 + y[t]^2 + z[t]^2)^3),
12.             z''[t] == -((f*(m1 + m2) *z[t])/(x[t]^2 + y[t]^2 + z[t]^2)^3)};
13.     rul1 = {Subscript[r, j_] -> Sqrt[
14.         Subscript[x, j][t]^2 + Subscript[y, j][t]^2 +
15.          Subscript[z, j][t]^2]};
16.     initial = {x[0] == 1, x'[0] == 0,
17.        y[0] == 0, y'[0] == 4,
18.        z[0] == 0, z'[0] == 3};
19. sol = NDSolve[eq1, initial, {x[t], y[t], z[t]}, {t, 0, 50}];
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.

Top