In this research paper, we suggest a novel analytical method to obtain the approximate solutions of linear and nonlinear differential equations of arbitrary order. The proposed technique is a good combination of Kharrat–Toma transform and q-homotopy analysis method. Additionally, we also find out the solution of Cauchy reaction–diffusion equation associated with regularized version of Hilfer–Prabhakar fractional derivative. The Cauchy reaction–diffusion equation is broadly used to describe the dynamical processes in physics, geology, biology, ecology, etc. Some characteristic properties of the considered model and existence as well as the uniqueness of the solutions are also discussed. The outcomes of the considered model are exhibited graphically to show the accuracy and efficiency of the suggested method.