This model allows for the simulation of a parasitic infection that requires an intermediate host for transmission

simulate_parasites_ode(
Sh = 1000,
Ih = 1,
E = 1,
Si = 0,
Ii = 0,
tmax = 120,
bi = 0.01,
be = 0.01,
m = 0,
n = 0,
g = 0,
w = 0,
p = 0.01,
c = 0.001
)

## Arguments

Sh : initial number of susceptible definitive hosts : numeric : initial number of infected definitive hosts : numeric : initial number of pathogens in the environment : numeric : initial number of susceptible intermediate hosts : numeric : initial number of infected intermediate hosts : numeric : maximum simulation time : numeric : rate of transmission from infected intermediate host to susceptible host : numeric : rate of transmission from environment to susceptible intermediate host : numeric : the rate of births of intermediate hosts : numeric : the rate of natural intermediate hosts : numeric : the rate at which infected hosts recover/die : numeric : the rate at which host immunity wanes in host : numeric : rate at which infected host shed the pathogen in the environment : numeric : rate at which the pathogen decays in the environment : numeric

## Value

This function returns the simulation result as obtained from a call to the deSolve ode solver.

## Details

A compartmental ID model with several states/compartments is simulated as a set of ordinary differential equations. The function returns the output from the odesolver as a matrix, with one column per compartment/variable. The first column is time.

## Warning

This function does not perform any error checking. So if you try to do something nonsensical (e.g. negative values or fractions > 1), the code will likely abort with an error message.

  # To run the simulation with default parameters just call the function:
plot(result$ts[ , "time"],result$ts[ , "Sh"],xlab='Time',ylab='Number Susceptible',type='l')