---

A meshless method namely, discrete least square method (DLSM), is presented in the paper for the solution of free surface seepage problem. In this method computational domain is discredited by some nodes and then the set of simultaneous equations are built using moving least square (MLS) shape functions and least square technique. The proposed method does not need any background mesh therefore it is a truly meshless method. Several numerical two dimensional examples of Poisson partial differential equations (PDEs) are presented to illustrate the performance of the present DLSM. And finally a free surface seepage problem in a porous media is solved and results are presented.