This article presents a new method of solving partial differential equations. The method is an improvement of the previously reported compact finite difference quasilinearization method (CFDQLM) which is a combination of compact finite difference schemes and quasilinearization techniques. Previous applications of compact finite difference (FD) schemes when solving parabolic partial differential equations has been solely on discretizing the spatial variables and another numerical technique used to discretize temporal variables. In this work we attempt, for the first time, to use the compact FD schemes in both space and time. This ensures that the rich benefits of the compact FD schemes are carried over to the time variable as well, which improves the overall accuracy of the method. The proposed method is tested on four nonlinear evolution equations. The method produced highly accurate results which are portrayed in tables and graphs.