m_set_by_tan3.m

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

function M_set_by_tan3; 
% [分形地图]引入正切的M集 
% Author's email:[e-mail] 
close all 
t=linspace(-1.5,1.5,200); 
[x,y]=meshgrid(t); 
z=x+y*i; 
% a=2*rand(size(z)); 
a=3; 
c=(0.55-0.42i)*ones(size(z)); 
N=zeros(size(z)); 
M=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