Verificare dacă n este număr prim

Dacă se consideră un număr N se pune problema să se vadă dacă este sau nu număr prim. Acest lucru se realizează prin împărțiri succesive la numere impare mai mici decât radicalul numărului N. Dar ca să nu fie efort făcut fără rost se verifică mai întâi dacă N este număr par. Subprogramul returnează IK=0 dacă numărul N nu este număr prim și IK=1 dacă numărul N este număr prim.

         SUBROUTINE COMPAR (N,IK)
         IK=0

         IF(MOD(N,2).EQ.0) GO TO 20

         A=N

         M=SQRT(A)
         DO 10 I=2,M

         IF(MOD(N,I).EQ.0) GO TO 20

10     CONTINUE

         IK=1

20     CONTINUE

         RETURN
         END

(26 decembrie 2017)