| 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 |