在Fortran中，可以使用matmul函数来执行矩阵乘法操作。为了测试不同乘法尺寸下matmul函数的时间，可以编写一个简单的程序来计算不同尺寸矩阵乘法的时间。

以下是一个示例程序，用于计算不同乘法尺寸下matmul函数的时间：

program matrix_multiplication_time
  implicit none

  integer, parameter :: N = 1000 ! 最大矩阵尺寸
  integer :: i, j, k
  real :: start_time, end_time, elapsed_time
  real, dimension(N, N) :: A, B, C

  ! 随机生成矩阵 A 和 B
  call random_seed()
  call random_number(A)
  call random_number(B)

  do i = 100, N, 100 ! 递增矩阵尺寸
    start_time = omp_get_wtime() ! 获取开始时间

    ! 执行矩阵乘法
    C = matmul(A(1:i, 1:i), B(1:i, 1:i))

    end_time = omp_get_wtime() ! 获取结束时间
    elapsed_time = end_time - start_time ! 计算经过的时间

    write(*, '(A, I5, A, F6.4, A)') "矩阵尺寸:", i, " 时间:", elapsed_time, " 秒"
  end do

end program matrix_multiplication_time

在上面的示例程序中，我们使用随机数生成了两个大小为N×N的矩阵A和B。然后，我们使用matmul函数执行矩阵乘法操作，并记录开始时间和结束时间。最后，我们计算经过的时间，并将其打印到屏幕上。

请注意，在上面的示例程序中，我们使用了omp_get_wtime()函数来获取当前时间。这个函数是OpenMP库提供的，用于计算时间。如果你的编译器不支持OpenMP，你可以使用其他方法来计算时间，比如使用系统时间函数。