I just thought I'd post this here. It will help you solve a matrix by getting it into row-echelon form (well, close to it anyway). I avoided global variables and tried to make it as expandable as possible.
#include <iostream>
using namespace std;
void output(int height, int width, long...