docs: private double joint(int k, int y[]) {
          double p = 1.0; //priors[k];
          for (int j = 0; j < y.length; j++) {
              if (j != k && y[j] == 1) {
                  p *= Conditional[k][j];
              }
          }
          return p;
      }