Source code for frds.measures._option_price

from typing import Tuple
import numpy as np
from scipy.stats import norm




def blsprice(
    S: float, K: float, r: float, T: float, sigma: float, q=0.0
) -> Tuple[float, float]:
    """European option values from Black-Scholes model.
    See :doc:`/measures/option_prices`.

    Args:
        S (float): Current price of the underlying asset.
        K (float): Strike (exercise) price of the option.
        r (float): Annualized continuously compounded risk-free rate of return over the life of the option, expressed as a positive decimal number
        T (float): Time to expiration of the option, expressed in years.
        sigma (float): Annualized asset price volatility (i.e., annualized standard deviation of the continuously compounded asset return), expressed as a positive decimal number.
        q (float, optional): Annualized continuously compounded yield of the underlying asset over the life of the option, expressed as a decimal number. Defaults to 0.0.

    Returns:
        Tuple[float, float]: Prices of European call and put options
    """
    d1 = (np.log(S / K) + (r - q + sigma**2 / 2) * T) / (sigma * np.sqrt(T))
    d2 = d1 - sigma * np.sqrt(T)
    call = S * np.exp(-q * T) * norm.cdf(d1) - K * np.exp(-r * T) * norm.cdf(d2)
    put = K * np.exp(-r * T) * norm.cdf(-d2) - S * np.exp(-q * T) * norm.cdf(-d1)
    return call, put