sw/src/f4/PRIM.FT

   1 C HACHTI'S PRIMZAHLEN PROGRAMM, 14. FEBRUAR 2005
   2 C
   3       INTEGER START, END, MODULO, NOPRIM
   4       INTEGER BUFFER
   5       DIMENSION XVA(10)
   6       DIMENSION YVA(10)
   7       DIMENSION XPOS(1)
   8       DIMENSION YPOS(1)
   9       DIMENSION BUFFER(2000)
  10 C
  11 C
  12 C      DATA XVA/0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0/
  13 C      DATA YVA/0.5,0.5,0.5,1.0,1.0,1.0,0.3,0.3,0.3,0.3/
  14 C
  15        CALL CLRPLT(2000,BUFFER)
  16 C
  17  1000  WRITE (4, 99)
  18  0099  FORMAT (42H PHILIPP'S PRIMZAHLENPROGRAMM - 14.2.2005   )
  19        WRITE (4,199)
  20        READ (4,399) START
  21        WRITE (4,299)
  22        READ (4,499) END
  23        WRITE (4,599)
  24  0199  FORMAT (11H START(I4):)
  25  0299  FORMAT (11H  ENDE(I4):)
  26  0399  FORMAT (I4)
  27  0499  FORMAT (I4)
  28  0599  FORMAT (H0)
  29 C
  30       XSTART=FLOAT(START)
  31       XEND=FLOAT(END)
  32       CALL SCALE(XSTART,0.0,XEND,1.0)
  33       XPOS(1)=XEND
  34       XPOS(2)=XSTART
  35       YPOS(1)=0.4
  36       YPOS(2)=0.4
  37       CALL PLOT(2,XPOS,YPOS)
  38       YPOS(1)=0.5
  39 C
  40 0010  DO 100 NOUT = START , END
  41 C
  42       XPOS(1)=FLOAT(NOUT)
  43       CALL PLOT(1,XPOS,YPOS)
  44 C
  45       NOPRIM=0
  46  20   DO 200 NIN  = 2 , NOUT
  47 C
  48       IF (NIN.EQ.NOUT) GO TO 200
  49       I=(NOUT/NIN)
  50       MODULO=NOUT-I*NIN
  51       IF (MODULO.NE.0) GO TO 200
  52       NIN=NOUT
  53       NOPRIM=1
  54 0200  CONTINUE
  55       IF (NOPRIM.EQ.0) WRITE (4,22) NOUT
  56  0022 FORMAT (I5)
  57  0100 CONTINUE
  58       PAUSE 7
  59 C
  60        END
  61 $0