“…Moreover, Figure B shows very good matching between the exact solution ( α = 1.0), the proposed optimized solutions, and the predictor‐corrector (PC) method published in Odibat and Momani for different values of α = 0.4, 0.6, 0.8, and 1.0 in the interval t ∈ [0, 5]. In addition, Table shows a comparison between the proposed results with the results reported in three recent papers for α = 0.4 and 0.8 in the interval t ∈ [0, 1]. Therefore, the proposed method provides similar or better accuracy results than the methods published before in previous studies …”