In this study, the models for a double-diffusive thermohaline system heated from below were developed and the governing equations were established taking account of the density variation with time. The six order Runge-Kutta method was used for the solution of the simultaneous governing differential equations and the temperature and salt concentration distributions and height of each layer within the system were predicted. As the result of this study, it was found that the predicted values with the convective layer growing proportionality constant of 0.18 showed a good agreement with available experimental data. It was also found that the effect of density change with time on the temperature profile in the bottom convective layer could not be negligible.