library(deSolve)
## Chaos in the atmosphere

Q=c(1,2,3,4,5,6,7,8,9,10)
print(Q[4])
Lorenz <- function(t, state, parameters) {
with(as.list(c(state, parameters)), {
#dX <- a * X + Y * Z*t
#dY <- b * (Y - Z)
r=
dY = 10*c[t]
#dZ <- -X * Y + c * Y - Z
list(c(dY))
})
}
parameters <- c(a = -8/3, b = -10)
state <- c(Y = 0)
times <- seq(0, 5, by = 1)
out <- ode(y = state, times = times, func = Lorenz, parms = parameters)
print(out)