The flow of MHD nanofluid towards a permeable stretching/shrinking sheet with chemical reaction, partial slip and heat source effects at the stagnation point has been examined. The state of this partial slip boundary condition and heat source in flow configuration has devoted a considerable attention on account of its widespread application in industrial and chemical engineering. Slip velocity is included to the boundary condition as a part of entrenched stretching/shrinking velocity. The leading partial differential equations are reduced to a set of ordinary differential equation systems with the help of similarity transformations and then it has been clarified utilizing the boundary value problem solver (bvp4c) package in MATLAB. The skin friction coefficient, local Nusselt number, and local Sherwood number, as well as velocity, temperature, and nanoparticle volume fraction profiles for certain values of the leading parameters, namely partial slip parameter and heat source parameter, are calculated numerically in the study. Results indicates that nanoparticle and temperature volume fraction rise as the heat source variable enhances while the velocity remain constant. The size of the skin friction coefficient, Nusselt number, and Sherwood number are all significantly lowered in the first solution when heat source and slip factors rise. For a given range of the stretching/shrinking variable, it is determined that there are two different solutions persist.