A family of numerical method, based on rational approximation to matrix exponential function is developed for solving parabolic partial differential equation. The method is forth order accurate in space and time. In this method second order spatial derivatives is approximated by forth order finite difference approximation. Complex arithmetics are not required in proposed method. Algorithm is developed and accuracy of method is analyzed with some examples. Finally the results are compared with existing ones, which shows that our results are more accurate