//example 2.5 //page 57 clc; funcprot(0); //initialisation of variable Gamma=9810; ha=950*1.5-500*1; ha=ha/1000;//m of H2O; Pa=ha*Gamma; disp(Pa/1000,"Pressure at A (kPa)="); clear