IMPLICIT NONE
INTEGER::M1,M2,I,J
WRITE(*,*) '请输入两个整数'
READ(*,*) M1, M2
CALL PD()
CONTAINS
!---------
SUBROUTINE PD()
IC: DO I=M1,M2
DO J=2,INT(SQRT(1.0*I))
IF(MOD(I,J)==0) CYCLE IC
ENDDO
WRITE(*,*) I
ENDDO IC
END SUBROUTINE
END
! 程序已经调试通过,保证没有问题,请放心使用。