Mathematical models for polygamous mating systems