% BASINS12.M
% Calls NEWTBASN.M to generate Newton basins 
% for the 12th degree polynomial poly12 below
% with roots [-3 -1 1 3], [-3 -1 1 3] +/- 2i

% C. Henry Edwards  June 1998

poly12 = [1  0 -14  0  183  0 -612  0 ...
             -2209  0 -35374  0  38025]; 


newtmap12 = [ 0.00  0.00  0.00   % black
              0.50  0.40  0.00   % brown
              0.00  0.65  0.25   % green
              0.80  0.00  0.10   % red
              0.00  0.87  0.95   % light blue
              0.00  0.24  0.87   % blue
              0.95  0.71  0.00   % tan

              0.51  0.00  0.51   % purple

              1.00  0.48  0.00   % orange

              0.87  0.48  0.48   % pink

              0.50  0.00  0.00   % crimson

              0.00  0.00  0.48   % dark blue

              0.87  0.87  0.32   % yellow

              0.00  1.00  0.00   % bright green

              0.48  0.48  0.48   % gray

              0.87  0.87  0.87]; % white



% UN-REMARK THE WINDOW YOU WANT TO SEE
wind12 = [-4 4 -3 3];                      % The big picture
% wind12 = [-0.48 0.48 -0.36 0.36];        % The flower
% wind12 = [ 0.09  0.21  0.225  0.315];    % The bud

count = 50;          % Decrease for rough draft
flops0 = flops;
t0 = clock;
w = newtbasn(poly12,wind12,newtmap12,count);
t1 = etime(clock,t0);
flops1 = flops - flops0;
rate = flops1/t1           
% About 6 million flops/sec with a 266 MHz Pentium II
% About 25 million flops per iteration, 
% and hence about 4 seconds per iteration

image(w), colormap(newtmap12), axis off