c_232.m

分形与图形设计

clear all;
clc;
A=1;B=0;C=1;
for a=-1:0.005:1
    for b=-1:0.005:1
        x0=a+b*i;
        y=newton(x0);
        if (abs(y-A)<1.0e-6)
            plot(a,b,'r');hold on;
        elseif (abs(y-B)<1.0e-6)
            plot(a,b,'g');hold on;
        elseif (abs(y-C)<1.0e-6)
            plot(a,b,'b');hold on;
        end
    end
end