Система уравнений

uses crt;
var
d,d1,d2,d3:real;
a:array[1..3,1..3] of real;
b:array[1..3] of real;
x1,x2,x3:real;
begin
clrscr;
//Матрица
a[1,1]:=2; a[1,2]:=1; a[1,3]:=-1;
a[2,1]:=-2; a[2,2]:=2; a[2,3]:=-3;
a[3,1]:=1; a[3,2]:=-1; a[3,3]:=4;
//Свободные члены
b[1]:=1; b[2]:=2; b[3]:=3;
//Определитель матрицы
d:=a[1,1]*a[2,2]*a[3,3]-a[1,1]*a[2,3]*a[3,2]-a[1,2]*a[2,1]*a[3,3]+a[1,2]*a[2,3]*a[3,1]+a[1,3]*a[2,1]*a[3,2]-a[1,3]*a[2,2]*a[3,1];
//Определитель матрицы с заменой первого столбца столбцом свободных членов
d1:=b[1]*a[2,2]*a[3,3]-b[1]*a[2,3]*a[3,2]-a[1,2]*b[2]*a[3,3]+a[1,2]*a[2,3]*b[3]+a[1,3]*b[2]*a[3,2]-a[1,3]*a[2,2]*b[3];
//Определитель матрицы с заменой второго столбца столбцом свободных членов
d2:=a[1,1]*b[2]*a[3,3]-a[1,1]*a[2,3]*b[3]-b[1]*a[2,1]*a[3,3]+b[1]*a[2,3]*a[3,1]+a[1,3]*a[2,1]*b[3]-a[1,3]*b[2]*a[3,1];
//Определитель матрицы с заменой третьего столбца столбцом свободных членов
d3:=a[1,1]*a[2,2]*b[3]-a[1,1]*b[2]*a[3,2]-a[1,2]*a[2,1]*b[3]+a[1,2]*b[2]*a[3,1]+b[1]*a[2,1]*a[3,2]-b[1]*a[2,2]*a[3,1];
//Корни системы
x1:=d1/d; x2:=d2/d; x3:=d3/d;
writeln('x1=',x1:6:4 ,'x2=':6,x2:6:4,'x3=':6,x3:6:4);
readln
end.

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