为什么循环的结果都一样

flamesmnm

G=6.67*10^(-11);
r=100;
h=500;
6 ?5 y9 I, f/ X4 H! N  c7 Cm=0.1*(4/3)*pi*r^3;
3 A" O& z6 w. n+ h- l$ u5 Ng=zeros(41,41);
for x=-1000:50:1000
' k1 w' x" ^# n' f' q/ y0 O    for y=-1000:50:1000
( c& ?% n( l% f6 b; {* }        g(1:41,：)=G*m*h/((x^2+y^2+h^2)^(3/2));
    end
end

modengxian111

你的赋值语句不对，循环中每次赋值都是给g矩阵的整体赋值，所以结果就是最后一次循环给g整体赋值的结果
x=-1000:50:1000;
% V) e/ r9 g6 L% s  |6 O& M# f" {/ Uy=-1000:50:1000;
" l$ e" \% N  Y. M# l* q. {1 x* [for i=1:41
    for j=1:41
; s0 C4 q# r% I/ \; |        g(i,j)=G*m*h/((x(i)^2+y(j)^2+h^2)^(3/2));
    end
+ A9 r+ l, v# m. E7 o% tend
复制代码

5 U8 S, M: y9 N! |! b: o试试这个