Lineares Gleichungssystem lösen

[Jens01]

Ich habe das einfach nur in den Editor getippt. Probier mal aus. Ich habe schon lange keine Matrizen mehr lösen müssen.

Delphi-Quellcode:

var
  a, b : TDoubleDynArray;
  m1, m2, mres : IMatrix;
begin
  a := [1,2,3];
  b := [1,2,3,4,5,6,7,8,9];
  
  m1 := TDoubleMatrix.Create(a, 1, 3);
  m2 := TDoubleMatrix.Create(b, 3, 3);  
  
  try
    mres := m2.SolveLinEQ(m1);
  except
    // hier kommt was, wenn die Matrix nicht funktioniert
  end;
  
  write(mres[0,0].tostring);
  write(mres[0,1].tostring);
  write(mres[0,2].tostring);
end;