为什么循环的结果都一样

flamesmnm

G=6.67*10^(-11);
& v# ?8 R* j: b" u0 I9 Br=100;
h=500;
% M1 B$ }! i8 z0 z: A4 @- m6 @m=0.1*(4/3)*pi*r^3;
g=zeros(41,41);
8 C, F3 M( G0 b- W- gfor x=-1000:50:1000
* T  L6 O# V' d% _% n; x- E3 A    for y=-1000:50:1000
: e  ^( j8 Z$ a2 i) q. q7 `, C        g(1:41,：)=G*m*h/((x^2+y^2+h^2)^(3/2));
    end
end
2 E9 L5 ?0 D* V& z- m2 a* J

modengxian111

你的赋值语句不对，循环中每次赋值都是给g矩阵的整体赋值，所以结果就是最后一次循环给g整体赋值的结果
0 F( Q7 U0 L1 I, e6 px=-1000:50:1000;
y=-1000:50:1000;
  a8 i/ |0 U+ ]3 b3 }# ^; }for i=1:41
8 }# Z9 L4 J6 D" T, M8 l! g    for j=1:41
        g(i,j)=G*m*h/((x(i)^2+y(j)^2+h^2)^(3/2));
    end
end
6 w; e! p' n$ y3 D/ V复制代码

试试这个