G=6.67*10^(-11);/ A! m C( }/ D
r=100;3 b- [% n1 j3 h3 y }) H8 I
h=500;- ^5 c* T7 F* [( \
m=0.1*(4/3)*pi*r^3;& C# m% b# }) j, a2 W9 ^
g=zeros(41,41); 8 H; ]" T3 s- ]8 P [4 m) [ }for x=-1000:50:10001 _. ^: {- o- ]8 Z Z
for y=-1000:50:1000 : Y [4 P, p @& d g(1:41,:)=G*m*h/((x^2+y^2+h^2)^(3/2));0 h7 T( q9 v6 C$ r
end ) T _/ D* o% P5 I! l. Eend ) O |9 v+ `( z1 N6 k2 Q5 l+ N$ V