Numerical implementation of Daftardar-Gejji and Jafari method to the quadratic Riccati equation

The solution of quadratic Riccati differential equations can be found by classical numerical methods like Runge-Kutta method and the forward Euler method. Batiha et al.  [7] applied variational iteration method (VIM) for the solution of General Riccati Equation. In the paper of El-Tawil et al. [19] they used the Adomian decomposition method (ADM) to solve the nonlinear Riccati equation. In [3] Abbasbandy applied Iterated He's homotopy perturbation method for solving quadratic Riccati differential equation. In [2] Abbasbandy used the Homotopy perturbation method to get an analytic solution of the quadratic Riccati differential equation, and a comparison with Adomian's decomposition method was presented. In [1] Abbasbandy employed VIM to find the solution of the quadratic Riccati equation by using Adomian's polynomials. Tan and Abbasbandy [30] employed the Homotopy Analysis Method (HAM) to find the solution of the quadratic Riccati equation. Batiha [5] used the multistage variational iteration method (MVIM) to solve the quadratic Riccati differential equation.