--> //initializations

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

   0.0015


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

   300.


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

   12.


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

   4.5


--> B=150 //m (width)
 B  = 

   150.


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

   0.3


--> //calculations

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

   0.15


--> q1=(2*%pi*k*t*d)/log((R^2)/(B*rw)) //(discharge)
 q1  = 

   0.0613618


--> q2=(2*%pi*k*t*d)/log(R/rw) //(discharge for no interference)
 q2  = 

   0.0669576


--> Pr=((q2-q1)*100)/q2 //(percentage reduction in discharge)
 Pr  = 

   8.3571622


--> printf ('Percentage reduction in discharge=%f',Pr) //%
Percentage reduction in discharge=8.357162
