Error: Illegal qualifier

 uses graph, LABA4, crt;
 const
  shag_setki=7;   {kolichestvo liniy v sisteme koordinat}
  otst=80;
  k=8.62e-5;
  n=1.45e10;
 var
  Nd,Na,a,T,b,sh:real;
  U,Tabs: array [0..25] of real;
  i,  r : integer;
  grDriver: integer;
  grMode: integer;
  ymax, xmax:integer;
  sh_y,sh_x,                   {shag setki po y & x}
  dx,dy:integer;               {delta x & delta y}
  T_max,T_min,U_max,U_min,     {maximalnie & minimalnie znacheniya n & E}
  dT,dU,                       {delta n & delta E}
  kx,ky:real;                  {koeficienti perevoda v pixeli}
  y,x:array [0..50] of integer;{massivi dlya postroeniya tochek}
  Tstr,Ustr: string;           {peremennie dlya preobrazovaniya znacheniy n & E v stroki}
begin
  clrscr;
  writeln('Enter the minimum temperature value in Kelvins ranging from 100 to 500');
  readln(a);
  writeln('Enter the maximum temperature value in Kelvins ranging from 100 to 500');
  readln(b);
  writeln('Enter step');
  readln(sh);
  writeln('Enter legitima impurity concentration in the range of from 5e13 to 1e18');
  readln(Na);
  Nd:=Na;
  T:=a;
  r:=trunc((b-a)/sh);
  for i:=0 to r do
    begin
      Tabs[i] := T;
      U[i] := UrovenFermi(T);
      T:=T+sh;
    end;
    Tablica(r,U,tabs);
    begin
    U[i]:=k*T*ln((Na*Nd)/(n*n));
    end;
    readln;
                  {nahozhdenie maximalnih & minimalnih znacheniy T & U}
   begin
    for i:=0 to r do
     begin
      if T[r] <= T[i] then T_max:= T[i] else T_min:=T[i];
      if T[0] < T[i] then T_min:= T[0] else T_max:= T[0];
      if U[0] <= U[i] then U_max:= U[i] else U_min:=U[i];
      if U[r] < U[i] then U_min:= U[0] else U_max:= U[0];
     end;
    dT:=T_max-T_min;
    dU:=U_max-U_min;
   end;
                          {vkluchenie graficheskogo rezhima}
  grDriver:=VGA;
  grMode:= VGAHi;
  InitGraph(grDriver,grMode,'');
                            {sozdanie koordinatnoy setki}
  ymax:= GetMaxY;
  xmax:= GetMaxX;
  sh_y:= round((ymax-otst-otst)/shag_setki);
  sh_x:= round((xmax-otst-otst)/shag_setki);
  SetLineStyle(DottedLn,2,NormWidth);
  SetColor(9);
    for i:=0 to shag_setki do
     begin
      line(otst,ymax-otst-sh_y*i,xmax-otst,ymax-otst-sh_y*i);
      line(otst+sh_x*i,ymax-otst,otst+sh_x*i,ymax-otst-sh_y*shag_setki);
     end;
                                {sozdanie ramki}
  SetLineStyle(SolidLn,0,NormWidth);
  SetColor(15);
   Rectangle(otst,ymax-otst-sh_y*shag_setki,otst+sh_x*shag_setki,ymax-otst);
                            {sozdanie podpisey na osyah}
  Setcolor(10);
   SetTextJustify(1,2);
    OutTextXY(otst+sh_x*shag_setki,ymax-otst+20, 'T, K');
   SetTextJustify(1,0);
    OutTextXY(otst,ymax-otst-sh_y*shag_setki-10, 'Fi, B');
                         {preobrazovanie znacheniy v stroki}
  SetColor(2);
  for  i:=0 to shag_setki do
     begin
      str(T[0]+dT/shag_setki*i:6,Tstr);
       SetTextJustify(1,1);
        OutTextXY(otst+sh_x*i,ymax-otst+15,Tstr);
      str(U[0]-dU/shag_setki*i:2:4,Ustr);
       SetTextJustify(2,1);
        OutTextXY(otst-2,ymax-otst-sh_y*shag_setki+sh_y*i,Ustr);
     end;
                           {sozdanie nazvaniya grafika}
  SetColor(3);
  SetTextJustify(1,2);
   OutTextXY(round(xmax/2),round(otst/2),'A graph of the built-in potential temperature');
                            {sozdanie tochek na grafike}
  dx:=sh_x*shag_setki;
  dy:=sh_y*shag_setki;
  kx:=dx/dT;
  ky:=dy/dU;
  SetColor(2);
   for i:=0 to r do
    begin
     x[i]:=trunc(otst+(T[i]-T_min)*kx);
     y[i]:=trunc(ymax-otst-sh_y*shag_setki+(U_max-U[i])*ky);
     circle(x[i],y[i],3);
    end;
                                   {sozdanie liniy}
  setcolor(15);
    r:=r-1;
    for i:=0 to r do
     begin
      line (x[i],y[i],x[i+1],y[i+1]);
     end;
   readln;
   closegraph;
end.

Unit LABA4;
  Interface
    function UrovenFermi(Tt: real): real;
    procedure Tablica(rr:integer; Uu:array of real; Ttabs: array of real);
  Implementation
    function UrovenFermi(Tt: real): real;
      const
       k=8.62e-5;
       Na=7e15;
       n=1.45e10;
      begin
       UrovenFermi:=k*Tt*ln((Na*Na)/(n*n));
      end;
    procedure Tablica(rr:integer; Uu:array of real; ttabs:array of real);
      var j:integer;
      begin
           writeln('Temperature?, [K] The built-in potential, [B]');
           for j:=0 to rr do
              begin
              writeln('     ',Ttabs[j]:5:1 ,'                  ',Uu[j]:3:2);
              end;
      end;
  end.