%SLIKA  narise sliko Newtonove iteracije
%       v kompleksnem za funkcijo f(z)=z^3+1
%SLIKA(obmx,obmy,nx,ny,n,tol)
%      obmx=[levix,desnix] je obmocje v x smeri
%      obmy=[spodnjiy,zgornjiy] je obmocje v y smeri
%      nx je stevilo tock v x smeri
%      ny je stevilo tock v y smeri
%      n stevilo iteracij Newtonove metode
%      tol je tolearnca za dolocitev limitne tocke
%
function slika(obmx,obmy,nx,ny,n,tol)
  clf;
  x=linspace(obmx(1),obmx(2),nx);
  y=linspace(obmy(1),obmy(2),ny);
  I=sqrt(-1); %zavarujemo konstanto i
  A=4*ones(ny,nx);%vse barve v matriki so 4 (crne)
  y=fliplr(y); %obrnemo y, da slika ni postavljena na glavo
  z1=-1;%nicle funkcije z^3+1=0 
  z2=1/2*(1+I*sqrt(3));
  z3=z2';
  for i=1:ny
    for j=1:nx
      z=newtonC('fun1','dfun1',x(j)+I*y(i),n);
      if norm(z-z1)<=tol%namesto norm lahko uporabimo abs
        A(i,j)=1;%dolocimo barvo zacetnega priblizka
      elseif norm(z-z2)<=tol
        A(i,j)=2;
      elseif norm(z-z3)<=tol
        A(i,j)=3;
      end
    end
  end
  M=[1 0 0;0 1 0;0 0 1;0 0 0]; %matrika barv za RGB sistem
  colormap(M);
  image(A);  
