clear;
n=1000;
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=toc;

tic;
C=A*B;
t2=toc;

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

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