A steady laminar flow over a vertical stretching sheet with the existence of viscous dissipation, heat source/sink, and magnetic fields has been numerically inspected through a shooting scheme based Runge—Kutta–Fehlberg-integration algorithm. The governing equation and boundary layer balance are expressed and then converted into a nonlinear normal system of differential equations using suitable transformations. The impact of the physical parameters on the dimensionless velocity, temperature, the local Nusselt, and skin friction coefficient are described. Results show good agreement with recent researches. Findings reveal that the Nusselt number at the sheet surface augments, since the Hartmann number, stretching velocity ratio A, and Hartmann number Ha increase. Nevertheless, it reduces with respect to the heat generation/absorption coefficient δ.