%此程序為計算重力場模型的大地水準面和重力異常的階次誤差和累積誤差程序，將重力場模型按照（gfc，LMCSsigmaCsigmaS）的順序排
%列好并命名為“input.txt”進行計算。最后輸出的結果依次為大地水準面階次誤差及累積誤差和重力異常階次誤差和累積誤差（均從第2階次開始不為零），程序可能存在諸多問題，后續改進。
[a1LMCSsigmaCsigmaS]=textread（‘input.txt‘‘%s?%n?%n?%f?%f?%f?%f‘）;
T=sigmaC.*sigmaC+sigmaS.*sigmaS;
A=[LMT];
[ab]=size（A）;
for?i=1:a
????B（int16（A（i1））+1int16（A（i2））+1）=A（i3）;
end
[xy]=size（B）;
for?h=1:x
????R=6378136.46;
????GM=398600441500000;
????Q=0;
????for?m=1:h
????????Q=Q+B（hm）;
????end
????P（h1）=R*sqrt（Q）;
????G（h1）=GM*（h-2）/（R^2）*sqrt（Q）;
end
P2=P.*P;
G2=G.*G;
Pl（21）=0Gl（21）=0;
for?i2=3:x
????Pl（i21）=Pl（i2-11）+P2（i21）;
????Gl（i21）=Gl（i2-11）+G2（i21）;
????PL（i21）=sqrt（Pl（i21））;
????GL（i21）=sqrt（Gl（i21））;
end
RA=[PPLGGL];
save?ra.txt?-ascii?RA;