function v = BS_CND(x)
   
syms t;
v = single(int(exp(- t ^ 2 / 2) / sqrt(2 * pi), -inf, x));