// Program To calculate properties of dry air from 223 K to 373 K
clear all
clc
printf("Temp \t\t Density \t Dynamic.Viscosity \tKinematic.Viscosity \tThermal conductivity \tspecific heat \t Thermal dissusivity")
for t=-50:0.5:100
tfk=t+273         //temperature in kelvin
pf1=351.99/tfk
pf2=344.84/(tfk**2)
pf=pf1+pf2  //Density
x=3/2
y=5/2
dyvisc1=1.4592
dyvisc11=tfk**x
dyvisc13=dyvisc1*dyvisc11
dyvisc2=109.10+tfk
dyvisc3=dyvisc13/dyvisc2
dyvisc=dyvisc3*(1e-6)   //Dynamic viscosity
ka1=2.3340e-3*(tfk**1.5)
ka2=164.54+tfk
ka=ka1/ka2          //Thermal conductivity
cpa1=1030.5-(0.19975*tfk)
cpa2=3.9734e-4*(tfk**2)
cpa=cpa1+cpa2       //Specific heat
kvisc1=tfk**x
kvisc11=2.4090e8/kvisc1
kvisc2=tfk**y
kvisc22=2.6737e10/kvisc2
kvisc3=kvisc11+kvisc22
kvisc=1/kvisc3        //Kinematic viscosity
alph1=(1.5556e-4)*(tfk**2)
alph2=(4.1190e-2)*(tfk)
alph3= -4.3274+alph1+alph2
alpha=alph3*(1e-6)       //Thermal Diffusivity
printf("\t\n%0.1f C",t)
printf("\t\t%0.3f Kgm^3",pf)
printf("\t%0.9f Nsm^2",dyvisc)
printf("\t%0.9f m^2/s",kvisc)
printf("\t%f W/mK",ka)
printf("\t\t%d J/KgK",cpa)
printf("\t%0.9f m^2/s",alpha)
end