subroutine fmult(a,b,c,n)
implicit none
real*8 :: a(*)
real*8 :: b(*)
real*8 :: c(*)
integer :: n, i
do i =1,n
  c(i) = a(i) * b(i)
enddo
end