function y = forBLAS(A,x)

[M,N] = size(A);

if N ~= length(x)
  disp('Dimensioni della matrice e del vettore incompatibili')
  return
end

y = zeros(N,1);
for i=1:M
    y(i) = A(i,[i:N])*x([i:N]);
end