matlab画均匀球体重力异常、水平导数、垂向导数、垂向二阶导数

G =  6.67.*10.^(-11)
M = 200%kg；剩余质量为正
D = 10%m
x = linspace(-100,100,200)

%画的曲线是什么都在title里
figure(1)
delta_g = (G*M*D)./((x.^2 + D.^2).^(3/2))
plot(x,delta_g,'b', 'Linewidth', 2, 'MarkerSize', 8)
xlabel('x')
ylabel('delta-g')
text(0.766*D,(G*M*D)./(((0.766*D).^2 + D.^2).^(3/2)),'\leftarrow X1/2,0.766D')
text(1.035*D,(G*M*D)./(((1.035*D).^2 + D.^2).^(3/2)),'\leftarrow X1/3,1.035D')
title('I类 delta-g of sphere')

figure(2)
delta_g = ((-3)*G*M*D*x)./((x.^2 + D.^2).^(5/2))
plot(x,delta_g,'b', 'Linewidth', 2, 'MarkerSize', 8)
xlabel('x')
ylabel('H-gradient of delta-g')
title('II类 delta-g(H-gradient) of sphere')

figure(3)
delta_g = (G*M*(2*D.^2-x.^2))./((x.^2 + D.^2).^(5/2))
plot(x,delta_g,'b','Linewidth', 2, 'MarkerSize', 8)
hold on;
y = 0*x
plot(x,y,'g'