*4.1. Impact of Variation on θ(ζ)*

The MATLAB software is utilize to interpret the effects of variation *β*1, *β*2, *α*, Pr, *γ* and *Nt* on temperature profile. Figure 10 shows the solution plots and AE analysis plots for the temperature profile by varying of Prandtl fluid number, flexible number, ratio parameter,

Prandtl number, Biot number and thermophoresis number. The absolute error values for the following variations are 10−<sup>8</sup> to 10<sup>−</sup>3, 10−<sup>8</sup> to 10<sup>−</sup>3, 10−<sup>7</sup> to 10<sup>−</sup>3, 10−<sup>7</sup> to 10<sup>−</sup>3, 10−<sup>7</sup> to 10−<sup>4</sup> and 10−<sup>8</sup> to 10<sup>−</sup>3, respectively.

One may witness that:

