sph 2D PRO - plot(x,gz) hold on grid on xlabel('x in meter')

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
G=6.67e-8; R=50; rho=2.8; z=[50:5:500]'; x=[-1000:10:1000]'; k1=4*pi*R^3*rho*G/3; for k=1:length(x) d=(z(k).^2+x.^2).^(3/2); gz=k1*z(k)./d;
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: plot(x,gz) hold on grid on xlabel('x in meter') ylabel('anomaly') title('gravity anomaly in 2-D') end...
View Full Document

Ask a homework question - tutors are online