code: public static Matrix identity(int N) {
          Matrix I = new Matrix(N, N);
          for (int i = 0; i < N; i++) {
              I.data[i][i] = 1;
          }
          return I;
      }
docs: create and return the N-by-N identity matrix