SAS中文论坛

标题: 请教一个关于IML的问题，大家帮我看看我的程序哪里错了 [打印本页]

作者: shiyiming 时间: 2008-5-30 21:49
标题: 请教一个关于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;

欢迎光临 SAS中文论坛 (https://mysas.net/forum/)