--> //initializations

--> G=2.67 //(specific gravity)
 G  = 

   2.67


--> w=0.15 //(water content)
 w  = 

   0.15


--> Yw=1.0 //g/ml (unit weight of water)
 Yw  = 

   1.


--> Yd=1.78 //g/ml (dry density)
 Yd  = 

   1.78


--> //calculations

--> s=(Yd*w*G*100)/((G*Yw)-Yd) // (degree of saturation)
 s  = 

   80.1


--> na=(1-((Yd+Yd*w*G)/(G*Yw)))*100 // (air voids)
 na  = 

   6.6333333


--> yd=(G*Yw)/(1+w*G) // (dry density)
 yd  = 

   1.906462


--> //results

--> printf('degree of saturation=%f',s) //%
degree of saturation=80.100000
--> printf('air voids=%f',na) //%
air voids=6.633333
--> printf('value of dry density on the zero air void line=%f',yd) //g/ml
value of dry density on the zero air void line=1.906462
