! normality_of_pi.f90 ! Take N digits of pi and check for fraction of digits that are '3' character*100 line, dig*1 open(8,file='digits_pi2.txt') 10 write(6,'("Enter digit to search for")') read(5,*)n if((n.lt.0).or.(n.gt.9))then write(6,'("Digit must be in 0 to 9")') goto 10 endif ! Convert from number to character dig=char(n+48) write(6,'(a1)')dig icount=0 itot=0 do i=1,1000 read(8,'(a100)',err=100,end=100)line ! In a line, find its length ilen=len_trim(line) itot=itot+ilen do j=1,ilen if(line(j:j).eq.dig)icount=icount+1 enddo enddo 100 close(8) write(6,'(i6," digits out of ",i7," were 3")')icount,itot-1 stop end