This study investigates the convective heat transfer in a magnetized flow of non-Newtonian nanofluid, incorporating multiple slip effects and the impact of nonuniform heat source/sink distributions. The model is developed by extending traditional boundary layer equations to account for complex interactions between the magnetic field, fluid flow, and thermal variations, along with activation energy fluctuations. The governing partial differential equations are transformed into ordinary differential equations (ODEs) using appropriate transformations, and the resulting system is then solved by using the shooting method with a fourth-order Runge-Kutta (RK-4) update. The effects of key parameters such as magnetic field strength (M), heat source/sink parameter (H), radiation parameter (R), and Prandtl number (Pr) on velocity, temperature, and concentration profiles are explored. The key findings include a reduction in velocity with increased values of M and lambda; an increase in temperature with higher values of M,L,Nt,epsilon,At,As, and R; and a decrease in temperature with higher values of Pr,lambda,Nb, and M. Additionally, concentration increases with rising values of M,E, and lambda, while decreasing with higher values of Le,Kr,n,So,delta, and Nb. The results are compared with existing literature and visualized by using MATLAB 2023 software. This work providing valuable insights for applications in heat exchangers, electronics cooling, and various industrial systems involving non-Newtonian nanofluids.