

function s = pwell_ex(pw,qw,rw,K,r)
  
  pi = 4.*atan(1);
  s = pw + qw/(2*pi*K)*log(r/rw);

endfunction  



function s = qwell_ex(qw,r)
  
  pi = 4.*atan(1);
  s = qw/(2*pi*r);

endfunction  











  
  

