cyl 2D pro - d=z(k)^2+x.^2; gz=k1*z(k)./d; plot(x,gz) hold...

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

View Full Document Right Arrow Icon
Program:- clear all close all clc G=6.67e-8; z=[50:10:1000]'; rho=0.7; A=40; x=[-1000:10:1000]'; k1=2*pi*G*A^2*rho; for k=1:length(z)
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: d=z(k)^2+x.^2; gz=k1*z(k)./d; plot(x,gz) hold on end grid on xlabl('x in meters') ylabl('anomaly')...
View Full Document

This note was uploaded on 05/09/2010 for the course EARTH SCIE APPLIED GE taught by Professor Es during the Spring '09 term at IIT Bombay.

Ask a homework question - tutors are online