The production of heavy flavours at HERA has been investigated with a new Monte Carlo generator, which describes the heavy quark production in e-p scattering via boson-gluon fusion. The program has been developed in the framework of the HERWIG Monte Carlo and, different from the existing generator AROMA from the Lund package, takes account of the initial state radiation associated with higher order QCD corrections. The bbBAR production at HERA energies has been studied in the two different approaches. As expected, large differences for the particle and energy flows have been found.