|
|
楼主

楼主 |
发表于 2008-5-30 21:49:28
|
只看该作者
请教一个关于IML的问题,大家帮我看看我的程序哪里错了
刚刚接触IML,实在是摸不着头脑,涉及到模块的嵌套,我就是想求数值积分,那个fun(k)里又有一个变量i,对于不同的i值求fun(k)的数值积分,可能有好多错误,请大家指正,谢谢!
proc iml;
start quad(i);
start fun(k);
ds=cos(k*i)*exp(-(1.1**1.2)*k**1.2*cos(1.2*atan(80*k))/cos(0.6*3.1415926535)+(1.1)**1.2*(1/80)**1.2/cos(0.6*3.1415926535));
return(ds);
finish;
a={0 .p};
call quad(z,'fun',a) peak=0;
finish quad;
T=quad({0,0.1,0.2,0.3});
print T; |
|