# parameters
a <--8/3;b<--10;c<-28
# Initial conditions
yini <- c(X = 1, Y = 1, Z = 1)
## Lorenz function:
Lorenz <- function (t, y, parms) {
with(as.list(y), {
dX <- a * X + Y * Z
dY <- b * (Y - Z)
dZ <- -X * Y + c * Y - Z
list(c(dX, dY, dZ)) })
}













