探空火箭优化设计

箭头60kg，总冲量400KN·s，总推进剂180kg，总质量500kg，请编程设计二级探空火箭并优化设计一二级分配方案，使火箭飞行高度最大。

程序编的比较粗糙，仅供参考。

clear
clc

Cd=0.2;%阻力系数
Cl=0.1;%升力系数

alpha=0; %攻角

S=0.1; %参考面积，m^2

g=9.8; %重力加速度
row=1; %大气密度

Dt=0.001;
n1=1;

for m1=40:20:110  %一级推进剂质量，kg
    m2=180-m1; %二级推进剂质量，kg
    for I1=100000:20000:300000   %一级总比冲
        I2=400000-I1;  %二级总比冲
        for F1=(500*g+100):500:20000
            t1=I1/F1;
            for F2=((500-m1-20)*g+100):500:20000
                t2=I2/F2;
                
                x=0;
                y=0;
                v=35;
                t=0;
                theta=90*pi/180; %弹道倾角
                flag=0;