The current study proposes an e cient numerical method to evaluate the e ects of the variable heat ux, viscous dissipation, and slip velocity on the viscous Casson Heat Transfer (CHT), considering the unsteady stretching sheet that took into account the e ect of heat generation or absorption. In this respect, Finite Element Method (FEM) was employed to solve the resulting ODEs that described the problem. The e ects of the factors governing the HT such as unsteadiness parameter, slip velocity parameter, Casson parameter, local Eckert number, heat generation parameter, and Prandtl number were studied. In addition, the local skin friction coe cient and local Nusselt number on the stretching sheet were also computed. Finally, the obtained solutions con rmed that the proposed procedure could be an easy and e cient tool for nding a solution to such uid models.