为什么循环的结果都一样

flamesmnm

G=6.67*10^(-11);
r=100;
) N$ l  E0 f$ m# Lh=500;
m=0.1*(4/3)*pi*r^3;
g=zeros(41,41);
for x=-1000:50:1000
$ E8 l* `4 h$ b$ J& |" A  o( X    for y=-1000:50:1000
        g(1:41,：)=G*m*h/((x^2+y^2+h^2)^(3/2));
    end
end
3 P8 @6 B+ f+ T* N" ]% l7 H

modengxian111

你的赋值语句不对，循环中每次赋值都是给g矩阵的整体赋值，所以结果就是最后一次循环给g整体赋值的结果
x=-1000:50:1000;
. x. A- \( `, I% M% U# A" }y=-1000:50:1000;
. Y4 X2 V; {6 J2 vfor i=1:41
    for j=1:41
        g(i,j)=G*m*h/((x(i)^2+y(j)^2+h^2)^(3/2));
    end
end
复制代码

6 ]' r" _- q( P. A- B8 Q- a试试这个