image thumbnail
from PlotMeTheGreeks by Yazann Romahi
VisualizationTool for Option Greeks in 3 and 4 dimensions.

VarVega=bsVarVega(CallPutFlag, S,X,T,r,r_f, sigma)
function VarVega=bsVarVega(CallPutFlag, S,X,T,r,r_f, sigma) 
    %sigma=sigma+0.0000000001;
    %T=T+0.0000000001;
    d1=(log(S./X)+(r-r_f+(sigma.^2)./2).*T)./(sigma.*sqrt(T));
    nd1=(1./sqrt(2.*pi))*exp(-d1.^2/2);
    %vega=S.*exp(-r_f.*T).*nd1.*sqrt(T)
    VarVega=(S.*sqrt(T)./(2*sigma)) .* nd1;
end

% bsVarVega(CallPutFlag, S, X, T, r, r_f, sigma)
% Returns Garman-Kohlhagen Variance-Vega for a European currency option

%   version: v1.0
%   author: Yazann Romahi [yazann@romahi.com]

Contact us at files@mathworks.com