//example 1.5 //page 14 clc; funcprot(0); //initialisation of variable T=68+460;//degree R R=1716; kappa=1.4//value from table c=(kappa*T*R)^(0.5); disp(c,"velocity of sound(ft/s)="); clear