Решение нелинейного уравнения и системы методом Гаусса

uses crt;
function F(x:real):real;
begin
   F:=x+sqrt(x)+exp(ln(x)/3)-2.5;
end;
var eps,f1,f2,x,x1,x2,a,b,c:real;
    n,i:integer;
begin
clrscr;
writeln('Нахождение корня уравнения x+x^1/2+x^1/3-2.5':55);
writeln('на интервале [0.4;1]':40);
writeln('методами итераций и дихотомии':45);
writeln;
writeln;
a:=0.4;
b:=1;
eps:=0.0001;
{метод итераций}
n:=round((b-a)/eps);
x1:=a; f1:=F(x1);
for i:=1 to n do
  begin
   x2:=x1+(b-a)/n;
   f2:=F(x2);
   if f1*f2 <= 0 then writeln('Метод итераций  x = ', (x2 + x1)/2:0:5);
   x1:=x2;
   f1:=f2;
  end;
{метод деления пополам-дихотомии}
repeat
 c:=(a+b)/2;
 if F(a)*F(c)<=0 then b:=c
 else a:=c;
until b-a<eps;
x:=(a+b)/2;
writeln('Метод дихотомии x = ',x:0:5);
readln
End.

uses crt;
const n=3;
var a1,b1,c1,//первое уравнение
    a2,b2,c2,//второе
    x1,x2:real;//корни
begin
clrscr;
a1:=3;b1:=4;c1:=-1;
a2:=2;b2:=-5;c2:=30;
{приведение матрицы к треугольному виду}
{обнуляем коэффициент при х1 во 2 уравнении}
b2:=b1+b2*(-a1/a2); c2:=c1+c2*(-a1/a2);
{ищем корни, начиная с последнего}
x2:=c2/b2;
x1:=(c1-b1*x2)/a1;
writeln('Корни системы:');
write('x1=',x1:0:2,'  x2=',x2:0:2);
readln;
end.