Квадратное уравнение с комплексными числами

uses crt;
type tcomplex=record
              re,im:real;
              end;
var a,b,c,d,f,g,h,x1,x2:tcomplex;
begin
clrscr;
writeln('Ââåäèòå êîýôôèöèåíòû óðàâíåíèÿ:');
write('äåéñòâèòåëüíàÿ ÷àñòü A=');
readln(a.re);
write('ìíèìàÿ ÷àñòü A=');
readln(a.im);
write('äåéñòâèòåëüíàÿ ÷àñòü B=');
readln(b.re);
write('ìíèìàÿ ÷àñòü B=');
readln(b.im);
write('äåéñòâèòåëüíàÿ ÷àñòü C=');
readln(c.re);
write('ìíèìàÿ ÷àñòü C=');
readln(c.im);
clrscr;
writeln('Óðàâíåíèå:');
if a.im<0 then write('(',a.re:0:2,a.im:0:2,'i)')
else write('(',a.re:0:2,' +',a.im:0:2,'i)');
write('x^2');
if b.re<0 then write('-(',-b.re:0:2)else write('+(',b.re:0:2);
if b.im<0 then write(b.im:0:2,'i')else write('+',b.im:0:2,'i)');
write('x');
if c.re<0 then write('-(',-c.re:0:2)else write('+(',c.re:0:2);
if c.im<0 then write(c.im:0:2,'i)')else write('+',b.im:0:2,'i)');
f.re:=b.re*b.re-b.im*b.im;{b^2}
f.im:=b.re*b.im+b.re*b.im;
g.re:=4*a.re*c.re-4*a.im*c.im;{4ac}
g.im:=4*a.re*c.im+4*a.re*c.im;
d.re:=f.re-g.re;{discr.}
d.im:=f.im-g.im;
writeln('=0');
readln
end.

uses crt;
var Rea,Ima,Reb,Imb,Rec,Imc,Rek,Imk,Rel,Iml,Rem,Imm,Ren,Imn,
    Red,Imd,Red1,Imd1,Red2,Imd2,Reo1,Imp1,Reo2,Imp2,Reo3,Imp3,Reo4,Imp4,Rex1,Rex2,Rex3,Rex4,
    Imx1,Imx2,Imx3,Imx4,sgnImn,q1,q11,q2,q22,q3,q33,q4,q44,w1,w11,w2,w22,w3,w33,w4,w44:real;
begin
clrscr;
write('Vvedit Re a:');readln(Rea);
write('Vvedit Im a');readln(Ima);
write('Vvedit Re b:');readln(Reb);
write('Vvedit Im b');readln(Imb);
write('Vvedit Re c:');readln(Rec);
write('Vvedit Im c');readln(Imc);
Rek:=Reb*Reb-Imb*Imb;
Imk:=2*Reb*Imb;
Rel:=4*(Ima*Imc-Rea*Rec);
Iml:=-4*(Rea*Imc+Rec*Ima);
Ren:=Rek+Rel;
Imn:=Imk+Iml;
if Imn>=0 then sgnImn:=1
else sgnImn:=-1;
Red1:= sqrt((sqrt(Ren*Ren+Imn*Imn)+Ren)/2);
Red2:=-sqrt((sqrt(Ren*Ren+Imn*Imn)+Ren)/2);
Imd1:=sgnImn*sqrt((sqrt(Ren*Ren+Imn*Imn)-Ren)/2);
Imd2:=-sgnImn*sqrt((sqrt(Ren*Ren+Imn*Imn)-Ren)/2);
Reo1:=Red1-Reb;
Imp1:=Imd1-Imb;
Reo2:=Red2-Reb;
Imp2:=Imd2-Imb;
Reo3:=-(Red1+Reb);
Imp3:=-(Imd1+Imb);
Reo4:=-(Red2+Reb);
Imp4:=-(Imd2+Imb);
w1:=Imp1*Rea/(2*(Rea*Rea+Ima*Ima));
w11:=Reo1*Ima/(2*(Rea*Rea+Ima*Ima));
q1:=Reo1*Rea/(2*(Rea*Rea+Ima*Ima));
q11:=Imp1*Ima/(2*(Rea*Rea+Ima*Ima));
Rex1:=q1+q2;
Imx1:=w1-w2;
q2:=Reo2*Rea/(2*(Rea*Rea+Ima*Ima));
q22:=Imp2*Ima/(2*(Rea*Rea+Ima*Ima));
Rex2:=q2-q22;
w2:=Imp2*Rea/(2*(Rea*Rea+Ima*Ima));
w22:=Reo2*Ima/(2*(Rea*Rea+Ima*Ima));
Imx2:=w2-w22;
q3:=Reo3*Rea/(2*(Rea*Rea+Ima*Ima));
q33:=Imp3*Ima/(2*(Rea*Rea+Ima*Ima));
Rex3:=q3+q33;
w3:=Imp3*Rea/(2*(Rea*Rea+Ima*Ima));
w33:=Reo3*Ima/(2*(Rea*Rea+Ima*Ima));
Imx3:=w3-w33;
q4:=Reo4*Rea/(2*(Rea*Rea+Ima*Ima));
q44:=Imp4*Ima/(2*(Rea*Rea+Ima*Ima));
Rex4:=q4-q44;
w4:=Imp4*Rea/(2*(Rea*Rea+Ima*Ima));
w44:=Reo4*Ima/(2*(Rea*Rea+Ima*Ima));
Imx4:=w4-w44;
if Imx1>0 then writeln('X1= ',Rex1:0:2,' +',Imx1:0:2,'i')
else writeln('X1= ',Rex1:0:2,' ',Imx1:0:2,'i');
if Imx2>0 then writeln('X2= ',Rex2:0:2,' +',Imx2:0:2,'i')
else writeln('X2= ',Rex2:0:2,Imx2:0:2,'i');
if Imx3>0 then writeln('X3= ',Rex3:0:2,' +',Imx3:0:2,'i')
else writeln('X3= ',Rex3:0:2,Imx3:0:2,'i');
if Imx4>0 then writeln('X4= ',' ',Rex4:0:2,' +',Imx4:0:2,'i')
else writeln('X4= ',Rex4:0:2,' ',Imx4:0:2,'i');
readln
end.