We explore the various Mathematica functions for performing Gaussian elimination, specifically when the linear system contains input parameters. On the one hand, we want the solution, perhaps heuristically found, to be in terms of the input parameters, and as general as possible. On the other hand, we must adhere to all necessary conditions and domain specifications. Demonstrating that no Mathematica functions handle all such cases, we build our own function.