class BjerksundStenslandModel {

  protected double[] getLambdaAdjoint(double gamma, double r, double b, double sigmaSq) {
    double[] res = new double[5];
    double temp = 0.5 * gamma * (gamma - 1);
    res[0] = -r + gamma * b + temp * sigmaSq; // lambda
    res[1] = b + (gamma - 0.5) * sigmaSq; // gammaBar
    res[2] = -1.0; // rBar
    res[3] = gamma; // bBar
    res[4] = temp; // sigmasqBar
    return res;
  }

}