n = 16000;
M = rand(n,n);
v = rand(n);
Times = [];
sizeM = [];
count = 1;

for k = 50:50:n
  Mk = M(1:k, 1:k);
  vk = v(1:k)';
  tic;
  % Here is the main work: Matrix-vector multiplication:
  res = Mk * vk;
  t = toc;
  Times = [Times, t];
  sizeM = [sizeM, k];
  count = count + 1;
end

plot(sizeM, Times, 'o');
save Largescaling Times sizeM;