%   time_factor.m  measure time needed to factor numbers
N=50;
  t=zeros(1,N);
  for k=4:N
     for j=1:30
        tic;
        factor(randi(2^k));
        t(k)=toc+t(k);
      end
  end
  plot(log(t/N))

