# what to do when array indicies exceed for , FOR loop comparision for below , trying to plot density of atmosphere, how to do the comparison ?

1 view (last 30 days)
SL on 24 Sep 2021
Commented: Walter Roberson on 25 Sep 2021
h(1) = 0;
T(1) = 273.15;
p(1) = 101.2;
for t = 1:328084
if h(t)<[11000]
T(t+1) = 15.04 - 0.00649.*h(t);
p(t+1) = 101.29 .*((T(t) + 273.15) /288.08).^5.256;
elseif [11000]> h(t) && h(t)<=[25000]
T(t) = -56.46;
p(t+1) = 22.65.*e^(1.73-0.000157.*h(t));
elseif h(t)>[25000]
T(t) = -131.21+0.00299.*h(t);
p(t+1) = 2.488 .* ((T(t) + 273.15)/216.6).^-11.388;
end
rho_h(t) = p(t) ./(0.2869 * (T(t) + 273.15))
end
figure(1)
plot(T,h)
grid on
legend on
figure(2)
plot(rho_h, h)
grid on
legend on

Walter Roberson on 24 Sep 2021
h(1) = 0;
That is the only place you assign to h
if h(t)<[11000]
but you need h to be 328084 elements by the end of the loop.
SL on 25 Sep 2021
Edited: SL on 25 Sep 2021
perfect thank you very much highly appreciate it
Walter Roberson on 25 Sep 2021
please recheck rho_h(1) which is not assigned a value.

### Categories

Find more on MATLAB in Help Center and File Exchange

R2021b

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!