The convection of a Maxwell fluid over a stretching porous surface with a heat source/sink in the presence of nanoparticles is investigated. The Lie symmetry group transformations are used to convert the boundary layer equations into coupled nonlinear ordinary differential equations. The ordinary differential equations are solved numerically by the Bvp4c with MATLAB, which is a collocation method equivalent to the fourth-order mono-implicit Runge-Kutta method. Furthermore, more attention is paid to the effects of the physical parameters, especially the parameters related to nanoparticles, on the temperature and concentration distributions with consideration of permeability and the heat source/sink.