--> //initializations

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

   1850.


--> Yw=9.81 //(unit weight of water)
 Yw  = 

   9.81


--> g=9.81
 g  = 

   9.81


--> h=8 //m (height)
 h  = 

   8.


--> d=5 //m (depth)
 d  = 

   5.


--> //calculations

--> ha=0.42*h // (height of base of soil prism)
 ha  = 

   3.36


--> Yd=(D*g*(10^-3))-g //(dry unit weight)
 Yd  = 

   8.3385


--> F=(Yd*d)/(Yw*ha) //(factor of safety)
 F  = 

   1.264881


--> //results

--> printf('factor of safety=%f',F) 
factor of safety=1.264881
