Abstract: A boundary fitted coordinate Finite Analytical Method (FAM) is described for the numerical solution of steady and unsteady seepage flow with pressure in foundation of sluice and dams. It has many advantages, such as adaptability to boundary, scantling of input data and small memory capacity required. The calculation shows that the results are reliable.