Three-dimensional magnetohydromagnetic slip flow of chemically reacting fluid over a variable thickness stretching sheet with space and temperature dependent heat source/sink was analysed numerically. Runge-Kutta and Newton’s methods are employed for solving the reduced ordinary differential equations with the help of similarity variables. Plots are demonstrated and examined for several parameters of concern. Also the effect of the same parameters on skin friction coefficient, heat and mass transfer rates are presented in tabular form. We found a good agreement of the present results by comparing with the published results. It is observed that space and temperature dependent heat source/sink parameters acts like controlling parameters of heat transfer. Slip effects reduce the development of concentration and thermal boundary layers.