m_set_by_tan2.m

来自「分形树生成程序」· M 代码 · 共 21 行

function M_set_by_tan2; 
%分形魔鬼 引入正切的M集 
% Author's email:[e-mail] 
t=linspace(-1.5,1.5,200); 
[x,y]=meshgrid(t); 
z=x+y*i; 
% a=2*rand(size(z)); 
a=1.2; 
c=(0.3-0.709i)*ones(size(z)); 
N=zeros(size(z)); 
M=1.2; 
for k=1:200; 
   z=tan(a.*real(z).^2-imag(z).^2+real(c)).^2+... 
       i*(2*imag(z.^2)+imag(c)); 
   c(abs(z)>M)=0; 
   N(abs(z)>M)=k; 
   z(abs(z)>M)=0; 
end 
pcolor(t,t,N'); 
colormap(prism(256)) 
shading interp