The present paper considers a numerical study of the non-Darcy natural convection over a vertical flat plate in a fluid-saturated porous medium. Forchheimer extension is considered in the flow equations. The nondimensional governing equations are solved by the finite element method. The resulting nonlinear integral equations are linearized and solved using the Newton− Raphson iteration. The finite element implementations are prepared by using the MATLAB® software packages. Numerical results for the details of the stream function, velocity, and temperature contours and profiles, as well as heat transfer rate in terms of Nusselt number, which are shown on graphs, have been presented.