function mul(a,b) result(c)
  implicit none
  integer, intent(in) :: a, b
  integer             :: c

  c = a * b
end function mul