*/

var('x0');var('y0');var('z0');var('x1');var('y1');var('z1');var('M');var('m');
P0 = vector([x0,y0,z0]);
P1 = vector([x1,y1,z1]);
c = G*M*m;
mag_r = sqrt(((x1-x0)**2)+((y1-y0)**2)+((z1-z0)**2));
r = vector([x1-x0,y1-y0,z1-z0]);
force_mag = c/(mag_r)**2;
G = (6.67*(10**-8));
force_dir = -(r/mag_r);

/* code_ends

Kreyszig-8.4-ex3 (last edited 2010-12-15 12:56:03 by 172)