function y = forfor(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
  for j=i:N
    y(i) = y(i) + A(i,j)*x(j);
  end
end