为什么循环的结果都一样

flamesmnm

G=6.67*10^(-11);
r=100;
( p% a/ e# L5 g; m+ i/ }h=500;
m=0.1*(4/3)*pi*r^3;
$ L: _( x% E) U, c1 S5 J  d1 Dg=zeros(41,41);
for x=-1000:50:1000
2 j) M3 X: ]: t' Z( z    for y=-1000:50:1000
9 O% z/ R* ]5 A; P. x- {        g(1:41,：)=G*m*h/((x^2+y^2+h^2)^(3/2));
    end
4 B8 m6 e/ i: `1 Qend

modengxian111

你的赋值语句不对，循环中每次赋值都是给g矩阵的整体赋值，所以结果就是最后一次循环给g整体赋值的结果
x=-1000:50:1000;
y=-1000:50:1000;
for i=1:41
    for j=1:41
        g(i,j)=G*m*h/((x(i)^2+y(j)^2+h^2)^(3/2));
* I1 [$ i% J2 K, c) h# A    end
end
复制代码

试试这个