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