Массив простых чисел

!Без деления задаем границу массива
!но при этом явно указываем шаг просеивания
forall (i=4 :100:2) a(i) = 0
forall (i=9 :100:3) a(i) = 0
forall (i=25:100:5) a(i) = 0
forall (i=49:100:7) a(i) = 0

program eratosfen
implicit none
integer :: a(2:100), j, n

forall (j=2:100) a(j) = j

	n = 2	
	do while (n*n < 100)
	
	   forall (j = n*n : 100 : n) a(j) = 0
	   n = n + 1
		
	      do while (a(n)==0)
	         n = n + 1
	      end do
	end do

  n = count (a/=0); n = n + 1

  a(2:n) = pack (a, a/=0)

print 100, a(2:n)

100 format (10i5)
end program eratosfen