davil2000 发表于 2012-12-17 14:13 
我的PC反应有些慢了
您好 能帮我走一下这个程序么 我想看一下得到的结果
毕业论文急用啊 但是安装几个版本的sas软件都安装不了 心急如焚啊
多谢多谢了 教授!
proc iml
cov={{1.0000 0.4015 0.4193
0.4015 1.0000 0.3639
0.4193 0.3639 1.0000}}
h={0.07,0.25,0.12};
T={[10000]0}
V={[10000]0}
r=0.03
value=0
rv=rannor(repeat(1000,10000,10))
Y=rv*(root(cov))
i=1
j=1
do while(i<=10000);
do while(j<=10)
x=cdf('Normal',Y[i,j])
Y[i,j])=-(log(1-x))/h[j]
j=j+1
end
T
=min(Y[i,])
if T<6 then v=50000*exp(-r*T)
value=value+v
j=1
i=i+1
end
value=value/10000
print value
run
quit