circular.m

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

function Circular;
% Circular Infinity
% copyritht: zjliu
% email:[e-mail]
r=1;z=0;N=5;
close all;
hold on;
[z,r]=ccular(z,r,N);
axis square;
function [z,r]=ccular(z,r,N);
zp=z+r*exp(linspace(0,pi*2,100)*i);
plot(zp,'color',[0.9,0.8,0]);
if N>0;
    [z,r]=ccular(z-r/2,r/2,N-1);
    r=r*2;z=z+r;
    [z,r]=ccular(z,r/2,N-1);
    r=r*2;z=z-r/2;
end