NB. ************************************* NB. Lorenz partial difference equation NB. and 3D graphics NB. original script is written Kiyoshi Yamashita NB. modified & graphics by Toshio Nishikawa NB. partial modified by M. Shimura NB. ************************************** NB. Usage: 5000 run_frac '' NB. ------------------------------------- require 'plot numeric trig' dt=: 0.005 NB. 0.005 s=: 10 NB. 10 r=: 50 NB. 50 b=: 8r3 NB. 8r3 init=: 5 8 10 lz=: 3 : 0 'X0 Y0 Z0'=: y XX=. X0 + dt*s*(Y0-X0) YY=. Y0 + dt*((r*X0)-(Y0+X0*Z0)) ZZ=. Z0+dt*((X0*Y0)-b*Z0) XX,YY,ZZ ) run_frac=: 4 : 0 NB. Usage: 5000 run_frac '' NB. x repeat times e.g. 5000 DAT=. <"1 |: lz ^:(i.x ) init pd 'reset' pd 'noaxes' pd DAT pd 'show' )