“…A comprehensive parametric study is carried out and numerical results are reported in this section through graphs. The numerical computation is performed for each parameter by assigning fixed values to other related parameters as ξ = π 2 [13,28], ω = 0.2 [28], Rb = 0.1 [28],Nt = 0.1 [13,28], Nb = 0.3 [13,28], Sc = 10 [13,28], Pe = 1 [28], Sb = 1 [28], Nr = 0.5 [13,28], Pr = 6.8 [28], δ = 0.5 [4], K c = 2 [18], Bi = 5 [13,28], K m = 2, M = 0.5 [4,17], f w [18]= 0.5, Rd = 0.5 [17], θ f = 1.1 [16], Λ = 1 [17], ε = 1 [28], β = 0.5 [37], Q = 0.2 [37], α = 0.3 [39], E = 2 [18], p = 0.2 [18], Ω = 0.1 [17], λ = 5 [13,28], and Nn = 0.2 [29]. To corroborate the correctness of the present numerical method, the results for the heat transfer rate −θ ′ (0) were contrasted against those disclosed by Merkin [15], Nazar et al [1], Rashad et al [13], Rashad et al [28], and Zahar et al [16] in Table 1.…”