--> //initializations

--> e0=1.50 // (initial void ratio)
 e0  = 

   1.5


--> e1=1.44 // (void ratio reduce)
 e1  = 

   1.44


--> Ho=5 //m (thickness of clay stratum)
 Ho  = 

   5.


--> dR=120 //KN/m^2 (effective pressure)
 dR  = 

   120.


--> //calculations

--> de=e0-e1 //(void ratio after reduction)
 de  = 

   0.06


--> mv=-(de/(dR*(1+e0))) // (coefficient of volume change)
 mv  = 

  -0.0002


--> printf('NOTE: negative sign indicates reduction in volume')
NOTE: negative sign indicates reduction in volume
--> Mv=-(mv) //(reduction in volume)
 Mv  = 

   0.0002


--> sf=Mv*dR*Ho*10^3 // (final settlement)
 sf  = 

   120.


--> //results

--> printf('coefficient of volume compressibility =%f',Mv) //m^2/KN
coefficient of volume compressibility =0.000200
--> printf('final settlement of the stratum=%f',sf) //mm
final settlement of the stratum=120.000000
