Nothing
sir <- function() {
sigma <- 10.0
R <- 28.0
b <- 8.0 / 3.0
initial <- function(t = 0, pars = NULL) {
c(10, 1, 1)
}
derivs <- function(t, y, .) {
y1 <- y[[1L]]
y2 <- y[[2L]]
y3 <- y[[3L]]
list(c(sigma * (y2 - y1),
R * y1 - y2 - y1 * y3,
-b * y3 + y1 * y2))
}
list(derivs = derivs, initial = initial, t = c(0, 15))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.