In this study, the general partial differential equation of unsteady groundwater flow through non-homogeneous and anisotropic aquifer was solved by the fully implicit finite difference method with control volume approach. The calibration of permeability and storage coefficients as model parameters was performed by a transformation technique and application of Powell's non-linear optimization me...