subroutine fmult(a,b,c,n)
implicit none
real*8 :: a(n)
real*8 :: b(n)
real*8 :: c(n)
integer :: n, i
!f2py intent(hide), depend(a) :: n=len(a) 
!f2py real*8               :: a(n) 
!f2py real*8               :: b(n) 
!f2py real*8, intent(out)  :: c(n) 
do i =1,n
  c(i) = a(i) * b(i)
enddo
end