function?gg=qtmxpm（）
%??????均勻球體重力異常--剖面?????%
%球心埋深D（m），球的半徑R（m），剩余密度p（kg/m^3）
%剩余質量M=4/3*pi*R^3*p
%重力異常g=G*M*D/（x^2+D^2）^（3/2）??（m/s^2）
%G=6.67e-11（m^3/kg.s^2）?萬有引力常數
%pi??圓周率
%數據保存在‘均勻球體剖面重力異常.txt’

clear
clc

%球的參數%
d=200;??????%球心埋深D
r=100;???????%球的半徑R
p=2*10^3;????????%剩余密度p
G=6.67e-11;

%計算異常%
x=（-1000:10:1000）;???
m=（4/3）*pi*r^3*p;
n=x.^2+d^2;
gg=G*m*d./（n.^（3/2））*10^5;??%單位mGal

%成圖%
plot（xgg）
xlabel（‘‘）
ylabel（‘‘）
title（‘均勻球體剖面重力異常‘）

%數據生成文本%
y=[x;gg];
fid=fopen（‘均勻球體剖面重力異常.txt‘‘wt‘）;?%wt以文本格式寫入
fprintf（fid‘%4.2f?%.8e\n‘y）;
fclose（fid）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????711??2015-12-01?21:10??qtmxpm.m