  This article is devoted to exact solutions of the Boussinesq equation that
models nonlinear shallow water waves. For this we use the Hirota bilinear
method and differential constrains. Out solutions describe in particular the
motion of the wave packets, waves on soliton and "dancing" waves. We present a
simple method for multiplication of solutions of the Hirota equation which
allow us to generate more complex structures from the waves.
