Start generating basis
Add test case
Free vars
Split Gaussian elimination into full and partial pivoting
Use partial Gaussian in Gauss-Jordan
Changes to LU
Changes to probset
New basis
Generate Vandermonde matrix