We deal with a multiple contact process, that is a modification of the contact process. This system contains N kinds of species (N = 10) on a finite-sized square lattice. Simulations are carried out using two different methods: local and global interactions. It is found that the waiting time to the extinction of the first species follows a power law. Moreover, we find that local interaction promotes the coexistence of multiple species.