cantor.m

来自「matlab的分形程序：康托集」· M 代码 · 共 16 行

function cantor(ax,ay,bx,by)
c=0.2;d=2;
if (bx-ax)>c
    x=[ax,bx];y=[ay,by];hold on;
    plot(x,y,'LineWidth',5);hold off;
    cx=ax+(bx-ax)/3;
    cy=ay-d;
    dx=bx-(bx-ax)/3;
    dy=by-d;
    ay=ay-d;
    by=by-d;
    cantor(ax,ay,cx,cy);
    cantor(dx,dy,bx,by);
end