The shape of Bose-Einstein (or HBT) correlation functions is determined for the case when particles are emitted from a stable source, obtained after convolutions of large number of elementary random processes. The two-particle correlation function is shown to have a stretched exponential shape, characterized by the Lévy index of stability 0 < α ≤ 2 and the scale parameter R. The normal, Gaussia...