Assumptions

• 1. The number of cockroach has reached the highest value in stable stage
• 2. Ignore natural birth and death rates in our system
• 3. Infectious individuals can not recover
• 4. Other factors that may affect the experiment are ignored

Natural condition

In natural condition indoors, due to environmental resistance like food, water, space, the population of cockroaches is more likely to follow a S-shaped growth curve (sigmoid growth curve), which can be formalized mathematically by logistic function.

With infection

Our model was constructed based on SIR epidemic model (Susceptible, Infectious, Recovered) , following are some basic properties:

• 1. Naturally all cockroaches are susceptible individuals, they can infect by M.anisopliae becoming infectious individuals.
• 2. The number of individual being infected in a contact between a susceptible and an infectious subject is simulate by standard incidence .
• 3. The transition rate between Infectious and dead is ��, its reciprocal (1/��) determines the average infectious period, which is estimate by experiment data.

Parameters

Model was simulate during 30 days, with total number of 60.

The system without so-called vital dynamics (birth and death) described above can be expressed by the following set of ordinary differential equations:

This system is non-linear, and the analytic solution does not exist, but we can compute the numerical solution by MATLAB. (see results)

Results

The following curves show dynamics number change of each kinds of individuals. We see that the infectious individuals grow fast before first 6 day, and then began to drop. The total number of cockroaches continuously going down. We specify the median lethal time (LT50), which in this condition is 11.1 days.

Sensitivity Analysis

We use sensitivity analysis to analyze the impacts of some important parameter values (��, ��) on our model outcomes (LT50). The figures below show the tendency of dead number with respect to each parameter change.