In this research, we intend to develop a dynamical system for the magnetohydrodynamic (MHD) flow of an electrically conducting Casson nanofluid on exponentially shrinking and stretching surfaces, in the presence of a velocity and concertation slip effect, with convective boundary conditions. There are three main objectives of this article, specifically, to discuss the heat characteristics of flow, to find multiple solutions on both surfaces, and to do stability analyses. The main equations of flow are governed by the Brownian motion, the Prandtl number, and the thermophoresis parameters, the Schmid and Biot numbers. The shooting method and three-stage Lobatto IIIa formula have been employed to solve the equations. The ranges of the dual solutions are f w c 1 ≤ f w and λ c ≤ λ , while the no solution ranges are f w c 1 > f w and λ c > λ . The results reveal that the temperature of the fluid increases with the extended values of the thermophoresis parameter, the Brownian motion parameter, and the Hartmann and Biot numbers, for both solutions. The presence of dual solutions depends on the suction parameter. In order to indicate that the first solution is physically relevant and stable, a stability analysis has been performed.