function [ v ] = BS_Gamma(S, X, T, v, r)

[d1, ] = BS_D1D2(S, X, T, v);

v = exp(-(d1 ^ 2) / 2) / sqrt(2 * pi) * exp(- r * T) / (S * v * sqrt(T));