円環状カオスの式を改変してレンダリングしました
【式】
改変円環状カオス W=(x(n),y(n)) (0≦n≦9600)
x(n+1)=y(n)+Ax(n)+K/(1+x(n)^2)+L+Me^(-y(n)^2)
y(n+1)=-Bx(n)
A=0.1405 , B=1 , K=7 , L=0 , M=0.0001 , x(0)=0 , y(0)=0
プロット範囲 4＜ｎ≦９６００
参考までに POV-Rayによる
改変円環状カオス の描画方法を載せておきます
//*****************************************************************************
#macro Modified_chaos2C(Nmax,A,B,K,L,M,P0,Vctl,Lr)
#local X=P0.x; #local Y=P0.y;
#local N=0; #while(N<=Nmax)
//
#local rC=(sin(mod(Vctl.x*Y,2*pi))+1.2)/2.2;
#local gC=(cos(mod(Vctl.y*X,2*pi))+1.2)/2.2;
#local bC=(sin(mod(Vctl.z*X*Y,2*pi))+1.2)/2.2;
//
#if(N>4) sphere { <X,Y> Lr pigment { color rgb <rC,gC,bC> } } #end
//
#local Xn=X;
#local X=Y+A*X+K/(1+X*X)+L+M*exp(-Y*Y); #local Y=-B*Xn;
#local N=N+1; #end
#end // macro end
//*****************************************************************************
#local Sc=6.5; #local Tr=<-0.2,1.1>*Sc; #local Lw=0.5/Sc;
#local Vctl=<4,2,1>/15;
#local A=0.1405; #local B=1; #local K=7; #local L=0; #local M=0.0001;
union { Modified_chaos2C(9600,A,B,K,L,M,<0,0>,Vctl,Lw)
scale Sc translate Tr }
