Calculos dan NaN

[Instru]

Hola.
Tengo un problema. Tengo una función para hacer un Gauss -Jordan.
Pero por alguna razon al final la matriz termina con puros NaN(not a number).
alguien sabe por que sucede esto?

Código:

function gaussJordan()
{
	var N, M, i, j, k, p;
	var mat;
	var a=0.0;
	
	mat=[
			[0.8875, , 0.7875, 0],
			[-0.1125, 0.2125, 0]
		];
	M=mat.length;
	N=M+1;
	alert("M="+M);
	
	//For each row
	for(i=0; i<M; i++)
	{
		//Divide the row by its pivot
		for(j=0; j<N; j++)
		{
			a=mat[i][i];
			mat[i][j]/=a;
			
		}
		//Make the variables above and under the pivot 0
		for(k=0; k<M; k++)
		{
			for(p=0; p<N; p++)
			{
				if(k!=i)
				{
					alert(mat);
					a=mat[k][p]-mat[i][p]*mat[k][i];
					mat[k][p]=a;
				}
			}
		}
	}
	return mat;
}