clear;

%for k=1:11
%  n=2^k;
  %n=1000;
nlist=50:50:2000;
for k=1:length(nlist)
  n=nlist(k);
  A=rand(n);
  B=rand(n);
  C=zeros(n);
  D=zeros(n);
  u=rand(n,1);
  v=zeros(n,1);

  tic;
  v=A*u;
  t1(k)=toc;

  tic;
  C=A*B;
  t2(k)=toc;

  tic;
  for j=1:n
    D(:,j)=A*u;
  end 
  t3(k)=toc;

  disp('------------------------------------');
  disp(sprintf('Size: n=%f',n));
  disp(sprintf('   Mat-vec elapsed time: %8.4f',t1(k)));
  disp(sprintf('   Mat-mat elapsed time: %8.4f',t2(k)));
  disp(sprintf('n Mat-vecs elapsed time: %8.4f',t3(k)));
  disp(sprintf(' n*Mat-vec elapsed time: %8.4f',n*t1(k)));
  disp('------------------------------------');
  N(k)=n;
end

loglog(N,t1,N,t2,N,t3);
legend('m-v','m-m','n-m-v');
grid on;