rtMortalityStableSeek: seeks day1 mortality rates for F&M that give popn stability
In AndySouth/rtsetse: simulating tsetse fly populations

Description Usage Arguments

rtMortalityStableSeek seeks a age1 mortality rate for F that leads to population stability and for M that leads to the user supplied ratio of males to females. It does this by testing successive mortality rates, starting at very low values and increasing until target values are met. In each case age1 mortality rates are converted into age dependent mortalities using other user inputs. The resulting mortalities are used to fill an age structure from age1. For F, user supplied parameters are used to simulate larval production. The target is to produce the same number of larvae as age1 pupae - thus leading to a stable population. Once a stable female age structure has been achieved, the function searches in the same way for the age1 male mortality rate that generates a male age structure to satisfy the to the user supplied ratio of males to females.

rtMortalityStableSeek(iMaxAge = 120, iTargetPopAge0 = 100,
  fInterval = 0.001, fHighestTestMort = 0.3, fMperF = 0.5,
  pMortLarva = 0.05, iMortMinAgeStartF = 10, iMortMinAgeStopF = 60,
  fMortMinPropF = 0.2, fMortOldPropF = 0.3, iMortMinAgeStartM = 10,
  iMortMinAgeStopM = 40, fMortMinPropM = 0.2, fMortOldPropM = 0.3,
  propMortAdultDD = 0.25, pMortPupa = 0.25, iPupDurF = 26,
  iPupDurM = 28, iFirstLarva = 16, iInterLarva = 10, plot = TRUE,
  verbose = FALSE)

`iMaxAge`	max age of fly allowed in model
`iTargetPopAge0`	target pop of both sexes at age 0
`fInterval`	the interval between mortalities used to test
`fHighestTestMort`	the highest mortality tested
`fMperF`	desired numbers of males per female, default=0.5
`pMortLarva`	larval mortality per period
`iMortMinAgeStartF`	Age at which min death rates start.
`iMortMinAgeStopF`	Age at which min death rates stop.
`fMortMinPropF`	What proportion of the maximum death rate on day 0 is the minimum death rate.
`fMortOldPropF`	What proportion of the maximum death rate on day 0 is the death rate after iDeathMinAgeStop.
`iMortMinAgeStartM`	Age at which min death rates start.
`iMortMinAgeStopM`	Age at which min death rates stop.
`fMortMinPropM`	What proportion of the maximum death rate on day 0 is the minimum death rate.
`fMortOldPropM`	What proportion of the maximum death rate on day 0 is the death rate after iDeathMinAgeStop.
`propMortAdultDD`	proportion of adult mortality that is density dependent
`pMortPupa`	pupal mortality per period
`iPupDurF`	days it takes pupa(F) to develop
`iPupDurM`	days it takes pupa(M) to develop
`iFirstLarva`	Age that female produces first larva
`iInterLarva`	Inter-larval period
`plot`	whether to output plots
`verbose`	whether to output progress text #' @return float mortality probability