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

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

v = PC * exp(- r * T) * BS_CND(PC * d1);