Code covered by the BSD License  

Highlights from
Exercises in Advanced Risk and Portfolio Management

from Exercises in Advanced Risk and Portfolio Management by Attilio Meucci
text and comments on solutions available at http://symmys.com/node/170

NormalCopulaPDF(u,Mu,Sigma)
function F_U = NormalCopulaPDF(u,Mu,Sigma)

N=length(u);
s=sqrt(diag(Sigma));

x=norminv(u,Mu,s);

Numerator = (2*pi)^(-N/2) * (  (det(Sigma))^(-.5) ) * exp(-.5*(x-Mu)'*inv(Sigma)*(x-Mu));

fs=normpdf(x,Mu,s);
Denominator = prod(fs);

F_U = Numerator/Denominator;


Contact us at files@mathworks.com