Can anyone help me Thanks
- To: mathgroup at smc.vnet.net
- Subject: [mg75695] Can anyone help me Thanks
- From: 佳玲chialing <chia20040825 at yahoo.com.hk>
- Date: Wed, 9 May 2007 04:17:03 -0400 (EDT)
My question is as follow,
'S1,s1,S2,s2 is constant, fn is a function with variable A1, A2
<< Statistics`MultinormalDistribution`;
rr = {{1, 0.75}, {0.75, 1}};
ndist2 = MultinormalDistribution[{0, 0}, rr];
pf2[A1_, A2_] : = {PDF[ndist2,{(Log[A1/S1]-(s1^2/2))/s1,((Log[A2/S2]-(s2^2/2))/s2}]}[[1]];
z = NIntegrate[1/(A1 A2)*(fn[A1, A2]pf2[A1, A2]), {A1, 0, Infinity}, {A2, 0, Infinity}, MinRecursion -> 5, MaxRecursion -> 10];
I want to extend the 36 dimensional MultinormalDistribution, could you help me?
i.e.
<< Statistics`MultinormalDistribution`;
36r = {{1, 0.75,........................}, {0.75, 1,..............}.........................};
ndist36 = MultinormalDistribution[{0, 0,...................}, 36r];
pf36[A1_, A2_,A3_...............] : = {PDF[ndist36,{(Log[A1/S1]-(s1^2/2))/s1,((Log[A2/S2]-(s2^2/2))/s2},..........]}[[1]];
z = NIntegrate[1/(A1 A2 .........)*(fn[A1, A2,...........]pf36[A1, A2,.........]), {A1, 0, Infinity}, {A2, 0, Infinity},................ MinRecursion -> 5, MaxRecursion -> 10];
Thank you again.
Sincerely,
Alice