为什么循环的结果都一样

flamesmnm

G=6.67*10^(-11);
r=100;
% ?8 L7 u* R+ }# }1 P+ I7 nh=500;
, \; y1 t4 X8 T2 um=0.1*(4/3)*pi*r^3;
, C+ p  x; u  Mg=zeros(41,41);
9 ^! N" P' v4 Z" [2 E" _8 W# M+ Ifor x=-1000:50:1000
    for y=-1000:50:1000
        g(1:41,：)=G*m*h/((x^2+y^2+h^2)^(3/2));
    end
/ m- N% c, J1 ]0 y5 z  cend

modengxian111

你的赋值语句不对，循环中每次赋值都是给g矩阵的整体赋值，所以结果就是最后一次循环给g整体赋值的结果
$ S% }4 L  K" V( |  z0 ?) ox=-1000:50:1000;
y=-1000:50:1000;
6 P5 f/ T9 v2 W% |' [for i=1:41
    for j=1:41
7 e2 `* t% e' t( L: O! ]  p+ v& W0 B        g(i,j)=G*m*h/((x(i)^2+y(j)^2+h^2)^(3/2));
    end
end
复制代码

" {2 p7 t! V, W4 j2 `试试这个