--> //initializations

--> c=10 //KN/m^2 (cohession)
 c  = 

   10.


--> Ys=20 //KN/m^3 (unit weight in saturation)
 Ys  = 

   20.


--> i=10 //(inclination to horizontal)
 i  = 

   10.


--> x=25 //(angle)
 x  = 

   25.


--> Yw=10 //KN/m^3 (unit weight of water)
 Yw  = 

   10.


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

   5.


--> //calculations

--> Fs=(c+(Ys-Yw)*H*(cosd(i))^2*tanm((x*%pi)/180))/(Ys*H*cosd(i)*sinm((i*%pi)/180)) //(factor of safety)
 Fs  = 

   1.907042


--> //results

--> printf('factor of safety=%f',Fs)
factor of safety=1.907042
