Penetrative convection due to purely internal heating in a horizontal ferrofluid-saturated porous layer is examined by performing linear stability analysis. Four different types of heat supply functions are considered. The Darcy model is used to incorporate the effect of the porous medium. Numerical solutions are obtained by using the Chebyshev pseudospectral method, and the results are discussed for all three boundary conditions: when both boundaries are impermeable and conducting; when both boundaries are conducting with lower boundary impermeable and free upper boundary; and when both boundaries are impermeable with lower boundary conducting and upper with constant heat flux. The effect of the Langevin parameter, width of ferrofluid layer, permeability parameter, and nonlinearity of the fluid magnetization has been observed at the onset of penetrative convection for water-and ester-based ferrofluids. It is seen that the Langevin parameter, width of ferrofluid layer, and permeability parameter have stabilizing effects on the onset of convection, while the nonlinearity of the fluid magnetization advances the onset of convection.