Позволяет решить неоднородную систему n линейных алгебраических уравнений с n неизвестными, если определитель основной матрицы не равен нулю.

Code

{$N+} { Задаем порядок системы уравнений } Const n = 3; Type   { Тип, описывающий матрицу системы (включая свободные члены !!!) }   Equation = Array[1 .. n, 1 .. Succ(n)] Of Double;   matrix =    Array[1 .. n, 1 .. n] Of Double; Const   a: Equation =   ((2, 1, 3, 9),    (1, -2, 1, -2),    (3, 2, 2, 7)); { Процедура замены очередного столбца матрицы на свободные члены } Procedure GetMatrix(wout: Integer; Var m: matrix);   Var i, j: Integer;   Begin    For i := 1 To n Do     For j := 1 To n Do      If j <> wout Then m[i, j] := a[i, j]      Else m[i, j] := a[i, Succ(n)]   End; { Нахождение определителя } Function Det(a: matrix; n: integer): Double;   Var i, j, k: Integer; d: Double;   Const    Eps = 10E-6;   Begin    For i := 1 To Pred(n) Do     Begin      If Abs(a[i, i]) < Eps Then       Begin        Det := 0.0; Exit       End;      For j := Succ(i) To n Do       Begin        d := a[j, i] / a[i, i];        For k := i To n Do         a[j, k] := a[j, k] - d * a[i, k];       End;     End;    d := 1.0;    For i := 1 To n Do     d := d * a[i, i];    Det := d   End; Var   i: Integer;   mx: matrix;   Determ: Double; begin   GetMatrix(Succ(n), mx);   Determ := Det(mx, n);   If Abs(Determ) < 1E-6 Then    Writeln( 'Определитель исходной матрицы = 0' )   Else    For i := 1 To n Do     Begin      GetMatrix(i, mx);      WriteLn( 'x(', i, ') = ', (Det(mx, n) / Determ):7:4 )     End

Это оно?
Кстати очень интересный форум называется Всё о Паскале
Только это кажись паскаль... а не С++