--> //initializations

--> T=10 //m (thickness of silty clay layer)
 T  = 

   10.


--> D=1864.64 //kg/m^3 (bulk density)
 D  = 

   1864.64


--> h=12.5 //m (gravel is under artesian pressure)
 h  = 

   12.5


--> d=2 //m (excavate a foundation trench)
 d  = 

   2.


--> p=98.1 //KN/m^2 (pressure)
 p  = 

   98.1


--> g=9.81
 g  = 

   9.81


--> //calculations

--> F=(D*g*(10^-3)*(T-d))/(g*h) //(factor of safety)
 F  = 

   1.1933696


--> //After Constructions

--> f=((D*g*(10^-3)*(T-d))+p)/(g*h) //(factor of safety)
 f  = 

   1.9933696


--> //results

--> printf('factor of safety=%f',F)
factor of safety=1.193370
--> printf('factor of safety due to uniform pressure=%f',f)
factor of safety due to uniform pressure=1.993370
