--> //initializations

--> e1=1.0 //(void ratio at 1)
 e1  = 

   1.


--> e2=0.8 //(void ratio at 2)
 e2  = 

   0.8


--> R1=85 //KN/m^2 (pressure at 1)
 R1  = 

   85.


--> R2=465 //KN/m^2 (pressure at 2)
 R2  = 

   465.


--> R3=600 //KN/m^2 (pressure at 3)
 R3  = 

   600.


--> //calculations

--> Cc=((e1-e2)*2.303)/(log(R2/R1)) //(coefficient of compressibility)
 Cc  = 

   0.271039


--> e3=e1-((Cc*log(R3/R1))/2.303) //(void ratio at 3)
 e3  = 

   0.7700018


--> //results

--> printf('void ratio for pressure R3=%f',e3)
void ratio for pressure R3=0.770002
