简单弹道方程MATLAB代码

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 <  32146
%     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 +