二元一次方程 高斯消元法計算函數

  TRS=array[1..2] of Real;

{二元一次方程 高斯消元法計算函數}
function X1Y1_2(const ia1,ib1,ic1,ia2,ib2,ic2:Real):TRS;
const //爲數組的行數與列數
  Rows = 2;
  Cols = 3;
var
  C: array[1..Rows, 1..Cols] of Real; //係數及常數項
  R1, R2: Real;
  m,n,i: Integer;
begin
  C[1,1]:=ia1;
  C[1,2]:=ib1;
  C[1,3]:=ic1;
  C[2,1]:=ia2;
  C[2,2]:=ib2;
  C[2,3]:=ic2;

  //高斯消元
  for n := 1 to Rows do
  begin
    R1 := C[n,n];
    for m := 1 to Cols do C[n,m] := C[n,m] / R1; //使每個方程對角線上的係數爲1
    for i := 1 to Rows do
      if i <> n then //依次對其它方程的第n項進行消元
      begin
        R2 := C[i,n];
        for m := 1 to Cols do C[i,m] := C[i,m] - C[n,m]*R2;
      end;
  end;
  Result[1]:=C[1,3];
  Result[2]:=C[2,3];
end;