Advertisement
Guest User

Untitled

a guest
Jul 22nd, 2017
48
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1.  
  2.     program PrimZahlen
  3.  
  4.     implicit none
  5.    
  6.     integer     :: obereGrenze
  7.     integer     :: counter
  8.     logical     :: jawohl
  9.     integer     :: i,j, at
  10.     real        :: secnds,t0,t1
  11.     integer, allocatable, dimension(:) :: prime
  12.    
  13.  
  14.  
  15.     print *, "Wie viele Zahlen sollen es denn berechnet werden ? "
  16.     read (*,*) obereGrenze
  17.    
  18.     allocate(prime(obereGrenze))
  19.    
  20.     prime(1) = 2
  21.     at = 2
  22.     counter  =1
  23.    
  24.     t0 = secnds(0.0)
  25.    
  26.     do
  27.         jawohl = .true.
  28.         do i = 1, counter
  29.             if (modulo(at, prime(i)) == 0) then
  30.                 jawohl = .false.
  31.                 at = at + 1
  32.                 continue
  33.             end if
  34.         end do  ! i
  35.         counter = counter + 1
  36.         prime(counter) = at
  37.         !print *, at
  38.         at = at + 1
  39.         if (counter == obereGrenze) then
  40.             exit
  41.         endif
  42.     end do
  43.    
  44.    
  45.     t1 = secnds(t0)
  46.    
  47.     write(*,"(A20,1X,I8)") "Anzahl Primzahlen = ", counter
  48.     write(*,"(A,1X,I10)") "Die groesste Primzahl = ", at - 1
  49.     write (*,"(A,F12.9,A)") "Man hat ", t1, "s gebraucht"
  50. 99  format(I,A)
  51.    
  52.     print*, "Zum Schliessen beliebige Taste druecken ..."
  53.     read(*,*)
  54.    
  55.     end program PrimZahlen
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement