Программа неправильно решает уравнение

uses MyCub;
 
var
  a, b, c:double;
  f:byte;
begin
  repeat
    writeln('Vvedite 1 dlya resh. uravn ili 0 dlya vihoda');
    readln(f);
    if f=1 then
    begin
      writeln('Vvedite A, B, C');
      readln(a, b, c);
      SolveCub(a, b, c);
    end;
  until(f<>1);
end.

unit MyCub;
 
interface
 
procedure SolveCub(a, b, c:real);
 
implementation
 
function sgn(x:real):integer;
begin
  if x>0 then sgn:=1;
  if x<0 then sgn:=-1;
  if x=0 then sgn:=0;
end;
 
procedure SolveCub(a, b, c:real);
var
  q, p:real;       //коеф. приведен. уравнения
  d:real;          //дискриминант прив. уравн.
  x1, x2, x3:real; //корни уравн
  y1, y2:real;     //вспогательные переменные для вычисления дейст. и мнимой части
  Re, Im:real;     //действит. и мнимая часть комплекного корня
  F:real;          //аргумент угла
begin
  p:=(3*b-a*a)/3;
  q:=(2*a*a*a-9*a*b+27*c)/27;
  d:=(q*q)/4+(p*p*p)/27;
  y1:=0; y2:=0; F:=0;
  if d>0 then
  begin
    y1:=sgn(-q/2+Sqrt(d))*exp(ln(abs(-q/2+Sqrt(d)))/3);
    y2:=sgn(-q/2-Sqrt(d))*exp(ln(abs(-q/2-Sqrt(d)))/3);
    
    x1:=y1+y2-a/3;
    Re:=-(y1+y2)/2-a/3;
    Im:=(y1-y2)*Sqrt(3)/2;
    
    writeln('x1 = ', x1:6:5); 
    writeln('x2 = ', Re:6:5,' + i*', Im:6:5);
    writeln('x2 = ', Re:6:5,' - i*', Im:6:5);
  end;
  
  if d<0 then
  begin
    if q<0 then F:=Arctan(-2*Sqrt(-d)/q);
    if q>0 then F:=Arctan(-2*Sqrt(-d)/q)+Pi;
    if q=0 then F:=Pi/2;
    
    x1:=2*Sqrt(-p/3)*Cos(F/3)-a/3;
    x2:=2*Sqrt(-p/3)*Cos((F+2*Pi)/3)-a/3;
    x3:=2*Sqrt(-p/3)*Cos((F+4*Pi)/3)-a/3;
    if q=0 then x3:=-a/3;
    
    writeln('x1 = ', x1:6:5);
    writeln('x2 = ', x2:6:5);
    writeln('x3 = ', x3:6:5);
  end;
  
  if d=0 then
  begin
    y1:=sgn(-q/2)*exp(ln(abs(-q/2))/3);
    
    x1:=2*y1-a/3;
    x2:=-y1-a/3;
    x3:=-y1-a/3;
    
    writeln('x1 = ', x1:6:5);
    writeln('x2 = ', x2:6:5);
    writeln('x3 = ', x3:6:5);
  end;
end;
 
end.