Lyapunov functions for two-species cooperative systems

Lyapunov functions for two-species mutualism models subject to proportional harvesting are introduced. Global stability of the unique coexisting equilibrium is thereby established.