“…Among them, the Adams-type predictor-corrector method (El-Sayed et al 2007), the iterative method (Bhalekar and Daftardar-Gejji 2012), the optimal homotopy analysis method (Mohamed 2014), the weighted Mittag-Leffler functions method (West 2015), and the spectral iterative method (Shoja et al 2016). Khan et al (2015) considered the above model and used the operational matrices of Bernstein polynomials (OMBP) for getting its approximate solutions on t ∈ [0, 1] at three different choices of r , λ, and the absolute errors were about 10 −5 , see Figs. 2, 4, and 6 in Khan et al (2015).…”