In this paper, the rate of heat transfer of the steady MHD stagnation point flow of Casson fluid on the shrinking/stretching surface has been investigated with the effect of thermal radiation and viscous dissipation. The governing partial differential equations are first transformed into the ordinary (similarity) differential equations. The obtained system of equations is converted from boundary value problems (BVPs) to initial value problems (IVPs) with the help of the shooting method which then solved by the RK method with help of maple software. Furthermore, the three-stage Labatto III-A method is applied to perform stability analysis with the help of a bvp4c solver in MATLAB. Current outcomes contradict numerically with published results and found inastounding agreements. The results reveal that there exist dual solutions in both shrinking and stretching surfaces. Furthermore, the temperature increases when thermal radiation, Eckert number, and magnetic number are increased. Signs of the smallest eigenvalue reveal that only the first solution is stable and can be realizable physically.