subroutine forMatMul(A, B, n, C)

      implicit none

      integer n
      double precision A(n,n)
      double precision B(n,n)
      double precision C(n,n)

!f2py intent(out) :: C
!f2py intent(in)     :: A
!f2py intent(in)     :: B
!f2py intent(in)     :: n

      C = matmul(A,B)

      return

      end subroutine forMatMul