fortran：判断一个数是否为素数

下面代码以2-30为数字区间，判断这个区间的哪些数字是素数
logical :: log = .false.
do n = 2, 30
        do i = 2, n-1
                if ( mod(n,i) == 0 ) then
                        log = .true.
                        print*, n, '=', i, '*', n/i
                        exit
                else
                        continue
                end if
                log = .false.

        end do

        if (.not.log) then
                print*, n, 'is a prime number'
        end if
        
end do