function?[?PV?]?=?orbitdynamics_runge_kutta（?PVh?）

x=P（1）;
y=P（2）;
z=P（3）;
vx=V（1）;
vy=V（2）;
vz=V（3）;

K11=FC1（vx）;
K12=FC2（vy）;
K13=FC3（vz）;
K14=FC4（xyz）;
K15=FC5（xyz）;
K16=FC6（xyz）;
?
x=P（1）+（h/2）*K11;
y=P（2）+（h/2）*K12;
z=P（3）+（h/2）*K13;
vx=V（1）+（h/2）*K14;
vy=V（2）+（h/2）*K15;
vz=V（3）+（h/2）*K16;
?
K21=FC1（vx）;
K22=FC2（vy）;
K23=FC3（vz）;
K24=FC4（xyz）;
K25=FC5（xyz）;
K26=FC6（xyz）;

x=P（1）+（h/2）*K21;
y=P（2）+（h/2）*K22;
z=P（3）+（h/2）*K23;
vx=V（1）+（h/2）*K24;
vy=V（2）+（h/2）*K25;
vz=V（3）+（h/2）*K26;
??
K31=FC1（vx）;
K32=FC2（vy）;
K33=FC3（vz）;
K34=FC4（xyz）;
K35=FC5（xyz）;
K36=FC6（xyz）;


x=P（1）+h*K31;
y=P（2）+h*K32;
z=P（3）+h*K33;
vx=V（1）+h*K34;
vy=V（2）+h*K35;
vz=V（3）+h*K36;

K41=FC1（vx）;
K42=FC2（vy）;
K43=FC3（vz）;
K44=FC4（xyz）;
K45=FC5（xyz）;
K46=FC6（xyz）;

P（1）?=?P（1）+h/6*（K11+2*K21+2*K31+K41）;
P（2）?=?P（2）+h/6*（K12+2*K22+2*K32+K42）;
P（3）?=?P（3）+h/6*（K13+2*K23+2*K33+K43）;
V（1）?=?V（1）+h/6*（K14+2*K24+2*K34+K44）;
V（2）?=?V（2）+h/6*（K15+2*K25+2*K35+K45）;
V（3）?=?V（3）+h/6*（K16+2*K26+2*K36+K46）;

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function?[K]=FC1（vx）
K=vx;
end
?
function?[K]=FC2（vy）
K=vy;
end
?
function?[K]=FC3（vz）
K=vz;
end
?
function?[K]=FC4（xyz）
u=398600.44;
r?=?sqrt（x^2+y^2+z^2）;
K?=?-u*x/r^3;
end
?
function?[K]=FC5（xyz）
u=398600.44;
r?=?sqrt（x^2+y^2+z^2）;
K?=?-u*y/r^3;
end
?
function?[K]=FC6（xyz）
u=398600.44;
r?=?sqrt（x^2+y^2+z^2）;
K?=?-u*z/r^3;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-11-14?11:56??matlab實現\
?????目錄???????????0??2018-03-07?23:29??matlab實現\離心率為0\
?????目錄???????????0??2018-03-07?23:29??matlab實現\離心率為0.6\
?????文件????????1527??2017-10-06?17:06??matlab實現\離心率為0.6\orbitdynamics_runge_kutta.m
?????文件????????1835??2017-10-06?17:26??matlab實現\離心率為0.6\orbitdynamic_main.m
?????文件?????8018551??2017-09-28?14:59??matlab實現\離心率為0.6\Satellite1?J2000?Position?Velocity.csv
?????文件?????3517019??2017-10-06?17:17??matlab實現\離心率為0.6\satelliteposition.mat
?????文件????????8348??2017-09-28?14:49??matlab實現\離心率為0.6\捕獲.PNG
?????文件????????1527??2017-10-06?17:23??matlab實現\離心率為0\orbitdynamics_runge_kutta.m
?????文件????????1570??2017-10-06?17:24??matlab實現\離心率為0\orbitdynamic_main.m
?????文件?????3820122??2017-09-28?10:57??matlab實現\離心率為0\STK_real.mat
?????文件?????5696752??2017-09-28?10:40??matlab實現\離心率為0\輸入初始位置速度.csv
?????目錄???????????0??2018-11-14?11:55??python實現\
?????文件?????????749??2018-11-08?14:56??python實現\draw?a?ball.py
?????文件????????2248??2018-11-08?14:56??python實現\runge_kutta.py
?????文件???????80136??2018-11-11?22:04??python實現\x_array.npy
?????文件???????80136??2018-11-11?22:04??python實現\y_array.npy
?????文件???????80136??2018-11-11?22:04??python實現\z_array.npy