<p>The current work examines magnetohydrodynamic boundary layer flow and heat transfer across a stretching sheet enmeshed in a radiative porous medium while accounting for velocity, thermal slips and suction or injection effects. The governing nonlinear partial differential equations are transformed into a system of coupled nonlinear ordinary differential equations through suitable similarity transformations and solved numerically using the MATLAB bvp4c solver. The effects of key physical parameters, including the magnetic parameter, permeability parameter, radiation parameter, Prandtl number, pressure gradient, velocity slip, thermal slip, and suction or injection parameter, on the velocity and temperature fields are examined in detail. The results reveal that increasing the magnetic parameter, velocity slip parameter, and pressure gradient suppresses the fluid velocity, whereas higher permeability enhances the flow. The temperature field decreases with increasing thermal slip parameter, Prandtl number, and suction parameter, while it increases with stronger magnetic field, thermal radiation, pressure gradient, and permeability effects. Furthermore, the skin friction coefficient is found to increase with the magnetic and suction parameters but decreases with increasing permeability and velocity slip. The local Nusselt number increases with the Prandtl number, indicating enhanced heat transfer. The numerical results for the skin friction coefficient and heat transfer rate show excellent agreement with previously published studies, confirming the accuracy and reliability of the present formulation. The study provides useful insights into the control of momentum and thermal boundary layers in engineering systems involving electrically conducting fluids, porous structures, and thermal radiation.</p>