Delphi-Quellcode:
procedure Inverse(const N: integer; var A: TMatrixGauss);
//...
A[I, J]:= A[I, J] / A[I, I]; // Matrix singulär, falls A[I, I] = 0
//...
end;
Der Kommentar ist falsch bzw. unsinnig/irreführend. Wenn A[I, I] gleich 0 ist, heißt das erstmal nur, daß man ohne Zeilenvertauschungen nicht weiter kommt. Selbstverständlich ist zB die folgende Matrix
nicht singulär: