  We investigate the four parameter family of bilateral Gamma distributions.
The goal of this paper is to provide a thorough treatment of the shapes of
their densities, which is of importance for assessing their fitting properties
to sets of real data. This includes appropriate representations of the
densities, analyzing their smoothness, unimodality and asymptotic behaviour.
