请问这段程序出了什么问题,本人菜鸟,clear;clc;X=0:1:180;y=0:1:100;r=2; %半径rd=2; %密度rm=(4/3)pir^3rd; %质量x1=75;y1=30;h1=5;x2=95;y2=60;h2=5;x3=85;y3=40;h3=8;x4=105;y4=70;h4=8;G=66.7;for m=1:1:101for n=1:1:181gg1(m,n)=Grm*h1

来源：学生作业帮助网编辑：作业帮时间：2024/05/05 19:21:53

请问这段程序出了什么问题,本人菜鸟,clear;clc;X=0:1:180;y=0:1:100;r=2; %半径rd=2; %密度rm=(4/3)*pi*r^3*rd; %质量x1=75;y1=30;h1=5;x2=95;y2=60;h2=5;x3=85;y3=40;h3=8;x4=105;y4=70;h4=8;G=66.7;for m=1:1:101for n=1:1:181gg1(m,n)=G*rm*h1
请问这段程序出了什么问题,本人菜鸟,
clear;
clc;
X=0:1:180;
y=0:1:100;
r=2; %半径
rd=2; %密度
rm=(4/3)*pi*r^3*rd; %质量
x1=75;
y1=30;
h1=5;
x2=95;
y2=60;
h2=5;
x3=85;
y3=40;
h3=8;
x4=105;
y4=70;
h4=8;
G=66.7;
for m=1:1:101
for n=1:1:181
gg1(m,n)=G*rm*h1/((x(n)-x1)^2+(y(m)-y1)^2+h1^2)^(3/2);
gg2(m,n)=G*rm*h2/((x(n)-x2)^2+(y(m)-y2)^2+h2^2)^(3/2);
gg3(m,n)=G*rm*h3/((x(n)-x3)^2+(y(m)-y3)^2+h3^2)^(3/2);
gg4(m,n)=G*rm*h4/((x(n)-x4)^2+(y(m)-y4)^2+h4^2)^(3/2);
gg(m,n)=gg1(m,n)+gg2(m,n)+gg3(m,n)+gg4(m,n);
end
end
contour(gg);