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

diff --git a/sw/src/pascal/QUAGL.PS b/sw/src/pascal/QUAGL.PS
new file mode 100644
index 0000000..35ae186
--- /dev/null
+++ b/sw/src/pascal/QUAGL.PS
@@ -0,0 +1,22 @@
+PROGRAM QUAGL(INPUT,OUTPUT);
+
+  VAR   A,B,C,DISKR,X1,X2,RE,IM: REAL;
+
+BEGIN
+    READ(A,B,C);
+    DISKR := B*B - 4*A*C;
+    IF DISKR>=0 THEN
+                    BEGIN
+                        X1 := (-B + SQRT(DISKR))/(2*A);
+                        X2 := (-B - SQRT(DISKR))/(2*A);
+                        WRITELN("X1=", X1:10:5);
+                        WRITELN("X2=", X2:10:5)
+                    END
+                ELSE
+                    BEGIN
+                        RE := -B/(2*A);
+                        IM := SQRT(-DISKR)/(2*A);
+                        WRITELN("X1=", RE:10:5, "     +", IM:10:5, " *J");
+                        WRITELN("X2=", RE:10:5, "     -", IM:10:5, " *J")
+                    END
+END.