为什么循环的结果都一样

flamesmnm

G=6.67*10^(-11);
- z  F5 P3 j2 U8 T- Vr=100;
% H, k0 L$ b9 K$ ^3 K* ih=500;
m=0.1*(4/3)*pi*r^3;
) w5 _! v+ [5 [9 n2 Lg=zeros(41,41);
for x=-1000:50:1000
    for y=-1000:50:1000
' E6 Y7 h8 b! [1 K1 y$ M1 K        g(1:41,：)=G*m*h/((x^2+y^2+h^2)^(3/2));
    end
9 ]' l+ S2 w) _; W/ Aend
' t2 ?5 `+ K$ G, V2 G0 L% i

modengxian111

你的赋值语句不对，循环中每次赋值都是给g矩阵的整体赋值，所以结果就是最后一次循环给g整体赋值的结果
x=-1000:50:1000;
y=-1000:50:1000;
1 m7 p8 c, H, N( A4 B1 G7 E: a' v6 Rfor i=1:41
# o1 b: V1 j/ C    for j=1:41
- x; k# v& \( ?/ p        g(i,j)=G*m*h/((x(i)^2+y(j)^2+h^2)^(3/2));
    end
end
. p& _5 k! `" V1 c9 f: a$ O& M# @复制代码

) |3 c) ~# s7 ]: {! U' z' v4 x; _试试这个