![a1|
!|a2| * [b1 b2 b3]
!|a3]

subroutine col_line_prod(n, a, b, prod)
real(4) a(n), b(n), prod(n,n)

do i=1,n
prod(i,:) =a(i)*b(:)
enddo

END
