The homotopy analysis method is employed to examine free convective heat and mass transfer in a steady two-dimensional magnetohydrodynamic fluid flow over a stretching vertical surface in porous medium. In this study thermal radiation and non-uniform magnetic field are taken into consideration. The two-dimensional boundary-layer governing partial differential equations are derived with considering Boussinesq and boundary-layer approximations, and the ordinary differential nonlinear forms of momentum, energy and concentration equations, obtained by the similarity solution, are solved analytically in the presence of buoyancy forces. The effects of different involved parameters such as magnetic field parameter, suction parameter, Prandtl number, buoyancy parameter, Schmidt number, Biot number and radiation parameter on velocity, temperature and concentration profiles are plotted and discussed in the paper.