Gillespie Method of SEIR on metapopulations