The hydrodynamic stability equations are formulated for buoyancy-induced flows adjacent to a vertical, planar, isothermal surface in cold pure water. The resulting stability equations, when reduced to ordinary differential equation by a similarity transformation, constitute a two-point boundary-value(eigenvalue) problem, which was numerically solved for various values of the density extremum parameter R=( $T_{m}$ - $T_.inf./) / ( $T_{o}$ - $T_.inf./). These stability equations have been solved using a computer code designed to accurately solve two-point boundary-value problems. The present numerical study includes neutral stability results for the region of the flows corresponding to 0.0.leq. R. leq.0.15, where the outside buoyancy force reversals arise. The results show that a small amount of outside buoyancy force reversal causes the critical Grashof number $G^*/ to increase significantly. A further increase of the outside buoyancy force reversal causes the critical Grashof number to decrease. But the dimensionless frequency parameter $B^*/ at $G^*/ is systematically decreased. When the stability results of the present work are compared to the experimental data, the numerical results agree in a qualitative way with the experimental data.erimental data.