--> //initializations

--> H=15 //m (height of water table)
 H  = 

   15.


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

   12.


--> h=9 //m (height of reduce lavel)
 h  = 

   9.


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

   300.


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

   0.0005


--> D=0.30 //m (diameter of the well)
 D  = 

   0.3


--> //calculations

--> rw=D/2 //(radius of the well)
 rw  = 

   0.15


--> q=(%pi*k*(2*t*H-t^2-h^2))/log(R/rw) //(discharge)
 q  = 

   0.027899


--> //results

--> printf('discharge=%f',q) //m^3/sec
discharge=0.027899
