SAS如何做三重积分？

webgu

1111

收藏 2013-04-04

悬赏 168 个论坛币已解决

例如：f(x, y, z)=b1*x +b2*y+b3*z;

在SAS可以做f(x,y,z)的三重积分吗？

如果可以，烦请给一个具体的实例。

谢谢。

最佳答案

ziyenano 查看完整内容

/*只能实现定积分*/ /*f(x, y, z)=b1*x +b2*y+b3*z*/ /*∫∫∫f(x,y,z)dxdydz 假设积分区间都定为[0,1]*/ proc iml; start fx(t) global(yv,zv,b1,b2,b3); f1=b1*t+b2*yv+b3*zv; return(f1); finish; start fy(v) global(yv,zv,b1,b2,b3,interval1); yv=v; call quad(f2,"fx",interval1); return(f2); finish; start fz(w) global(yv,zv,b1,b2,b3,interval2); zv=w; call quad(f3,"fy",interval2); return(f ...

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

ziyenano

2013-4-4 21:54:00

/*只能实现定积分*/
/*f(x, y, z)=b1*x +b2*y+b3*z*/
/*∫∫∫f(x,y,z)dxdydz 假设积分区间都定为[0,1]*/

proc iml;
start fx(t) global(yv,zv,b1,b2,b3);
f1=b1*t+b2*yv+b3*zv;
return(f1);
finish;

start fy(v) global(yv,zv,b1,b2,b3,interval1);
yv=v;
call quad(f2,"fx",interval1);
return(f2);
finish;

start fz(w) global(yv,zv,b1,b2,b3,interval2);
zv=w;
call quad(f3,"fy",interval2);
return(f3);
finish;

interval1={0 1};interval2={0 1};
interval3={0 1};
b1=1;b2=2;b3=3;

call quad(res,"fz",interval3);
print res;