(*$B- Algorithmus 2.1*)
PROGRAM quad;
VAR p,q,x1,x2,diskr,fak : real;
BEGIN
  writeln('p,q eingeben'); read(p,q);
  IF abs(p)>1 THEN
  BEGIN fak:=abs(p); diskr:=0.25-q/p/p END
  ELSE
  BEGIN fak:=1; diskr:=sqr(p/2)-q END;
  IF diskr < 0 THEN
  writeln('komplexe Loesungen',-p/2,'+-I*',fak*sqrt(-diskr))
  ELSE
  BEGIN
    x1:=abs(p/2)+fak*sqrt(diskr);
    IF p>0 THEN x1:=-x1;
    IF x1=0 THEN x2:=0 ELSE x2:=q/x1;
    writeln('reelle Loesungen',x1,x2)
  END
END.