重力场大地水准面及重力异常阶次误差和累积误差

%此程序为计算重力场模型的大地水准面和重力异常的阶次误差和累积误差程序，将重力场模型按照（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;