Метод гаусса

Есть код в паскале. Он вычисляет систему линейных уравнений методом Гаусса с выбором ведущего элемента по строке. Но когда система несовместима и не имеет решений или имеет бесконечное множество решений код пишет x=NaN. Нужно в этот код добавить строчки, что в таких случаях он не выводил х=NaN, а писал « Система не имеет решений» или « Система имеет бесконечное множество решений» Если код не понятный, можно написать новый, но а тем же смыслом. type matrix = Array [1..3, 1..4] of real; var A: matrix; i,j,k: integer; max: real; b: real; begin writeln('Введите коэффициенты и свободные члены уравнений:'); for i:=1 to 3 do for j:=1 to 4 do readln(A[i,j]); for k:=1 to 3 do begin max:=abs(A[k,k]); i:=k; for j:=k+1 to 3 do if abs(A[j,k])>max then begin max:=abs(A[j,k]); i:=j; end; if ik then for j:=k to 4 do begin b:=A[k,j]; A[k,j]:=A[i,j]; A[i,j]:=b; end; for i:=k+1 to 3 do begin b:=A[i,k]/A[k,k]; for j:=k to 4 do A[i,j]:=A[i,j]-b*A[k,j]; end; end; for i:=3 downto 1 do begin b:=A[i,i]; for j:=i to 4 do A[i,j]:=A[i,j]/b; for j:=i-1 downto 1 do begin b:=A[j,i]; for k:=i to 4 do A[j,k]:=A[j,k]-b*A[i,k]; end; end; writeln('Решения уравнений:'); for i:=1 to 3 do writeln('x', i, ' = ', A[i,4]:6:2); readln; end.