class BjerksundStenslandModel {

  protected double[] getKappaAdjoint(double gamma, double b, double sigmaSq) {
    double[] res = new double[4];
    double temp = 2 * b / sigmaSq;
    res[0] = temp + 2 * gamma - 1;
    res[1] = 2.0; // gammaBar
    res[2] = 2 / sigmaSq; // bBar
    res[3] = -temp / sigmaSq; // sigmasqBar
    return res;
  }

}