The purpose of this research is to look into the solution technique for obtaining MHD velocity and temperature profiles. In the presence of a changing heat flux, the unsteady laminar boundary layer flow and heat transfer of a viscous incompressible fluid across a stretching sheet are numerically investigated. The unsteadiness is thought to be generated by a sudden increase in the surface temperature and a time-dependent stretching velocity. The flow and heat transfer partial differential equations were numerically solved using an implicit finite difference scheme and a quasi-linearization technique. Both velocity and temperature rise with time and magnetic field, according to the findings. The computed results are compared to previous work that has been published. Variable heat flux (VHF) conditions have also been taken into account.