docs: private double ssqerr(int k0, int kmax, double[] logk, double[] log_kDist, double m, double t) {
      int k = kmax - k0;
      double result = 0;
      for(int i = 0; i < k; i++) {
        // double h = log_kDist[i] - (m * (logk[i] - logk[0]) + t); ???
        double h = log_kDist[i] - m * logk[i] - t;
        result += h * h;
      }
      return result;
    }