integer i, j, n
  real a(n,n), x(n),  y(n)

  do i = 1, n
    y(i) = 0.0
    do j = 1, n
      y(i) = y(i) + a(i,j) * x(j)
    end do
  end do