Lorenz attractor

LorenzAttractor.tq

let trajectory p dt n =
  let s, t, u = 10, 8/3, 28 in
  Array.unfold [(x, y, z), i ->
    if i=n then None else
    let p = x+s*(y-x)*dt, y+(x*(u-z)-y)*dt, z+(x*y-t*z)*dt in
    Some(p, (p, i+1))] (p, 0)

let () =
  trajectory (10, 0, 20) (1/200) 2000
  @ Array.map [x, y, z -> x, y]
  @ Chart.scatter "" "" ""