function [xstar,ystar] = linearSolve2(a,b,c,d,e,f)
if a*e - b*d == 0 
   xstar = Inf; ystar = Inf; %! No solution exists.
else
   xstar = (e*c - b*f) ./ (a*e - b*d);
   if e == 0
      ystar = c./b - xstar.*a./b;
   else
      ystar = f./e - xstar.*d./e;
   end
end
