In this paper, the authors derive a meshless algorithm using radial basis functions for quasilinear parabolic equations of second order. Numerical tests are performed for a one-dimensional drainage problem under evaporation, and it is demonstrated that the proposed algorithm may provide accurate numerical solutions. A specific farmland drainage problem further illustrates this application.