Fractional generalizations of the Poisson process and branching Furry process are
considered. The link between characteristics of the processes, fractional differential
equations and Levy stable densities are discussed and used for construction of the
Monte Carlo algorithm for simulation of random waiting times in fractional processes.
Numerical calculations are performed and limit distributions of the normalized variable
Z=N/<N> are found for both processes.