С помощью матрицы - значит, методом Крамера.
Delta = a11*a22 - a12*a21;
Deltax = b1*a22 - b2*a12;
Deltay = a11*b2 - a21*b1;
if ((Delta = 0) and (Deltax = 0) and (Deltay = 0)) then write "Решений бесконечно много";
if ((Delta = 0) and ((Deltax <> 0) or (Deltay <> 0)) then write "Решений нет";
if (Delta <> 0) then
begin
x = Deltax / Delta;
y = Deltay / Delta;
end;