docs: private static double QUARD_MINIMIZER(double qm, double u, double fu, double du, double v, double fv) {
          double a, d, gamma, theta, p, q, r, s;
          a = (v) - (u);
          (qm) = (u) + (du) / (((fu) - (fv)) / a + (du)) / 2 * a;
          return qm;
      }