X-Git-Url: http://gitweb.hachti.de/?a=blobdiff_plain;f=sw%2Fsrc%2Fpascal%2FPRIMZA.PS;fp=sw%2Fsrc%2Fpascal%2FPRIMZA.PS;h=424bfbd34cec830de02523009c76684fd04c25f2;hb=81e70d488b71bf995c459ca3a02c025993460ffa;hp=0000000000000000000000000000000000000000;hpb=07ec0278333ed187ac242dedcff13c56cf1b0b91;p=pdp8.git

diff --git a/sw/src/pascal/PRIMZA.PS b/sw/src/pascal/PRIMZA.PS
new file mode 100644
index 0000000..424bfbd
--- /dev/null
+++ b/sw/src/pascal/PRIMZA.PS
@@ -0,0 +1,30 @@
+PROGRAM PRIMZAHLEN(INPUT,OUTPUT);
+
+   CONST NMAX=1000;  WNMAX=32;
+
+   VAR   X,SQUARE,
+         I,K,LIM,N:  INTEGER;
+         PRIM: BOOLEAN;
+         P: ARRAY[1..NMAX] OF INTEGER;
+         V: ARRAY[1..WNMAX] OF INTEGER;
+
+ BEGIN
+   WRITE("DIE ERSTEN  N  PRIMZAHLEN:   N="); READ(N);
+   WRITELN;
+   P[1] := 2; WRITE(2 :5); X := 1; LIM := 1; SQUARE := 4;
+   FOR I:=2 TO N DO
+   BEGIN
+     REPEAT X := X + 2;
+        IF SQUARE <= X THEN
+            BEGIN V[LIM] := SQUARE;
+              LIM := LIM + 1;SQUARE := SQR(P[LIM])
+            END ;
+        K := 2; PRIM := TRUE;
+        WHILE PRIM AND (K<LIM) DO
+          BEGIN IF V[K]<X THEN V[K] := V[K] + P[K];
+             PRIM := (X<>V[K]); K := K+1
+          END
+     UNTIL PRIM ;
+     P[I] :=X; WRITE(X :7); IF I MOD 10 = 0 THEN WRITELN
+    END
+  END.