function xstar = linearSolve1(a,b)
if a == 0
   xstar = Inf*b; %! No solution exists.
else
   xstar = b./a;
end
