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

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