tylermorganwall

library(deSolve)
library(rayrender)
parameters = c(a = 0.2, b = 4, c = 8, d = 1)
state = c(X = 1, Y = 0, Z = 1)

Lorenz = function(t, state, parameters) {
  with(as.list(c(state, parameters)), {
    dX = -Y^2 - Z^2 - a*X + a*c
    dY = X * Y - b * X * Z - Y + d
    dZ = b * X * Y + X * Z - Z