function?[ty]?=?RK4（ftbegantendy1h）
n=（tend-tbegan）/h;?
y（:1）=y1;
t（1）=tbegan;
global?V;???????%導(dǎo)彈速度
global?qd;??????%λd?撞擊角度
global?td;??????%撞擊時間
global?flag;
flag=0;
i=1;
while（i<=n&&t（i）????K1=h*f（t（i）y（1:4i））;
????K2=h*f（t（i）+h/2y（1:4i）+K1/2）;
????K3=h*f（t（i）+h/2y（1:4i）+K2/2）;
????K4=h*f（t（i）+hy（1:4i）+K3）;
????t（i+1）=t（i）+h;
????y（1:4i+1）=y（1:4i）+（K1+2*K2+2*K3+K4）/6;
????%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
????y（5i+1）=td-t（i）-y（1i）;
????if（y（5i+1）<0.0000001）
???????flag=1;
????end
????r=y（1i）*V;
????lamuda=y（2i）+qd;
????y（6i+1）=-r*cos（lamuda）;
????y（7i+1）=-r*sin（lamuda）;
????y（8i+1）=lamuda;
????y（9i+1）=V*（（y（3i+1）-y（3i））/h-V*sin（y（3i））/r）;
????
????i=i+1;???
end

end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????798??2018-12-02?11:03??code\RK4.m

?????文件????????674??2018-12-02?11:00??code\solvefunc.m

?????文件???????1262??2018-12-02?11:03??code\solver.m

?????文件??????51128??2018-12-02?10:37??code\td=55s，λd改變，a變化曲線.fig

?????文件??????49072??2018-12-02?10:37??code\td=55s，λd改變，λ變化曲線.fig

?????文件??????50347??2018-12-02?10:37??code\td=55s，λd改變，σ變化曲線.fig

?????文件??????58204??2018-12-02?10:29??code\td=55s，λd改變，彈道圖.fig

?????文件??????55334??2018-12-02?10:56??code\λd=-120°，td改變，a變化曲線.fig

?????文件??????55292??2018-12-02?10:53??code\λd=-120°，td改變，λ變化曲線.fig

?????文件??????54307??2018-12-02?10:54??code\λd=-120°，td改變，σ變化曲線.fig

?????文件??????64458??2018-12-02?10:52??code\λd=-120°，td改變，彈道圖.fig

?????目錄??????????0??2019-09-02?11:24??code

-----------?---------??----------?-----??----

???????????????440876????????????????????12