c_2311.m

来自「分形与图形设计」· M 代码 · 共 23 行

r=2;
a=-0.11;b=0.65;
%a=-0.19;b=0.6557;
%a=-0.74543;b=0.11301;
%a=0;b=0.66;
%a=-0.27334;b=0.00742;
for x0=-1:0.01:1
    for y0=-1:0.01:1
        x=x0;y=y0;
        if x0^2+y0^2<1
            for n=1:80
                x1=x*x-y*y+a;
                y1=2*x*y+b;
                x=x1;y=y1;
            end
            if (x*x+y*y)<r
                plot(x0,y0);
            end
            hold on;
        end
    end
end