Найти комплексные корни уравнения

       program pr;
       uses crt;
         var t1:text;
         const
            a0=1E-5;
            ro1= 1000;
            ro2= 0.590;
            Cp1= 4216;
            //Cp2= 2034;
            lyam1= 0.680;
            gam= 1.3;
            l= 2257200;
            pi=3.14;
            p0= 1e5;
            T111=4967;
            p111=6.0716e+10;
            sigma= 0.05891;
            v= 279e-9;
            Rv=454;
            Rg=287;
            E=1e-8;
 
        var  mg0,mg00,rov,Cp2,Z,A,value1,vT,a00,pg0,lam_TT,Lam,Lam_R,Lam_VV,Lam_T,lam_RR,lam_MM,pv0,max,Ya,betta,Ts,Ts1,Nu,laym,laym1,value:real;
               i:integer;
 
          function f(a01,mg01:real):real;
          begin
               f:= (4*pi/3)*((pv0-p0)*a01*a01*a01-2*sigma*a01*a01)+mg01*Rg*Ts1;
          end;
 
 
          function Y(lam,a002:real):real;
          begin
               Y:= sqrt((lam*a002*a002*ro1*Cp1)/lyam1);
          end;
          
          function pg00(mg02,a002:real):real;
          begin
               pg00:=(3*mg0*Rg*Ts1)/(4*pi*a002*a002*a002);
          end;
          
          function FF(lam,a02,mg02:real):real;
          begin
              FF:= ro1*a02*a02*lam*lam+4*v*ro1*lam+(3*gam*pv0*Y(lam,a02)*Y(lam,a02)*(1+pg00(mg02,a02)/(ro2*l)))/
              (gam*Y(lam,a02)*Y(lam,a02)+betta*(1+Y(lam,a02)))+3*pg00(mg02,a02)-2*sigma/a02;
          end;
 
 
          function  pol_del1(mg_01:real):real;
           var a,b,c: real;
             begin
               a:= 0;
               b:= 6.6e-6;
               c:= (a+b)/2;
                 while (abs(b-a)>E) do
                       begin
                           if f(a,mg_01)*f(c,mg_01) < 0 then b:= c
                                    else a:= c;
                                 c:= (a+b)/2;
                       end;
                   value1:= (a+b)/2;
                   result:=value1;
             end;
 
          function  pol_del2(res,mg_02:real):real;
           var a,b,c: real;
             begin
              a:= -1e6;
               b:= 1e6;
               c:= (a+b)/2;
                 while (abs(b-a)>E) do
                       begin
                           if FF(a,res,mg_02)*FF(c,res,mg_02) < 0 then b:= c
                                    else a:= c;
                                 c:= (a+b)/2;
                       end;
                   Lam:= (a+b)/2;
                   result:=Lam;
             end;
 
 
        begin
        
            pv0:=p0+2*sigma/a0;
            vT:=lyam1/(ro1*Cp1);
            Ts1:=T111/ln(p111/pv0);
            mg00:=(128*sigma*sigma*sigma*pi)/(81*Rg*Ts1*(2*sigma/a0)*(2*sigma/a0));
            Cp2:=gam*Rv/(gam-1);
            betta:=3*(gam-1)*sqr(Cp2*Ts1/l)*((ro1*Cp1)/(ro2*Cp2));
 
 
            assign(t1,'D:\äàííûå\10 ìàÿ äëÿ óñòîé÷èâîãî.txt');
            rewrite(t1);
            writeln(t1,'   mg0                     value1                  Lam');
 
            begin
           
              for i:=0 to 1000 do begin
                        mg0:=i*mg00/1000;
                        value1:=pol_del1(mg0);
                        Lam:=pol_del2(value1,mg0);
 
                        writeln(t1,mg0,'    ',value1,'     ',Lam);
              end;
 
           close(t1);
          end;
          end.

type compl=record
                re,in:real;
                end;

Uses Complex;//это наш откомпилированный модуль, файл Complex.tpu, 
который должен быть в папке с программой