Numerical analysis of the free convection coupled heat and mass transfer is presented for non-Newtonian power-law fluids with the yield stress flowing over a two-dimensional or axisymmetric body of an arbitrary shape in a fluid-saturated porous medium. The governing boundary layer equations and boundary conditions are cast into a dimensionless form by the similarity transformation. The resulting system of equations is solved by a finite difference method. The parameters studied are the rheological constants, the buoyancy ratio, and the Lewis number. Representative velocity, temperature, and concentration profiles are presented and discussed. It is found that the results depend strongly on the values of the yield stress parameter and the power-law index of the non-Newtonian fluid.