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

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