public static double Trace(double[][] m) { double result = 0.0d; for (int i = 0; i < m.Length; i++) { result += m[i][i]; } return result; }