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. OK, I Understand
 
Top