# To run the simulation with default parameters:
result <- simulate_Characteristics_of_ID_ode()
# To choose values other than the standard one, specify them like this:
result <- simulate_Characteristics_of_ID_ode(S = 2000,P = 2,A = 0,I = 0,R = 0,D = 0)
# You can display or further process the result, like this:
plot(result$ts[,'time'],result$ts[,'S'],xlab='Time',ylab='Numbers',type='l')
print(paste('Max number of S: ',max(result$ts[,'S'])))
Run the code above in your browser using DataLab