“…To handle partial differential equations (PDEs), having order fraction is of physical importance, and effective, trustworthy, and appropriate numerical methods are required [12][13][14]. Several major strategies have been utilized in this regard, including the fractional operational matrix method (FOMM) [15], Elzaki transform decomposition method (ETDM) [16,17], homotopy analysis method (HAM) [18], homotopy perturbation method (HPM) [19,20], iterative Laplace transform method [21], and variational iteration method (FVIM) [22].…”