Generate a poisson process

 

    
main program 

    
sub routine