Need help on program

Hello,

I am trying to model the below conditions. Input h is given as 0:50:25000

rho value should not cross 1.23 (when h=0) and it must decrease as h increases. But I am getting rho value of 7+ which is incorrect (Refer graph).

Can someone help me on this?

,

Below is the script I have written for the above condition

for i = 1:length(h)

if h(i)<11000

Te = 15.04-0.00649.*h;

pe = 101.29.*((Te+273.1)./288.08).^5.256;

elseif h(i)>=11000 & h(i)<=25000

Te = -56.46;

pe = 22.65.*10.^(1.73-0.000157.*h);

else

Te = -131.21+(0.00299.*h);

pe = 2.488*((Te+273.1)/216.6).^(-11.388);

end

end

rhoe = pe./(0.2869.*(Te+273.1));

plot(h,rhoe,'Color',[0 0 1],'LineWidth',1.5);

Answer 1

[Torsten]

h = 0:50:50000;
for i = 1:length(h)
    if h(i)<11000
        Te(i) = 15.04-0.00649.*h(i);
        pe(i) = 101.29.*((Te(i)+273.1)./288.08).^5.256;
    elseif h(i)>=11000 && h(i)<=25000
        Te(i) = -56.46;
        pe(i) = 22.65.*exp(1.73-0.000157.*h(i));
    else
        Te(i) = -131.21+(0.00299.*h(i));
        pe(i) = 2.488*((Te(i)+273.1)/216.6).^(-11.388);
    end
end
rhoe = pe./(0.2869.*(Te+273.1));
plot(h,rhoe,'Color',[0 0 1],'LineWidth',1.5);

[Chinmayraj Doddarajappa]

Thank you