Сортировка строк в матрице

uses
  Printer;
Var x1,x2,dx,a,b,c,E:real;
    pa,pb,pc:string;
 
Function f(x:real):real;
Begin
    f:=(2*x)*sin(x)-cos(x);
end;
 
Function f1(x:real):real;
Begin
    f1:=x*cos(x)+sin(x);
end;
 
Begin{Основная программа}
 
{Метод Проб}
     Writeln('Введите [a,b], E=');
     readln(a,b,E);
     Writeln('Отрезок:');
     Writeln('[',a:5:1,',',b:5:1,']');
     Writeln('Точность вычисления:');
     Writeln(E:6:4);
     Repeat
     if f(a)>0 then pa:='+' else pa:='-';
     if f(b)>0 then pb:='+' else pb:='-';
     c:=(a+b)/2;
     if f(c)>0 then pc:='+' else pc:='-';
     if pc=pa then a:=c else b:=c;
     Until abs(b-a)<=E;
     Writeln('Метод Проб:','Корень=',c:5:4);
 
{Метод Ньютона}
     if pa='da' then x1:=b else x1:=a;
     Repeat
     dx:=(-(f(x1)/f1(x1)));
     x1:=x1+dx;
     Until abs(dx)<=E;
     Writeln('Метод Ньютона:', 'Корень=',x1:5:4);
 
{Метод Хорд}
     x1:=a; x2:=b;
     Repeat
     c:=(x1-x2)/(f(x1)-f(x2));
     if pa='da'then 
       Begin
        dx:=(-f(x2))*c;
        x2:=x2+dx;
       end
      else 
       Begin
        dx:=(-f(x1))*c;
        x1:=x1+dx;
       end;
     Until abs(dx)<=0.001;
     if pa='da' then Writeln('Метод Хорд:','Корень=',x2:5:4)
      else Writeln('Метод Хорд:','Корень=',x1:5:4);
  readln;
end.