为什么循环的结果都一样

flamesmnm

G=6.67*10^(-11);
3 e4 t0 y! E/ s) [. O5 r5 i' Jr=100;
h=500;
; M# l8 i; M! \1 o( a/ u/ [  C# Am=0.1*(4/3)*pi*r^3;
( H7 j; V  O, t+ B( U) H! |6 O1 A& Dg=zeros(41,41);
) i- W# K8 d  afor x=-1000:50:1000
    for y=-1000:50:1000
        g(1:41,：)=G*m*h/((x^2+y^2+h^2)^(3/2));
! x5 F) |# D7 z+ w    end
% X  P1 E' K4 d% D7 M, ^2 wend

modengxian111

你的赋值语句不对，循环中每次赋值都是给g矩阵的整体赋值，所以结果就是最后一次循环给g整体赋值的结果
( h! i7 `1 K) a4 Dx=-1000:50:1000;
* }( V2 M- C# X( b3 G4 Q/ ~0 _y=-1000:50:1000;
6 K* d2 |/ o4 ^. K+ c4 Ufor i=1:41
; d! @1 h0 h5 v) g    for j=1:41
        g(i,j)=G*m*h/((x(i)^2+y(j)^2+h^2)^(3/2));
    end
end
复制代码
4 Z3 F, X1 b( Z$ U
试试这个