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

d1 = (log(S / X) + (v ^ 2 / 2) * T) / (v * sqrt(T));
d2 = d1 - v * sqrt(T);