//example 1.1
//page 8
clc; funcprot(0);
//initialization of variable
P=50*144*47.88;//pressure
R=8323/29;
T=273+20;//temperature
rho=P/R/T;
disp(rho,"density (kg/m^3)=");
clear