*3.2. Numerical Solution*

The numerical (ND solve) solution of Equations (13)–(16) with boundary conditions (17) and (18) for different values of embedded parameters are calculated and compared with HAM in Tables 4–7.

**Table 4.** Comparison of HAM and numerical solution for velocity when *h* = −0.001, *Pr* = 10.6, *Nr* = *R* = *Sc* = *Sr* = τ = Δ = 0.5, *Gr* = β = *Mr* = 1.



**Table 5.** Comparison of HAM and numerical solution for microrotation angular velocity when *h* = −0.15, *Nr* = *R* = Δ = *Sc* = *Sr* = *τ* = 0.5, *Mr* = β = 1, *Gr* = 5, *Pr* = 10.6.

**Table 6.** Comparison of HAM and numerical solutions for temperature when *h* = −0.33, *Nr* = *R* = Δ = *Sc* = *Sr* = *τ* = 0.5, *Mr* = β = 1, *Gr* = 5, *Pr* = 10.6.


**Table 7.** Comparison of HAM and numerical solutions for concentration when *h* = −0.42, *Nr* = *R* = Δ = *Sc* = *Sr* = *τ* = 0.5, *Mr* = β = *Gr* = 1, *Pr* = 10.6.

