sw/src/pascal/QUAGL.PS

   1 PROGRAM QUAGL(INPUT,OUTPUT);
   2
   3   VAR   A,B,C,DISKR,X1,X2,RE,IM: REAL;
   4
   5 BEGIN
   6     READ(A,B,C);
   7     DISKR := B*B - 4*A*C;
   8     IF DISKR>=0 THEN
   9                     BEGIN
  10                         X1 := (-B + SQRT(DISKR))/(2*A);
  11                         X2 := (-B - SQRT(DISKR))/(2*A);
  12                         WRITELN("X1=", X1:10:5);
  13                         WRITELN("X2=", X2:10:5)
  14                     END
  15                 ELSE
  16                     BEGIN
  17                         RE := -B/(2*A);
  18                         IM := SQRT(-DISKR)/(2*A);
  19                         WRITELN("X1=", RE:10:5, "     +", IM:10:5, " *J");
  20                         WRITELN("X2=", RE:10:5, "     -", IM:10:5, " *J")
  21                     END
  22 END.