--> //initializations

--> H1=5*1000 //m (thickness of layer first)
 H1  = 

   5000.


--> H2=1*1000 //m (thickness of layer second)
 H2  = 

   1000.


--> H3=2.5*1000 //m (thickness of layer third)
 H3  = 

   2500.


--> kh1=3*(10^-2) //mm/sec (coefficient of permeability first layer)
 kh1  = 

   0.03


--> kh2=3*(10^-5) //mm/sec (coefficient of permeability second layer)
 kh2  = 

   0.00003


--> kh3=4*(10^-2) //mm/sec (coefficient of permeability third layer)
 kh3  = 

   0.04


--> //calculations

--> kh=(kh1*H1+kh2*H2+kh3*H3)/(H1+H2+H3) 
 kh  = 

   0.0294153


--> kv=(H1+H2+H3)/((H1/kh1)+(H2/kh2)+(H3/kh3))
 kv  = 

   0.0002533


--> //results

--> printf('coefficient of permeability in horizontal direction=%f',kh) //mm/sec
coefficient of permeability in horizontal direction=0.029415
--> printf('coefficient of permeability in vertical direction=%f',kv) //mm/sec
coefficient of permeability in vertical direction=0.000253
