The weakly nonlinear stability analysis of a convective flow in a planar vertical fluid layer is performed in this paper. The base flow in the vertical direction is generated by internal heat sources distributed within the fluid. The system of Navier-Stokes equations under the Boussinesq approximation and small-Prandtl-number approximation is transformed to one equation containing a stream function. Linear stability calculations with and without a small-Prandtl-number approximation lead to the range of the Prantdl numbers for which the approximation is valid. The method of multiple scales in the neighborhood of the critical point is used to construct amplitude evolution equation for the most unstable mode. It is shown that the amplitude equation is the complex Ginzburg-Landau equation. The coefficients of the equation are expressed in terms of integrals containing the linear stability characteristics and the solutions of three boundary value problems for ordinary differential equations. The results of numerical calculations are presented. The type of bifurcation (supercritical bifurcation) predicted by weakly nonlinear calculations is in agreement with experimental data.