m_set_by_tan1.m

来自「分形计算的几个程序原码」· M 代码 · 共 21 行

function M_set_by_tan1; 
% 引入正切的M集 
% Author's email:[e-mail] 
t=linspace(-1.5,1.5,200); 
[x,y]=meshgrid(t); 
c=x+y*i; 
a=rand(size(c)); 
z=zeros(size(c)); 
N=zeros(size(c)); 
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