--> //initializations

--> Q=100 //KN (vertical load)
 Q  = 

   100.


--> r=3 //m (radial distance)
 r  = 

   3.


--> z=3 //m (depth)
 z  = 

   3.


--> //calculations

--> nz=(Q/(%pi*(z^2)))*(1/(1+2*(r/z)^2))^(3/2) // (vertical stress at point P)
 nz  = 

   0.680653


--> //results

--> printf('vertical stress at point P=%f',nz) //KN/m^2
vertical stress at point P=0.680653
