function [ v ] = BS_TheoV(PC, S, X, T, v, r)
       
[d1, d2] = BS_D1D2(S, X, T, v);

v = S * BS_Delta(PC, S, X, T, v, r) - PC * X * exp(-r * T) * BS_CND(PC * d2);