net.finmath.functions

Class AnalyticFormulas

java.lang.Object

net.finmath.functions.AnalyticFormulas

public class AnalyticFormulas
extends Object

This class implements some functions as static class methods. It provides functions like

• the Black-Scholes formula,
• the inverse of the Back-Scholes formula with respect to (implied) volatility,
• the Bachelier formula,
• the inverse of the Bachelier formula with respect to (implied) volatility,
• the corresponding functions (versions) for caplets and swaptions,
• analytic approximation for European options under the SABR model,
• some convexity adjustments.

Version:

1.10

Author:

Christian Fries

Date:

27.04.2012

Method Summary

Modifier and Type	Method and Description
static double	bachelierGeneralizedOptionVega(double forward, double volatility, double optionMaturity, double optionStrike, double payoffUnit) Calculates the vega of a call, i.e., the payoff max(S(T)-K,0) P, where S follows a normal process with constant volatility, i.e., a Bachelier model \[ \mathrm{d} S(t) = r S(t) \mathrm{d} t + \sigma \mathrm{d}W(t) \]
static double	bachelierOptionDelta(double forward, double volatility, double optionMaturity, double optionStrike, double payoffUnit) Calculates the option delta of a call, i.e., the payoff max(S(T)-K,0), where S follows a normal process with constant volatility, i.e., a Bachelier model \[ \mathrm{d} S(t) = r S(t) \mathrm{d} t + \sigma \mathrm{d}W(t) \]
static double	bachelierOptionImpliedVolatility(double forward, double optionMaturity, double optionStrike, double payoffUnit, double optionValue) Calculates the Bachelier option implied volatility of a call, i.e., the payoff max(S(T)-K,0), where S follows a normal process with constant volatility.
static double	bachelierOptionValue(double forward, double volatility, double optionMaturity, double optionStrike, double payoffUnit) Calculates the option value of a call, i.e., the payoff max(S(T)-K,0), where S follows a normal process with constant volatility, i.e., a Bachelier model \[ \mathrm{d} S(t) = r S(t) \mathrm{d} t + \sigma \mathrm{d}W(t) \]
static RandomVariable	bachelierOptionValue(RandomVariable forward, RandomVariable volatility, double optionMaturity, double optionStrike, RandomVariable payoffUnit) Calculates the option value of a call, i.e., the payoff max(S(T)-K,0) P, where S follows a normal process with constant volatility, i.e., a Bachelier model \[ \mathrm{d} S(t) = r S(t) \mathrm{d} t + \sigma \mathrm{d}W(t) \]
static double	blackModelCapletValue(double forward, double volatility, double optionMaturity, double optionStrike, double periodLength, double discountFactor) Calculate the value of a caplet assuming the Black'76 model.
static double	blackModelDgitialCapletValue(double forward, double volatility, double periodLength, double discountFactor, double optionMaturity, double optionStrike) Calculate the value of a digital caplet assuming the Black'76 model.
static double	blackModelSwaptionValue(double forwardSwaprate, double volatility, double optionMaturity, double optionStrike, double swapAnnuity) Calculate the value of a swaption assuming the Black'76 model.
static double	blackScholesATMOptionValue(double volatility, double optionMaturity, double forward, double payoffUnit) Calculates the Black-Scholes option value of an atm call option.
static double	blackScholesDigitalOptionDelta(double initialStockValue, double riskFreeRate, double volatility, double optionMaturity, double optionStrike) Calculates the delta of a digital option under a Black-Scholes model
static double	blackScholesDigitalOptionRho(double initialStockValue, double riskFreeRate, double volatility, double optionMaturity, double optionStrike) Calculates the rho of a digital option under a Black-Scholes model
static double	blackScholesDigitalOptionValue(double initialStockValue, double riskFreeRate, double volatility, double optionMaturity, double optionStrike) Calculates the Black-Scholes option value of a digital call option.
static double	blackScholesDigitalOptionVega(double initialStockValue, double riskFreeRate, double volatility, double optionMaturity, double optionStrike) Calculates the vega of a digital option under a Black-Scholes model
static double	blackScholesGeneralizedOptionValue(double forward, double volatility, double optionMaturity, double optionStrike, double payoffUnit) Calculates the Black-Scholes option value of a call, i.e., the payoff max(S(T)-K,0) P, where S follows a log-normal process with constant log-volatility.
static RandomVariable	blackScholesGeneralizedOptionValue(RandomVariable forward, RandomVariable volatility, double optionMaturity, double optionStrike, RandomVariable payoffUnit) Calculates the Black-Scholes option value of a call, i.e., the payoff max(S(T)-K,0) P, where S follows a log-normal process with constant log-volatility.
static double	blackScholesGeneralizedOptionVega(double forward, double volatility, double optionMaturity, double optionStrike, double payoffUnit) Calculates the vega of a call, i.e., the payoff max(S(T)-K,0) P, where S follows a normal process with constant volatility, i.e., a Black-Scholes model \[ \mathrm{d} S(t) = r S(t) \mathrm{d} t + \sigma S(t)\mathrm{d}W(t) \]
static double	blackScholesOptionDelta(double initialStockValue, double riskFreeRate, double volatility, double optionMaturity, double optionStrike) Calculates the delta of a call option under a Black-Scholes model The method also handles cases where the forward and/or option strike is negative and some limit cases where the forward or the option strike is zero.
static RandomVariable	blackScholesOptionDelta(RandomVariable initialStockValue, RandomVariable riskFreeRate, RandomVariable volatility, double optionMaturity, double optionStrike) Calculates the delta of a call option under a Black-Scholes model The method also handles cases where the forward and/or option strike is negative and some limit cases where the forward or the option strike is zero.
static RandomVariable	blackScholesOptionDelta(RandomVariable initialStockValue, RandomVariable riskFreeRate, RandomVariable volatility, double optionMaturity, RandomVariable optionStrike) Calculates the delta of a call option under a Black-Scholes model The method also handles cases where the forward and/or option strike is negative and some limit cases where the forward or the option strike is zero.
static double	blackScholesOptionGamma(double initialStockValue, double riskFreeRate, double volatility, double optionMaturity, double optionStrike) This static method calculated the gamma of a call option under a Black-Scholes model
static RandomVariable	blackScholesOptionGamma(RandomVariable initialStockValue, RandomVariable riskFreeRate, RandomVariable volatility, double optionMaturity, double optionStrike) This static method calculated the gamma of a call option under a Black-Scholes model
static double	blackScholesOptionImpliedVolatility(double forward, double optionMaturity, double optionStrike, double payoffUnit, double optionValue) Calculates the Black-Scholes option implied volatility of a call, i.e., the payoff max(S(T)-K,0), where S follows a log-normal process with constant log-volatility.
static double	blackScholesOptionRho(double initialStockValue, double riskFreeRate, double volatility, double optionMaturity, double optionStrike) This static method calculated the rho of a call option under a Black-Scholes model
static double	blackScholesOptionTheta(double initialStockValue, double riskFreeRate, double volatility, double optionMaturity, double optionStrike) This static method calculated the vega of a call option under a Black-Scholes model
static double	blackScholesOptionValue(double initialStockValue, double riskFreeRate, double volatility, double optionMaturity, double optionStrike) Calculates the Black-Scholes option value of a call, i.e., the payoff max(S(T)-K,0), where S follows a log-normal process with constant log-volatility.
static double	blackScholesOptionValue(double initialStockValue, double riskFreeRate, double volatility, double optionMaturity, double optionStrike, boolean isCall) Calculates the Black-Scholes option value of a call, i.e., the payoff max(S(T)-K,0), or a put, i.e., the payoff max(K-S(T),0), where S follows a log-normal process with constant log-volatility.
static double	blackScholesOptionVega(double initialStockValue, double riskFreeRate, double volatility, double optionMaturity, double optionStrike) Calculates the vega of a call, i.e., the payoff max(S(T)-K,0) P, where S follows a normal process with constant volatility, i.e., a Black-Scholes model \[ \mathrm{d} S(t) = r S(t) \mathrm{d} t + \sigma S(t)\mathrm{d}W(t) \]
static double	huntKennedyCMSAdjustedRate(double forwardSwaprate, double volatility, double swapAnnuity, double optionMaturity, double swapMaturity, double payoffUnit) Calculate the adjusted forward swaprate corresponding to a change of payoff unit from the given swapAnnuity to the given payoffUnit using the Black-Scholes model for the swap rate together with the Hunt-Kennedy convexity adjustment.
static double	huntKennedyCMSFloorValue(double forwardSwaprate, double volatility, double swapAnnuity, double optionMaturity, double swapMaturity, double payoffUnit, double optionStrike) Calculate the value of a CMS strike using the Black-Scholes model for the swap rate together with the Hunt-Kennedy convexity adjustment.
static double	huntKennedyCMSOptionValue(double forwardSwaprate, double volatility, double swapAnnuity, double optionMaturity, double swapMaturity, double payoffUnit, double optionStrike) Calculate the value of a CMS option using the Black-Scholes model for the swap rate together with the Hunt-Kennedy convexity adjustment.
static double	margrabeExchangeOptionValue(double spot1, double spot2, double volatility1, double volatility2, double correlation, double optionMaturity) Calculates the value of an Exchange option under a generalized Black-Scholes model, i.e., the payoff \( max(S_{1}(T)-S_{2}(T),0) \), where \( S_{1} \) and \( S_{2} \) follow a log-normal process with constant log-volatility and constant instantaneous correlation.
static double	price(Date settlementDate, Date maturityDate, double coupon, double yield, double redemption, int frequency) Re-implementation of the Excel PRICE function (a rather primitive bond price formula).
static double	price(double timeToMaturity, double coupon, double yield, double redemption, int frequency) Re-implementation of the Excel PRICE function (a rather primitive bond price formula).
static double	sabrAlphaApproximation(double normalVolatility, double beta, double rho, double nu, double displacement, double underlying, double maturity) Return the parameter alpha (initial value of the stochastic vol process) of a SABR model using the to match the given at-the-money volatility.
static double	sabrBerestyckiNormalVolatilityApproximation(double alpha, double beta, double rho, double nu, double displacement, double underlying, double strike, double maturity) Return the implied normal volatility (Bachelier volatility) under a SABR model using the approximation of Berestycki.
static double	sabrHaganLognormalBlackVolatilityApproximation(double alpha, double beta, double rho, double nu, double underlying, double strike, double maturity) Calculated the approximation to the lognormal Black volatility using the standard SABR model and the standard Hagan approximation.
static double	sabrHaganLognormalBlackVolatilityApproximation(double alpha, double beta, double rho, double nu, double displacement, double underlying, double strike, double maturity) Calculated the approximation to the lognormal Black volatility using the standard SABR model and the standard Hagan approximation.
static double	sabrNormalVolatilityApproximation(double alpha, double beta, double rho, double nu, double displacement, double underlying, double strike, double maturity) Return the implied normal volatility (Bachelier volatility) under a SABR model using the approximation of Hagan.
static double	sabrNormalVolatilityCurvatureApproximation(double alpha, double beta, double rho, double nu, double displacement, double underlying, double maturity) Return the curvature of the implied normal volatility (Bachelier volatility) under a SABR model using the approximation of Berestycki.
static double	sabrNormalVolatilitySkewApproximation(double alpha, double beta, double rho, double nu, double displacement, double underlying, double maturity) Return the skew of the implied normal volatility (Bachelier volatility) under a SABR model using the approximation of Berestycki.
static double	volatilityConversionLognormalATMtoNormalATM(double forward, double displacement, double maturity, double lognormalVolatiltiy) Exact conversion of displaced lognormal ATM volatiltiy to normal ATM volatility.

Methods inherited from class java.lang.Object

clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Method Detail

blackScholesGeneralizedOptionValue

public static double blackScholesGeneralizedOptionValue(double forward,
                                                        double volatility,
                                                        double optionMaturity,
                                                        double optionStrike,
                                                        double payoffUnit)

Calculates the Black-Scholes option value of a call, i.e., the payoff max(S(T)-K,0) P, where S follows a log-normal process with constant log-volatility. The method also handles cases where the forward and/or option strike is negative and some limit cases where the forward and/or the option strike is zero.

Parameters:

forward - The forward of the underlying.

volatility - The Black-Scholes volatility.

optionMaturity - The option maturity T.

optionStrike - The option strike. If the option strike is ≤ 0.0 the method returns the value of the forward contract paying S(T)-K in T.

payoffUnit - The payoff unit (e.g., the discount factor)

Returns:

Returns the value of a European call option under the Black-Scholes model.

blackScholesGeneralizedOptionValue

public static RandomVariable blackScholesGeneralizedOptionValue(RandomVariable forward,
                                                                RandomVariable volatility,
                                                                double optionMaturity,
                                                                double optionStrike,
                                                                RandomVariable payoffUnit)

Calculates the Black-Scholes option value of a call, i.e., the payoff max(S(T)-K,0) P, where S follows a log-normal process with constant log-volatility. The model specific quantities are considered to be random variable, i.e., the function may calculate an per-path valuation in a single call.

Parameters:

forward - The forward of the underlying.

volatility - The Black-Scholes volatility.

optionMaturity - The option maturity T.

optionStrike - The option strike. If the option strike is ≤ 0.0 the method returns the value of the forward contract paying S(T)-K in T.

payoffUnit - The payoff unit (e.g., the discount factor)

Returns:

Returns the value of a European call option under the Black-Scholes model.

blackScholesOptionValue

public static double blackScholesOptionValue(double initialStockValue,
                                             double riskFreeRate,
                                             double volatility,
                                             double optionMaturity,
                                             double optionStrike)

Calculates the Black-Scholes option value of a call, i.e., the payoff max(S(T)-K,0), where S follows a log-normal process with constant log-volatility.

Parameters:

initialStockValue - The spot value of the underlying.

riskFreeRate - The risk free rate r (df = exp(-r T)).

volatility - The Black-Scholes volatility.

optionMaturity - The option maturity T.

optionStrike - The option strike. If the option strike is ≤ 0.0 the method returns the value of the forward contract paying S(T)-K in T.

Returns:

Returns the value of a European call option under the Black-Scholes model.

blackScholesOptionValue

public static double blackScholesOptionValue(double initialStockValue,
                                             double riskFreeRate,
                                             double volatility,
                                             double optionMaturity,
                                             double optionStrike,
                                             boolean isCall)

Calculates the Black-Scholes option value of a call, i.e., the payoff max(S(T)-K,0), or a put, i.e., the payoff max(K-S(T),0), where S follows a log-normal process with constant log-volatility.

Parameters:

initialStockValue - The spot value of the underlying.

riskFreeRate - The risk free rate r (df = exp(-r T)).

volatility - The Black-Scholes volatility.

optionMaturity - The option maturity T.

optionStrike - The option strike. If the option strike is ≤ 0.0 the method returns the value of the forward contract paying S(T)-K in T for the call and zero for the put.

isCall - If true, the value of a call is calculated, if false, the value of a put is calculated.

Returns:

Returns the value of a European call/put option under the Black-Scholes model.

blackScholesATMOptionValue

public static double blackScholesATMOptionValue(double volatility,
                                                double optionMaturity,
                                                double forward,
                                                double payoffUnit)

Calculates the Black-Scholes option value of an atm call option.

Parameters:

volatility - The Black-Scholes volatility.

optionMaturity - The option maturity T.

forward - The forward, i.e., the expectation of the index under the measure associated with payoff unit.

payoffUnit - The payoff unit, i.e., the discount factor or the anuity associated with the payoff.

Returns:

Returns the value of a European at-the-money call option under the Black-Scholes model

blackScholesOptionDelta

public static double blackScholesOptionDelta(double initialStockValue,
                                             double riskFreeRate,
                                             double volatility,
                                             double optionMaturity,
                                             double optionStrike)

Calculates the delta of a call option under a Black-Scholes model The method also handles cases where the forward and/or option strike is negative and some limit cases where the forward or the option strike is zero. In the case forward = option strike = 0 the method returns 1.0.

Parameters:

initialStockValue - The initial value of the underlying, i.e., the spot.

riskFreeRate - The risk free rate of the bank account numerarie.

volatility - The Black-Scholes volatility.

optionMaturity - The option maturity T.

optionStrike - The option strike.

Returns:

The delta of the option

blackScholesOptionDelta

public static RandomVariable blackScholesOptionDelta(RandomVariable initialStockValue,
                                                     RandomVariable riskFreeRate,
                                                     RandomVariable volatility,
                                                     double optionMaturity,
                                                     double optionStrike)

Calculates the delta of a call option under a Black-Scholes model The method also handles cases where the forward and/or option strike is negative and some limit cases where the forward or the option strike is zero. In the case forward = option strike = 0 the method returns 1.0.

Parameters:

initialStockValue - The initial value of the underlying, i.e., the spot.

riskFreeRate - The risk free rate of the bank account numerarie.

volatility - The Black-Scholes volatility.

optionMaturity - The option maturity T.

optionStrike - The option strike.

Returns:

The delta of the option

blackScholesOptionDelta

public static RandomVariable blackScholesOptionDelta(RandomVariable initialStockValue,
                                                     RandomVariable riskFreeRate,
                                                     RandomVariable volatility,
                                                     double optionMaturity,
                                                     RandomVariable optionStrike)

Calculates the delta of a call option under a Black-Scholes model The method also handles cases where the forward and/or option strike is negative and some limit cases where the forward or the option strike is zero. In the case forward = option strike = 0 the method returns 1.0.

Parameters:

initialStockValue - The initial value of the underlying, i.e., the spot.

riskFreeRate - The risk free rate of the bank account numerarie.

volatility - The Black-Scholes volatility.

optionMaturity - The option maturity T.

optionStrike - The option strike.

Returns:

The delta of the option

blackScholesOptionGamma

public static double blackScholesOptionGamma(double initialStockValue,
                                             double riskFreeRate,
                                             double volatility,
                                             double optionMaturity,
                                             double optionStrike)

This static method calculated the gamma of a call option under a Black-Scholes model

Parameters:

initialStockValue - The initial value of the underlying, i.e., the spot.

riskFreeRate - The risk free rate of the bank account numerarie.

volatility - The Black-Scholes volatility.

optionMaturity - The option maturity T.

optionStrike - The option strike.

Returns:

The gamma of the option

blackScholesOptionGamma

public static RandomVariable blackScholesOptionGamma(RandomVariable initialStockValue,
                                                     RandomVariable riskFreeRate,
                                                     RandomVariable volatility,
                                                     double optionMaturity,
                                                     double optionStrike)

This static method calculated the gamma of a call option under a Black-Scholes model

Parameters:

initialStockValue - The initial value of the underlying, i.e., the spot.

riskFreeRate - The risk free rate of the bank account numerarie.

volatility - The Black-Scholes volatility.

optionMaturity - The option maturity T.

optionStrike - The option strike.

Returns:

The gamma of the option

blackScholesOptionVega

public static double blackScholesOptionVega(double initialStockValue,
                                            double riskFreeRate,
                                            double volatility,
                                            double optionMaturity,
                                            double optionStrike)

Calculates the vega of a call, i.e., the payoff max(S(T)-K,0) P, where S follows a normal process with constant volatility, i.e., a Black-Scholes model \[ \mathrm{d} S(t) = r S(t) \mathrm{d} t + \sigma S(t)\mathrm{d}W(t) \]

Parameters:

initialStockValue - The initial value of the underlying, i.e., the spot.

riskFreeRate - The risk free rate of the bank account numerarie.

volatility - The Black-Scholes volatility.

optionMaturity - The option maturity T.

optionStrike - The option strike.

Returns:

The vega of the option

blackScholesOptionTheta

public static double blackScholesOptionTheta(double initialStockValue,
                                             double riskFreeRate,
                                             double volatility,
                                             double optionMaturity,
                                             double optionStrike)

This static method calculated the vega of a call option under a Black-Scholes model

Parameters:

initialStockValue - The initial value of the underlying, i.e., the spot.

riskFreeRate - The risk free rate of the bank account numerarie.

volatility - The Black-Scholes volatility.

optionMaturity - The option maturity T.

optionStrike - The option strike.

Returns:

The vega of the option

blackScholesOptionRho

public static double blackScholesOptionRho(double initialStockValue,
                                           double riskFreeRate,
                                           double volatility,
                                           double optionMaturity,
                                           double optionStrike)

This static method calculated the rho of a call option under a Black-Scholes model

Parameters:

initialStockValue - The initial value of the underlying, i.e., the spot.

riskFreeRate - The risk free rate of the bank account numerarie.

volatility - The Black-Scholes volatility.

optionMaturity - The option maturity T.

optionStrike - The option strike.

Returns:

The rho of the option

blackScholesOptionImpliedVolatility

public static double blackScholesOptionImpliedVolatility(double forward,
                                                         double optionMaturity,
                                                         double optionStrike,
                                                         double payoffUnit,
                                                         double optionValue)

Calculates the Black-Scholes option implied volatility of a call, i.e., the payoff

max(S(T)-K,0)

, where S follows a log-normal process with constant log-volatility. The admissible values for optionValue are between forward * payoffUnit - optionStrike (the inner value) and forward * payoffUnit.

Parameters:

forward - The forward of the underlying (which is equal to S(0) / payoffUnit, given the spot value S(0)).

optionMaturity - The option maturity T.

optionStrike - The option strike. If the option strike is ≤ 0.0 the method returns the value of the forward contract paying S(T)-K in T.

payoffUnit - The payoff unit (e.g., the discount factor), (which is equal to exp(-maturity * r), given the interest rate r).

optionValue - The option value. The admissible values for optionValue are between forward * payoffUnit - optionStrike (the inner value) and forward * payoffUnit.

Returns:

Returns the implied volatility of a European call option under the Black-Scholes model.

blackScholesDigitalOptionValue

public static double blackScholesDigitalOptionValue(double initialStockValue,
                                                    double riskFreeRate,
                                                    double volatility,
                                                    double optionMaturity,
                                                    double optionStrike)

Calculates the Black-Scholes option value of a digital call option.

Parameters:

initialStockValue - The initial value of the underlying, i.e., the spot.

riskFreeRate - The risk free rate of the bank account numerarie.

volatility - The Black-Scholes volatility.

optionMaturity - The option maturity T.

optionStrike - The option strike.

Returns:

Returns the value of a European call option under the Black-Scholes model

blackScholesDigitalOptionDelta

public static double blackScholesDigitalOptionDelta(double initialStockValue,
                                                    double riskFreeRate,
                                                    double volatility,
                                                    double optionMaturity,
                                                    double optionStrike)

Calculates the delta of a digital option under a Black-Scholes model

Parameters:

initialStockValue - The initial value of the underlying, i.e., the spot.

riskFreeRate - The risk free rate of the bank account numerarie.

volatility - The Black-Scholes volatility.

optionMaturity - The option maturity T.

optionStrike - The option strike.

Returns:

The delta of the digital option

blackScholesDigitalOptionVega

public static double blackScholesDigitalOptionVega(double initialStockValue,
                                                   double riskFreeRate,
                                                   double volatility,
                                                   double optionMaturity,
                                                   double optionStrike)

Calculates the vega of a digital option under a Black-Scholes model

Parameters:

initialStockValue - The initial value of the underlying, i.e., the spot.

riskFreeRate - The risk free rate of the bank account numerarie.

volatility - The Black-Scholes volatility.

optionMaturity - The option maturity T.

optionStrike - The option strike.

Returns:

The vega of the digital option

blackScholesDigitalOptionRho

public static double blackScholesDigitalOptionRho(double initialStockValue,
                                                  double riskFreeRate,
                                                  double volatility,
                                                  double optionMaturity,
                                                  double optionStrike)

Calculates the rho of a digital option under a Black-Scholes model

Parameters:

initialStockValue - The initial value of the underlying, i.e., the spot.

riskFreeRate - The risk free rate of the bank account numerarie.

volatility - The Black-Scholes volatility.

optionMaturity - The option maturity T.

optionStrike - The option strike.

Returns:

The rho of the digital option

blackModelCapletValue

public static double blackModelCapletValue(double forward,
                                           double volatility,
                                           double optionMaturity,
                                           double optionStrike,
                                           double periodLength,
                                           double discountFactor)

Calculate the value of a caplet assuming the Black'76 model.

Parameters:

forward - The forward (spot).

volatility - The Black'76 volatility.

optionMaturity - The option maturity

optionStrike - The option strike.

periodLength - The period length of the underlying forward rate.

discountFactor - The discount factor corresponding to the payment date (option maturity + period length).

Returns:

Returns the value of a caplet under the Black'76 model

blackModelDgitialCapletValue

public static double blackModelDgitialCapletValue(double forward,
                                                  double volatility,
                                                  double periodLength,
                                                  double discountFactor,
                                                  double optionMaturity,
                                                  double optionStrike)

Calculate the value of a digital caplet assuming the Black'76 model.

Parameters:

forward - The forward (spot).

volatility - The Black'76 volatility.

periodLength - The period length of the underlying forward rate.

discountFactor - The discount factor corresponding to the payment date (option maturity + period length).

optionMaturity - The option maturity

optionStrike - The option strike.

Returns:

Returns the price of a digital caplet under the Black'76 model

blackModelSwaptionValue

public static double blackModelSwaptionValue(double forwardSwaprate,
                                             double volatility,
                                             double optionMaturity,
                                             double optionStrike,
                                             double swapAnnuity)

Calculate the value of a swaption assuming the Black'76 model.

Parameters:

forwardSwaprate - The forward (spot)

volatility - The Black'76 volatility.

optionMaturity - The option maturity.

optionStrike - The option strike.

swapAnnuity - The swap annuity corresponding to the underlying swap.

Returns:

Returns the value of a Swaption under the Black'76 model

margrabeExchangeOptionValue

public static double margrabeExchangeOptionValue(double spot1,
                                                 double spot2,
                                                 double volatility1,
                                                 double volatility2,
                                                 double correlation,
                                                 double optionMaturity)

Calculates the value of an Exchange option under a generalized Black-Scholes model, i.e., the payoff \( max(S_{1}(T)-S_{2}(T),0) \), where \( S_{1} \) and \( S_{2} \) follow a log-normal process with constant log-volatility and constant instantaneous correlation. The method also handles cases where the forward and/or option strike is negative and some limit cases where the forward and/or the option strike is zero.

Parameters:

spot1 - Value of \( S_{1}(0) \)

spot2 - Value of \( S_{2}(0) \)

volatility1 - Volatility of \( \log(S_{1}(t)) \)

volatility2 - Volatility of \( \log(S_{2}(t)) \)

correlation - Instantaneous correlation of \( \log(S_{1}(t)) \) and \( \log(S_{2}(t)) \)

optionMaturity - The option maturity \( T \).

Returns:

Returns the value of a European exchange option under the Black-Scholes model.

bachelierOptionValue

public static double bachelierOptionValue(double forward,
                                          double volatility,
                                          double optionMaturity,
                                          double optionStrike,
                                          double payoffUnit)

Calculates the option value of a call, i.e., the payoff max(S(T)-K,0), where S follows a normal process with constant volatility, i.e., a Bachelier model \[ \mathrm{d} S(t) = r S(t) \mathrm{d} t + \sigma \mathrm{d}W(t) \]

Parameters:

forward - The forward of the underlying \( F = S(T) \exp(r T) \).

volatility - The Bachelier volatility \( \sigma \).

optionMaturity - The option maturity T.

optionStrike - The option strike K.

payoffUnit - The payoff unit (e.g., the discount factor)

Returns:

Returns the value of a European call option under the Bachelier model.

bachelierOptionValue

public static RandomVariable bachelierOptionValue(RandomVariable forward,
                                                  RandomVariable volatility,
                                                  double optionMaturity,
                                                  double optionStrike,
                                                  RandomVariable payoffUnit)

Calculates the option value of a call, i.e., the payoff max(S(T)-K,0) P, where S follows a normal process with constant volatility, i.e., a Bachelier model \[ \mathrm{d} S(t) = r S(t) \mathrm{d} t + \sigma \mathrm{d}W(t) \]

Parameters:

forward - The forward of the underlying \( F = S(T) \exp(r T) \).

volatility - The Bachelier volatility \( \sigma \).

optionMaturity - The option maturity T.

optionStrike - The option strike.

payoffUnit - The payoff unit (e.g., the discount factor)

Returns:

Returns the value of a European call option under the Bachelier model.

bachelierOptionImpliedVolatility

public static double bachelierOptionImpliedVolatility(double forward,
                                                      double optionMaturity,
                                                      double optionStrike,
                                                      double payoffUnit,
                                                      double optionValue)

Calculates the Bachelier option implied volatility of a call, i.e., the payoff

max(S(T)-K,0)

, where S follows a normal process with constant volatility.

Parameters:

forward - The forward of the underlying.

optionMaturity - The option maturity T.

optionStrike - The option strike. If the option strike is ≤ 0.0 the method returns the value of the forward contract paying S(T)-K in T.

payoffUnit - The payoff unit (e.g., the discount factor)

optionValue - The option value.

Returns:

Returns the implied volatility of a European call option under the Bachelier model.

bachelierOptionDelta

public static double bachelierOptionDelta(double forward,
                                          double volatility,
                                          double optionMaturity,
                                          double optionStrike,
                                          double payoffUnit)

Calculates the option delta of a call, i.e., the payoff max(S(T)-K,0), where S follows a normal process with constant volatility, i.e., a Bachelier model \[ \mathrm{d} S(t) = r S(t) \mathrm{d} t + \sigma \mathrm{d}W(t) \]

Parameters:

forward - The forward of the underlying \( F = S(T) \exp(r T) \).

volatility - The Bachelier volatility \( \sigma \).

optionMaturity - The option maturity T.

optionStrike - The option strike K.

payoffUnit - The payoff unit (e.g., the discount factor)

Returns:

Returns the value of a European call option under the Bachelier model.

huntKennedyCMSOptionValue

public static double huntKennedyCMSOptionValue(double forwardSwaprate,
                                               double volatility,
                                               double swapAnnuity,
                                               double optionMaturity,
                                               double swapMaturity,
                                               double payoffUnit,
                                               double optionStrike)

Calculate the value of a CMS option using the Black-Scholes model for the swap rate together with the Hunt-Kennedy convexity adjustment.

Parameters:

forwardSwaprate - The forward swap rate

volatility - Volatility of the log of the swap rate

swapAnnuity - The swap annuity

optionMaturity - The option maturity

swapMaturity - The swap maturity

payoffUnit - The payoff unit, e.g., the discount factor corresponding to the payment date

optionStrike - The option strike

Returns:

Value of the CMS option

huntKennedyCMSFloorValue

public static double huntKennedyCMSFloorValue(double forwardSwaprate,
                                              double volatility,
                                              double swapAnnuity,
                                              double optionMaturity,
                                              double swapMaturity,
                                              double payoffUnit,
                                              double optionStrike)

Calculate the value of a CMS strike using the Black-Scholes model for the swap rate together with the Hunt-Kennedy convexity adjustment.

Parameters:

forwardSwaprate - The forward swap rate

volatility - Volatility of the log of the swap rate

swapAnnuity - The swap annuity

optionMaturity - The option maturity

swapMaturity - The swap maturity

payoffUnit - The payoff unit, e.g., the discount factor corresponding to the payment date

optionStrike - The option strike

Returns:

Value of the CMS strike

huntKennedyCMSAdjustedRate

public static double huntKennedyCMSAdjustedRate(double forwardSwaprate,
                                                double volatility,
                                                double swapAnnuity,
                                                double optionMaturity,
                                                double swapMaturity,
                                                double payoffUnit)

Calculate the adjusted forward swaprate corresponding to a change of payoff unit from the given swapAnnuity to the given payoffUnit using the Black-Scholes model for the swap rate together with the Hunt-Kennedy convexity adjustment.

Parameters:

forwardSwaprate - The forward swap rate

volatility - Volatility of the log of the swap rate

swapAnnuity - The swap annuity

optionMaturity - The option maturity

swapMaturity - The swap maturity

payoffUnit - The payoff unit, e.g., the discount factor corresponding to the payment date

Returns:

Convexity adjusted forward rate

sabrHaganLognormalBlackVolatilityApproximation

public static double sabrHaganLognormalBlackVolatilityApproximation(double alpha,
                                                                    double beta,
                                                                    double rho,
                                                                    double nu,
                                                                    double underlying,
                                                                    double strike,
                                                                    double maturity)

Calculated the approximation to the lognormal Black volatility using the standard SABR model and the standard Hagan approximation.

Parameters:

alpha - initial value of the stochastic volatility process of the SABR model.

beta - CEV parameter of the SABR model.

rho - Correlation (leverages) of the stochastic volatility.

nu - Volatility of the stochastic volatility (vol-of-vol).

underlying - Underlying (spot) value.

strike - Strike.

maturity - Maturity.

Returns:

Implied lognormal Black volatility.

sabrHaganLognormalBlackVolatilityApproximation

public static double sabrHaganLognormalBlackVolatilityApproximation(double alpha,
                                                                    double beta,
                                                                    double rho,
                                                                    double nu,
                                                                    double displacement,
                                                                    double underlying,
                                                                    double strike,
                                                                    double maturity)

Calculated the approximation to the lognormal Black volatility using the standard SABR model and the standard Hagan approximation.

Parameters:

alpha - initial value of the stochastic volatility process of the SABR model.

beta - CEV parameter of the SABR model.

rho - Correlation (leverages) of the stochastic volatility.

nu - Volatility of the stochastic volatility (vol-of-vol).

displacement - The displacement parameter d.

underlying - Underlying (spot) value.

strike - Strike.

maturity - Maturity.

Returns:

Implied lognormal Black volatility.

sabrBerestyckiNormalVolatilityApproximation

public static double sabrBerestyckiNormalVolatilityApproximation(double alpha,
                                                                 double beta,
                                                                 double rho,
                                                                 double nu,
                                                                 double displacement,
                                                                 double underlying,
                                                                 double strike,
                                                                 double maturity)

Return the implied normal volatility (Bachelier volatility) under a SABR model using the approximation of Berestycki.

Parameters:

alpha - initial value of the stochastic volatility process of the SABR model.

beta - CEV parameter of the SABR model.

rho - Correlation (leverages) of the stochastic volatility.

nu - Volatility of the stochastic volatility (vol-of-vol).

displacement - The displacement parameter d.

underlying - Underlying (spot) value.

strike - Strike.

maturity - Maturity.

Returns:

The implied normal volatility (Bachelier volatility)

sabrNormalVolatilityApproximation

public static double sabrNormalVolatilityApproximation(double alpha,
                                                       double beta,
                                                       double rho,
                                                       double nu,
                                                       double displacement,
                                                       double underlying,
                                                       double strike,
                                                       double maturity)

Return the implied normal volatility (Bachelier volatility) under a SABR model using the approximation of Hagan.

Parameters:

alpha - initial value of the stochastic volatility process of the SABR model.

beta - CEV parameter of the SABR model.

rho - Correlation (leverages) of the stochastic volatility.

nu - Volatility of the stochastic volatility (vol-of-vol).

displacement - The displacement parameter d.

underlying - Underlying (spot) value.

strike - Strike.

maturity - Maturity.

Returns:

The implied normal volatility (Bachelier volatility)

sabrAlphaApproximation

public static double sabrAlphaApproximation(double normalVolatility,
                                            double beta,
                                            double rho,
                                            double nu,
                                            double displacement,
                                            double underlying,
                                            double maturity)

Return the parameter alpha (initial value of the stochastic vol process) of a SABR model using the to match the given at-the-money volatility.

Parameters:

normalVolatility - ATM volatility to match.

beta - CEV parameter of the SABR model.

rho - Correlation (leverages) of the stochastic volatility.

nu - Volatility of the stochastic volatility (vol-of-vol).

displacement - The displacement parameter d.

underlying - Underlying (spot) value.

maturity - Maturity.

Returns:

The implied normal volatility (Bachelier volatility)

sabrNormalVolatilitySkewApproximation

public static double sabrNormalVolatilitySkewApproximation(double alpha,
                                                           double beta,
                                                           double rho,
                                                           double nu,
                                                           double displacement,
                                                           double underlying,
                                                           double maturity)

Return the skew of the implied normal volatility (Bachelier volatility) under a SABR model using the approximation of Berestycki. The skew is the first derivative of the implied vol w.r.t. the strike, evaluated at the money.

Parameters:

alpha - initial value of the stochastic volatility process of the SABR model.

beta - CEV parameter of the SABR model.

rho - Correlation (leverages) of the stochastic volatility.

nu - Volatility of the stochastic volatility (vol-of-vol).

displacement - The displacement parameter d.

underlying - Underlying (spot) value.

maturity - Maturity.

Returns:

The skew of the implied normal volatility (Bachelier volatility)

sabrNormalVolatilityCurvatureApproximation

public static double sabrNormalVolatilityCurvatureApproximation(double alpha,
                                                                double beta,
                                                                double rho,
                                                                double nu,
                                                                double displacement,
                                                                double underlying,
                                                                double maturity)

Return the curvature of the implied normal volatility (Bachelier volatility) under a SABR model using the approximation of Berestycki. The curvatures is the second derivative of the implied vol w.r.t. the strike, evaluated at the money.

Parameters:

alpha - initial value of the stochastic volatility process of the SABR model.

beta - CEV parameter of the SABR model.

rho - Correlation (leverages) of the stochastic volatility.

nu - Volatility of the stochastic volatility (vol-of-vol).

displacement - The displacement parameter d.

underlying - Underlying (spot) value.

maturity - Maturity.

Returns:

The curvature of the implied normal volatility (Bachelier volatility)

volatilityConversionLognormalATMtoNormalATM

public static double volatilityConversionLognormalATMtoNormalATM(double forward,
                                                                 double displacement,
                                                                 double maturity,
                                                                 double lognormalVolatiltiy)

Exact conversion of displaced lognormal ATM volatiltiy to normal ATM volatility.

Parameters:

forward - The forward

displacement - The displacement (considering a displaced lognormal model, otherwise 0.

maturity - The maturity

lognormalVolatiltiy - The (implied) lognormal volatility.

Returns:

The (implied) normal volatility.

See Also:

Dimitroff, Fries, Lichtner and Rodi: Lognormal vs Normal Volatilities and Sensitivities in Practice

price

public static double price(Date settlementDate,
                           Date maturityDate,
                           double coupon,
                           double yield,
                           double redemption,
                           int frequency)

Re-implementation of the Excel PRICE function (a rather primitive bond price formula). The re-implementation is not exact, because this function does not consider daycount conventions.

Parameters:

settlementDate - Valuation date.

maturityDate - Maturity date of the bond.

coupon - Coupon payment.

yield - Yield (discount factor, using frequency: 1/(1 + yield/frequency).

redemption - Redemption (notional repayment).

frequency - Frequency (1,2,4).

Returns:

price Clean price.

price

public static double price(double timeToMaturity,
                           double coupon,
                           double yield,
                           double redemption,
                           int frequency)

Re-implementation of the Excel PRICE function (a rather primitive bond price formula). The re-implementation is not exact, because this function does not consider daycount conventions. We assume we have (int)timeToMaturity/frequency future periods and the running period has an accrual period of timeToMaturity - frequency * ((int)timeToMaturity/frequency).

Parameters:

timeToMaturity - The time to maturity.

coupon - Coupon payment.

yield - Yield (discount factor, using frequency: 1/(1 + yield/frequency).

redemption - Redemption (notional repayment).

frequency - Frequency (1,2,4).

Returns:

price Clean price.

bachelierGeneralizedOptionVega

public static double bachelierGeneralizedOptionVega(double forward,
                                                    double volatility,
                                                    double optionMaturity,
                                                    double optionStrike,
                                                    double payoffUnit)

Calculates the vega of a call, i.e., the payoff max(S(T)-K,0) P, where S follows a normal process with constant volatility, i.e., a Bachelier model \[ \mathrm{d} S(t) = r S(t) \mathrm{d} t + \sigma \mathrm{d}W(t) \]

Parameters:

forward - The forward of the underlying \( F = S(T) \exp(r T) \).

volatility - The Bachelier volatility \( \sigma \).

optionMaturity - The option maturity T.

optionStrike - The option strike.

payoffUnit - The payoff unit (e.g., the discount factor)

Returns:

Returns the vega of a European call option under the Bachelier model.

blackScholesGeneralizedOptionVega

public static double blackScholesGeneralizedOptionVega(double forward,
                                                       double volatility,
                                                       double optionMaturity,
                                                       double optionStrike,
                                                       double payoffUnit)

Calculates the vega of a call, i.e., the payoff max(S(T)-K,0) P, where S follows a normal process with constant volatility, i.e., a Black-Scholes model \[ \mathrm{d} S(t) = r S(t) \mathrm{d} t + \sigma S(t)\mathrm{d}W(t) \]

Parameters:

forward - The forward of the underlying \( F = S(T) \exp(r T) \).

volatility - The Black-Scholes volatility \( \sigma \).

optionMaturity - The option maturity T.

optionStrike - The option strike.

payoffUnit - The payoff unit (e.g., the discount factor)

Returns:

Returns the vega of a European call option under the Black-Scholes model.