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