clear;
clc;
global y;
global AMa;
global Aa;
global ACX;
global ACY;
global AP;
global APT;
global AXG;
global AXGT;
global AMZWZ1;
global AMZWZ2;
global AMZAF;
global AJZ;
global AJZT;

load -ascii ACX.txt;    %讀入數據
load -ascii ACY.txt;
load -ascii AP.txt;
load -ascii APT.txt;
load -ascii AJZ.txt;
load -ascii AJZT.txt;
load -ascii AXG.txt;
load -ascii AXGT.txt;
load -ascii Aa.txt;
load -ascii AMa.txt;
load -ascii AMZAF.txt;
load -ascii AMZWZ1.txt;
load -ascii AMZWZ2.txt;

AE=0.01745329;   		%角度化弧度
h=0.005;

fo1=fopen('result.dat','wt');
fprintf(fo1,'y8        y1         y2         y3         y4         y5         y6         y7      ALPHA  \n');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%賦初值
y(1)=20.0;				           %初速
y(2)=18*AE;                        %初始彈道傾角
y(3)=0;				           %初始轉速
y(4)=18*AE;                           %初始俯仰角
y(5)=0;				           %初始x坐標
y(6)=20;				           %初始y坐標
y(7)=52.38;	           %起飛質量
y(8)=0                            %初始時間


fprintf(fo1,'%10.10f %10.10f %10.10f %10.10f %10.10f %10.10f %10.10f %10.10f %10.10f\n',...
             y(8),y(1),y(2)*57.3,y(3)*57.3/360,y(4)*57.3,y(5),y(6),y(7),y(4)-y(2));

while y(6)>=0
    rk(8,h);
    fprintf(fo1,'%10.10f %10.10f %10.10f %10.10f %10.10f %10.10f %10.10f %10.10f %10.10f\n',y(8),y(1),y(2)*57.3,y(3)*57.3/360,y(4)*57.3,y(5),y(6),y(7),y(4)-y(2));
end

fclose(fo1);