double cutOff = 0.10;
double mu = 2.5;
double strike = 0.15;
double t = 2.366105247;
EuropeanVanillaOption option = new EuropeanVanillaOption(strike, t, true);
SABRFormulaData sabrData = new SABRFormulaData(alpha, beta, rho, nu);
double forward = 0.0404500579038675;
SABRExtrapolationRightFunction sabrExtrapolation = new SABRExtrapolationRightFunction(forward, sabrData, cutOff, t, mu);
double shift = 0.000001;
SABRFormulaData sabrDataAP = new SABRFormulaData(alpha + shift, beta, rho, nu);
SABRFormulaData sabrDataBP = new SABRFormulaData(alpha, beta + shift, rho, nu);
SABRFormulaData sabrDataRP = new SABRFormulaData(alpha, beta, rho + shift, nu);
SABRFormulaData sabrDataNP = new SABRFormulaData(alpha, beta, rho, nu + shift);
SABRExtrapolationRightFunction sabrExtrapolationAP = new SABRExtrapolationRightFunction(forward, sabrDataAP, cutOff, t, mu);
SABRExtrapolationRightFunction sabrExtrapolationBP = new SABRExtrapolationRightFunction(forward, sabrDataBP, cutOff, t, mu);
SABRExtrapolationRightFunction sabrExtrapolationRP = new SABRExtrapolationRightFunction(forward, sabrDataRP, cutOff, t, mu);
SABRExtrapolationRightFunction sabrExtrapolationNP = new SABRExtrapolationRightFunction(forward, sabrDataNP, cutOff, t, mu);
// Above cut-off strike