利用fortran语言验证“抛硬币”的概率

使用ivf编译器
Program main
    implicit none 
    integer(kind=8) :: nn = 1000000000  !.. 设置抛硬币的总次数
    integer(kind=8) :: m, n, i 
    real(kind=8)    :: a   !.. a > 0.5: 设为硬币朝上；反之朝下
    
    m = 0  !.. 朝上的次数
    n = 0  !.. 朝下的次数
    call random_seed
    do i = 1, nn 
        call random_number( a )
        if ( a > 5d-1 ) m = m + 1
        if ( a < 5d-1 ) n = n + 1
    end do
    print*, dble(m) / nn  !.. 朝上的概率
    print*, dble(n) / nn  !.. 朝下的概率
    
End program main

由得出的结果可以看到，随着实验次数nn的增大，朝上与朝下的概率几乎趋于1/2.