% Heat capacity of air in J/(mol K), where T is in K. The
% equation is valid only in the range 273 K < T < 1800 K

for T = 10:500:2500 % loop begins at T=10, increments by 500, ends after 2000

   T % echos value of T

   if T < 273
      disp('Temperature too low')
   elseif T > 1800
      disp('Temperature too high')
   else
      Cp = 28.09 + 0.1965e-2*T + 0.4799e-5*T^2 - 1.965e-9*T^3
      end

  end
