//example 1.10 //page 29 clc; funcprot(0); //initialisation of variable d=0.01;//diameter sigma=0.073;//from table delP=4*sigma/d; disp(delP,"excess pressure in (N/m^2)"); clear