Holiday greeting ! enjoy with beautiful 3D atracta
This is Lorenz like Rossler atracta M.Shimura [email protected] NB. -----------cut from fere----------------------- NB. Rossler Attractor - standalone demo, j602 NB. Original 3D script is written by Oleg Kobchenko NB. 19/12/08 M.Shimura NB. lorenz equation script is written by K.Yamashita require 'opengl' NB. coclass 'plorenz' coinsert 'jzopenglutil' gsetdefaults'' OPENGL=: 0 : 0 pc opengl; xywh 0 0 400 225;cc g isigraph opengl rightmove bottommove; pas 0 0;pcenter; rem form end; ) opengl_run=: 3 : 0 wd OPENGL ogl=: '' conew 'jzopengl' wd'pshow;' ) opengl_close=: 3 : 0 destroy__ogl'' ogl=: '' wd 'pclose' ) NB. ========================================================= NB. '`X Y Z'=: (0&{)`(1&{)`(2&{) NB. ' s r b'=: 10 , 28 , 8%3 NB. dx=: s*(Y - X) NB. dy=: (X * r - Z) - Y NB. dz=: (X*Y) - b*Z NB. dt=: 0.01 NB. I=: + dt * dx,dy,dz NB. --------------------------- dt=: 0.01 lz=:4 : 0 NB. Lorenz attracta by K.Yamashita 's r b'=. x NB. parameter 'xx yy zz'=: y NB. initial 5 8 10 or 1 1 1 X=: xx + dt*s*(yy-xx) Y=: yy + dt*((r*xx)-(yy+xx*zz)) Z=: zz + dt*((xx*yy)-b*zz) X,Y,Z ) rsr=:4 : 0 NB. Rossler attracta NB. by M.Shimura 'a b c'=. x NB. 0.2 0.2 14 'xx yy zz'=: y X=: xx + dt*-(yy+zz) Y=: yy + dt* xx + a*yy Z=: zz + dt* (b+ (xx*zz)-c*zz) X,Y,Z ) NB. --------------------------------------------------------- NB. COUNT=: 4000 NB. L=: 0.07*I^:(<COUNT) 0 1 0 NB. L=: 0.07* 10 28 8r3 lz ^:(i. 10000) 1 1 1 L=: 0.075* 0.2 0.2 14 rsr ^:(i. 20000) 1 1 1 T=: -([:-:>./+<./)"1|:L SPHERESIZE=: 0.025 NB. 0.03 'SPHERE MODEL'=: 2041+i.2 NB. GS_ROTXYZ=: 355.519 359.069 321.959 GS_ROTXYZ=: 351 353 321 GS_TRNXYZ=: 0.025 0.325 0 NB. --------------------------------------------------------- paint=: 3 : 0 if. gsinit GS_LIGHT do. quad=. {.>gluNewQuadric'' gsnewlist SPHERE gluSphere quad,SPHERESIZE,10 10 gsendlist'' gluDeleteQuadric quad init makelist MODEL '' end. glCallList MODEL gsfini'' ) init=: 3 : 0 gscolor LIME NB. GRAY glPushMatrix '' glTranslate T SPHERE list_at L glPopMatrix '' ) list_at=: 4 : 0 "1 glPushMatrix '' glTranslate y glCallList x glPopMatrix '' ) NB. ========================================================= opengl_g_paint=: paint opengl_g_char=: gschar opengl_default=: gsdefault NB. ========================================================= opengl_run'' NB. end of filr ---------------------------------------------------------------------- For information about J forums see http://www.jsoftware.com/forums.htm
