In the present study we consider the problem of combined buoyancy and thermocapillary convection in an upright cube with a top free surface. All other walls are considered to be solid and impermeable. The side walls are at uniform but different temperatures, whilst all other surfaces are adiabatic. The top surface deformation and interactions with the gaseous phase above are neglected. The Navier-Stokes, continuity and energy equations are cast in the velocity-vorticity formulation. The governing equations are discretized by using finite difference approximations. The solution procedure consists of a three level Alternating Direction Implicit (ADI) scheme. In order to save on computational cost, the equations were marched in time using the false transient technique. A typical 81 × 81 × 61 uniform staggered mesh has been used in the numerical computations. All computations were performed in double precision on a work-station. Results are presented for a typical fluid with a moderate Prandtl number (i.e. Pr = 7). The effects of positive and negative Marangoni number on the three-dimensional convection at different Rayleigh numbers will be considered and discussed.