--> //initializations

--> q=100 //KN (load)
 q  = 

   100.


--> r=2 //m (radius of circular area)
 r  = 

   2.


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

   2.5


--> //calculations

--> nz1=q*(1-((1/(1+(r/z)^2))^(3/2))) //(vertical stress at point P)
 nz1  = 

   52.386048


--> nz2=0.2*q // (vertical stress at point Q)
 nz2  = 

   20.


--> //results

--> printf('vertical stress at point P=%f',nz1) //KN/m^2
vertical stress at point P=52.386048
--> printf('vertical stress at point Q=%f',nz2) //KN/m^2
vertical stress at point Q=20.000000
