Abstract
The approximate analytical solution of fractional order, nonlinear, reaction differential equations, namely the nonlinear diffusion equations, with a given initial condition, is obtained by using the homotopy analysis method. As a demonstration of a good mathematical model, the present article gives graphical presentations of the effect of the reaction terms on the solution profile for various anomalous exponents of particular cases, to predict damping of the field variable. Numerical computations of the convergence control parameter, used to evaluate the convergence of approximate series solution through minimizing error, are also presented graphically for these cases.