function?pitching_function
clc;?close?all;
%
parameters;
%
param.x?=?0;%?
param.y?=?10000;
param.vel?=?250;
param.theta?=?0;
param.zeta?=?0;
param.omegaZ?=?0;
param.stepsize?=?0.005;
param.time?=?0;
%
param.theta?=?param.theta/180*pi;?%?度轉化為弧度
param.zeta?=?param.zeta/180*pi;
param.alpha?=?param.zeta?-?param.theta;
%
results.x?=?param.x;
results.y?=?param.y;
results.vel?=?param.vel;
results.theta?=?param.theta;
results.zeta?=?param.zeta;
results.omegaZ?=?param.omegaZ;
results.alpha?=?param.alpha;
results.time?=?param.time;
n?=?0;
t?=?cputime;
while?param.y?>??-1e-6
%?while?n?%?????param?=?Runge_Kutta（param）;
????param?=?Adams（param）;
????n?=?n?+?1
????results.x（end?+?1:）?=?param.x;
????results.y（end?+?1:）?=?param.y;
????results.vel（end?+?1:）?=?param.vel;
????results.theta（end?+?