%右函數定義的，用于單步龍格庫塔法計算。
%如果需要使用全局變量，要使用global定義。
% 這種格式的，也完全可以被ode45函數調用。

function dy = rightF(t,y) % right functions
global XNC VM;

dy = zeros(5,1);    % a column vector

dy(1) = XNC/VM;
dy(2) = -XNC*sin(y(1)); %過載垂直于theta（彈道傾角）
dy(3) =1*XNC*cos(y(1));

dy(4) =y(2);
dy(5) =y(3);

% 下面的，為方便調用程序，把stepRK.m的源代碼放在這

%函數說明，t表示自變量（一般為時間），y表示變量的初值（是一個行向量），h表示步長
% %運行的結果，也是一個行向量
% 
% function dy = stepRK(t,y,h)
% K1=rightF(t,y);
% K2=rightF(t+h/2,y'+h*K1/2);
% K3=rightF(t+h/2,y'+h*K2/2);
% K4=rightF(t+h,y'+h*K3);
% dy=y'+(K1+2*K2+2*K3+K4)*h/6;
% dy=dy';
%