cantor3.m

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

function Cantor3(n);
% n is the number of recursion
% 
% \copyright:zjliu
% Author's email: [e-mail]
if nargin==0;
    n=5;
end
hold on;
[x,L]=ccp(0,0.5,n,0);
axis off;
function [x,L]=ccp(x,L,n,y);
plot([x-L,x+L],[y,y],'k','linewidth',2);
if n>1;
    [x,L]=ccp(x-L/1.5,L/3,n-1,y-1);
    L=L*3;x=x+L/1.5;
    [x,L]=ccp(x+L/1.5,L/3,n-1,y-1);
    L=L*3;x=x-L/1.5;
end