--> //initializations

--> Ho=2 //m (thickness of clay layer)
 Ho  = 

   2.


--> Cr=0.05 // (coefficient of recompression)
 Cr  = 

   0.05


--> Cc=0.25 // (coefficient of compressibility)
 Cc  = 

   0.25


--> eo=1.4 // (void ratio at clay layer)
 eo  = 

   1.4


--> Rc=75 //KN/m^2 (preconsolidation pressure)
 Rc  = 

   75.


--> Ro=50 //KN/m^2 (overburden pressure)
 Ro  = 

   50.


--> dR=40 //KN/m^2 (increase in pressure)
 dR  = 

   40.


--> //calculations

--> sf=(((Cr*Ho)/(1+eo))*((log(Rc/Ro))/2.303))+(((Cc*Ho)/(1+eo))*((log((Ro+dR)/Rc))/2.303))
 sf  = 

   0.0238289


--> Sf=sf*1000 // (final settlement due to an increase in pressure)
 Sf  = 

   23.828935


--> //results

--> printf('final settlement due to an increase in pressure=%f',Sf) //mm
final settlement due to an increase in pressure=23.828935
