%% This is a code to calculate female's fecunduty based on the number of matings

function F=fecundity(c,matings)
F=(matings>0).*exp(-c*(matings-1));