The main purpose of this manuscript is to present a new approximation method to find the numerical solution of partial differential equations with Caputo fractional derivatives. We have also included the stability analysis of the proposed technique. To illustrate the effectiveness and reliability of the method, we have examined the method with several numerical examples. The high computational efficiency of the proposed method can be seen by observing the accurate approximate solutions obtained from the process. The suggested procedure is a potent tool and it can be used to solve other similar problems in modeling real-world problems.
