% function check_main()
clc
clear
close all

[AOC_t0,AOC_fai,AOC_thita,AOC_pusai,AOC_a0,AOC_e0,AOC_i0,AOC_OMG0,AOC_os0,AOC_u0,AOC_Tstep] = readInput();

AOC_dTime=86400;
AOC_Tin=0:AOC_Tstep:AOC_dTime;
Mu=3.986005e5;                      %地球引力常数
w0=(Mu/AOC_a0^3)^0.5;               %卫星轨道角速度
we=7.29211585530487e-5;             %地球自转角速度
Re=6371.23;                         %地球平均半径
Len=length(AOC_Tin);
Thita=linspace(-8.6,8.6,Len);
for num=1:Len
    Tin=AOC_Tin(num);
    AOC_u=AOC_u0+we*Tin;
    Rp_b=[sind(Thita(num)),0,cosd(Thita(num))].';
    [Tmp_Lon,Tmp_Lat]=SightPoint(Tin,AOC_fai,AOC_thita,AOC_pusai,AOC_a0,AOC_e0,AOC_i0,AOC_OMG0,AOC_os0,AOC_u,Rp_b);
    Lon(num,1)=Tmp_Lon;
    Lat(num,1)=Tmp_Lat;
end

saveData(Lon,)