大气密度函数的m文件

function Ro=density（r）
%大气密度叉乘公式
%注意单位为千米
%大气层高度设定为1000Km
%r 飞行器在大气中的高度
r=r/1000;
if r<0
 Ro=1.225*exp（-0.14*r）;
    return
else
    if r>=1000
    Ro=0;
    else
    H=[0.0000 0.5000 1.0000 1.5000 2.0000 2.5000 3.0000 3.5000 4.0000 4.5000 5.0000 5.5000 6.0000 6.5000 7.0000 7.5000 8.0000 8.5000 9.0000 9.5000...
   10.000 10.500 11.000 12.000 13.000 14.000 15.000 16.000 17.000 18.000 19.000 20.000 21.000 22.000 23.000 24.000 25.000 26.000 27.000 28.000...
   29.000 30.000 31.000 32.000 34.000 36.000 38.000 40.000 42.000 44.000 46.000 48.000 50.000 55.000 60.000 65.000 70.000 75.000 80.000 85.000...
   90.000 95.000 100.00 110.00 120.00 130.00 140.00 150.00 160.00 170.00 180.00 190.00 200.00 210.00 220.00 230.00 240.00 250.00 260.00 270.00...
   280.00 290.00 300.00 320.00 3