/*只能实现定积分*/
/*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;


雷达卡





京公网安备 11010802022788号







