This lab focuses on a metapopulation harvesting model, which will involve two separate controls. It is based on work done by Salinas, Lenhart, and Gross in [165]. A metapopulation is a population consisting of multiple local populations. Here, we wish to model a black bear population which is motivated by the scenario in the Great Smoky Mountain National Park. We consider the bear population in a “generic” park, in the surrounding national forest, and in the “outside” or “human-populated” areas (e.g., small towns and tourist areas). These areas have shared borders, and the density-dependent emigration and immigration of the bears is based on the connectivity of the patches.