This study is therefore aimed at developing classes of efficient numerical integration schemes, for direct solution of second-order Partial Differential Equations (PDEs) with the aid of the method of lines. The power series polynomials were used as basis functions for trial solutions in the derivation of the proposed schemes via collocation and interpolation techniques at some appropriately chosen grid and off-grid points the derivedschemes are consistent, zero-stable and convergent. the proposed methods perform better in terms of accuracy than some existing methods in the literature.