clear all

K = [1,2,3];

for k=K
k
N = 10^k

A = rand(N,N);
x = rand(N,1);

tic
y1 = forfor(A,x);
t1 = toc

tic
y2 = forBLAS(A,x);
t2 = toc

diff = norm(y1-y2,inf)

ratio = t1/t2
end