--> //initializations

--> qp=0.04 //m^3/sec (discharge)
 qp  = 

   0.04


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

   0.25


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

   200.


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

   0.004


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

   20.


--> W=9 //m (penetrates a confined aqifer)
 W  = 

   9.


--> //calculations

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

   0.125


--> C=(W/t)*(1+7*sqrt(rw/(2*W))*cosd((%pi*W*180)/(2*t*%pi))) //(correction factor)
 C  = 

   0.6496066


--> d=(qp*log(R/rw))/(2*%pi*k*t*C) //m (drawdown at the well)
 d  = 

   0.903783


--> //results

--> printf('drawdown at the well=%f',d) //m
drawdown at the well=0.903783
