The non-linear equations governing the convective flow of heat, beyond marginal stability, in a fluid sphere containing a uniform distribution of heat sources, and cooled at its surface, are expressed in the Boussinesq approximation as a perturbation of the steady-state conduction solution. A steady-state axisymmetric poloidal flow is assumed, and the equations for an approximate solution, involving the first spherical harmonic only, are obtained using a special case of the general evolution criterion of Glansdorff and Prigogine. The solution is found numerically using an iterative procedure, the results showing good agreement with those obtained by Stuart's ‘shape assumption’ for values of the Rayleigh number just above the critical value for the onset of marginal stability.