--> //initializations

--> u1=10 //KN/m^2 (final pore pressure)
 u1  = 

   10.


--> u2=-60 //KN/m^2 (initial pore pressure)
 u2  = 

  -60.


--> r1=100 //KN/m^2 (final cell pressure)
 r1  = 

   100.


--> r2=0 //KN/m^2 (initial cell pressure)
 r2  = 

   0.


--> up=-70 //KN/m^2 (pore pressure)
 up  = 

  -70.


--> rd=500 //KN/m^2 (deviator stress)
 rd  = 

   500.


--> //calculations

--> B=(u1-u2)/(r1-r2) //(pore pressure coefficient B)
 B  = 

   0.7


--> A=(up-u1)/(B*rd) //(pore pressure coefficient A)
 A  = 

  -0.2285714


--> //results

--> printf('pore pressure coefficient A=%f',A)
pore pressure coefficient A=-0.228571
--> printf('pore pressure coefficient B=%f',B)
pore pressure coefficient B=0.700000
