function cantr=cantor(l,a,x)
  if nargin <1
    tic
    l=0;
    N=1E7;
    cantr=ones(1,N);
  else
    cantr=a;
    l=l+1;
  endif
  if l==14
      return
  endif
  n=length(cantr);
  cantr(ceil(n/3)+1:2*ceil(n/3))=0;
  cantr(1:ceil(n/3))=cantor(l,cantr(1:ceil(n/3)));
  cantr(2*ceil(n/3)+1:end)=cantor(l,cantr(2*ceil(n/3)+1:end));
  if 1==0
    subplot (2,1,1)
    plot(cantr);
    imcan=zeros(1,length(cantr));
    for i=2:length(cantr)
      if cantr(i)
        imcan(i)=imcan(i-1)+dy;
      else
        imcan(i)=imcan(i-1);
      end
    end
    subplot (2,1,2)
    plot (imcan)
    toc
 endif
 endfunction