为什么循环的结果都一样

flamesmnm

G=6.67*10^(-11);
# Z+ H( I% l% }& C4 L  xr=100;
. f& f! n, K1 R- j& \. N/ m* Kh=500;
m=0.1*(4/3)*pi*r^3;
g=zeros(41,41);
for x=-1000:50:1000
- u  R% K6 v+ O1 B7 E) p, Q    for y=-1000:50:1000
        g(1:41,：)=G*m*h/((x^2+y^2+h^2)^(3/2));
    end
end

modengxian111

你的赋值语句不对，循环中每次赋值都是给g矩阵的整体赋值，所以结果就是最后一次循环给g整体赋值的结果
x=-1000:50:1000;
( T2 P$ t, @( ]& ~; q$ {+ my=-1000:50:1000;
for i=1:41
2 ?$ A' s2 ^: Y$ x- N- Z& k9 |0 `    for j=1:41
7 [% r( t2 E( E. Q# W' J: Z  H, }- |        g(i,j)=G*m*h/((x(i)^2+y(j)^2+h^2)^(3/2));
# }+ H& C5 w: O2 ?    end
end
' z3 Z; j9 t5 Y! {复制代码
; e% t: o9 B$ r+ B' C# ]( b! W
试试这个