The evidence of square patterns in B\'enard-Marangoni convection relatively far above threshold is demonstrated for the first time by means of a direct numerical integration of the three-dimensional Navier-Stokes equations. It is shown that regular squares are only obtained for a finite Prandtl number. Finally a qualitative explanation of the obtained sequence of instabilities based on a model of amplitude equations is discussed.