Outer product of two vectors in Fortran

Unfortunately, no outer product intrinsic function available in Fortran. Outer product return a matrix from two user-provided vector a and b,

where