为什么循环的结果都一样

flamesmnm

G=6.67*10^(-11);
r=100;
  h5 I1 ]1 p  Th=500;
; I3 Q0 T, D1 p( ~4 I0 B5 g* @; {m=0.1*(4/3)*pi*r^3;
3 \& l3 W5 }3 Ug=zeros(41,41);
for x=-1000:50:1000
    for y=-1000:50:1000
        g(1:41,：)=G*m*h/((x^2+y^2+h^2)^(3/2));
    end
/ A3 [& Z6 F" mend
+ k" h  r9 _) J1 w; M1 d

modengxian111

你的赋值语句不对，循环中每次赋值都是给g矩阵的整体赋值，所以结果就是最后一次循环给g整体赋值的结果
; b/ P/ ?; O5 a, _# b& u( Cx=-1000:50:1000;
y=-1000:50:1000;
for i=1:41
! N; a# H$ E- [" d    for j=1:41
        g(i,j)=G*m*h/((x(i)^2+y(j)^2+h^2)^(3/2));
% J2 h, E2 J1 h, O$ R* L: r. j8 Z    end
end
& q( P! t, O/ C4 n3 n& g+ `复制代码
) H2 F3 R4 _! W- B6 j& n$ y
) n8 \/ L  r8 c" M  P6 C试试这个