fortran 途中まで考えたのですが。。。。
エラトステネスのふるい（素数の倍数を除いていって残ったのが素数）で３桁の素数を求めて表示したいです。
私が考えたのは、
１・２～９９までの数を素数かどうか調べて、素数を配列に入れていく
２・１００～９９９まで素数の配列の中の数で割って、割り切れたらおしまい。割り切れなかったら表示していく
ということです。
しかし下のプログラムではうまく素数配列ができていないようなのです。
ここまでかなり時間がかかったのでこのプログラムに手をいれて
これ以外におかしくなるところもどこを直せばいいのか教えてくださるとうれしいです。
C
C q223.f
C
PROGRAM q223
C
IMPLICIT NONE
C
INTEGER N,i,K,s,l
REAL a(99),b(99),c(99),X,Y
C
real M
C
a(1)=2
a(2)=3
l=2
C
DO N=2,99,1
M=N**(0.5)
S=M
DO i=2,S,1
K=MOD(N,i)
IF(K ==0)THEN
exit
ELSE IF(K /=0)THEN
l=l+1
a(l)=N
ENDIF
ENDDO
ENDDO
C
do N=100,999
do l=1,99
X=a(l)
Y=N/X
if(Y==0)then
exit
else if(Y/=0)then
write(*,*)N
end if
end do
end do
c
end
よろしくおねがいします
お礼
ありがとうございます。