**4. Conclusion Remarks**

The magnetohydrodynamic (MHD) flow over a shrinking sheet and heat transfer with viscous dissipation is numerically studied. The governing equations of the momentum and energy are transformed into the ordinary differential equations by using similarity transformation. The resultant equations have been transformed into a system of fractional differential equations by using the Caputo derivative. Fractional differential boundary layer equations, based on Caputo operators, are solved numerically by the Adams-type predictor-corrector method. We compared our result with a past published article and found it in good agreement with the first solution. Further, our results of the second solution did not concur with the published results. This was because of the heat equation and its different parameter effects on the unstable (second) solution. On the other hand, there existed two different ranges, namely, no solution and dual solutions, which depended on the magnetic and suction parameters.

**Author Contributions:** Z.O. and L.A.L. formulated the problem. L.A.L. and S.O.A. transformed the problem and computed the results. I.K. and K.S.N. discussed the results physically. All the authors equally contributed in writing the manuscript.

**Funding:** This research received no external funding.

**Acknowledgments:** The authors would like to thank Deanship of Scientific Research, Majmaah University for supporting this work.

**Conflicts of Interest:** The authors declare no conflict of interest.
