template <typename TwoDArray>
void GaussJordan(const int& n, TwoDArray& a)
{
	for (int i=0; i<n; i++)
	{
		if(a[i][i] == 0.0)
		{
			cout << "Error, zero pivot values found!";
			exit(1);
		}
		
		for (int j=0; j<n; j++)
		{
			if (i!=j)
			{
				double ratio = a[j][i]/a[i][i];
				for (int k=0; k<n+1; k++)
				{
					a[j][k] = a[j][k] - ratio * a[i][k];
				}
			}
		}
	}
	
	for (int i=0; i<n; i++)
	{
		a[i][n] /= a[i][i];
		a[i][i] = 1.0;
	}
}