class BranchingProcess:
    def __init__(self):
        pass

    @staticmethod
    def branching_process(N, poisson_rate, t):
        """Branching_process display a phase transition from asymptotic extinction with probability 1 
        to asymptotic survival with a positive probability, depending on the average number of created particles.
        At the critical point, the avalanche size of a branching process is power-law distributed.
        
        https://www.nature.com/articles/s41598-020-69705-5

        Args:
            poisson_rate ([type]): [description]
            t ([type]): [description]
        """

        return NotImplementedError
