--> //initializations

--> k=1*10^-4 //m/sec (coefficient of permeability)
 k  = 

   0.0001


--> rw=0.1 //m (radius of well)
 rw  = 

   0.1


--> d=3 //m (depression head)
 d  = 

   3.


--> R=100 //m (radius of influence)
 R  = 

   100.


--> t=10 //m (thickness of confined aquifer)
 t  = 

   10.


--> //calculations

--> q1=2*%pi*k*rw*d //(discharge from a well with spherical flow)
 q1  = 

   0.0001885


--> q2=(2*%pi*k*t*d)/log(R/rw)//(discharge if the well had penetrated fully)
 q2  = 

   0.0027288


--> r=q1/q2 //(ratio of discharges)
 r  = 

   0.0690776


--> p=r*100  //(percentage)
 p  = 

   6.9077553


--> //results

--> printf('discharge from a well with spherical flow=%f',q1) //m^3/sec
discharge from a well with spherical flow=0.000188
--> printf('discharge if the well had penetrated fully=%f',q2) //m^2/sec
discharge if the well had penetrated fully=0.002729
--> printf('ratio of the discharge in spherical flow to in a fully penetrating well=%f',r)
ratio of the discharge in spherical flow to in a fully penetrating well=0.069078
--> printf('discharge in spherical flow is about %f of that in radial flow',p) //%
discharge in spherical flow is about 6.907755 of that in radial flow
