Numerical solution of a bifurcation problem for the Boussinesq equations at low Prandtl numbers by a multigrid method