scilab中求积分

scilab中求积分

===========================================================

• scilab中求一元积分函数：intg（a，b，f）

函数格式：[v,err]=intg(a,b,f [,ea [,er])

不过在scilab中求一元积分常用命令y＝intg（a，b，f）

这里的积分区间为（a，b）；

积分函数为 f；

例子：

deff('y=f(x);','y=x*x;'); //定义函数f（x）＝x*x

result=intg(1,2,f); //对函数f在（1，2）积分，

-->result //输入result，显示结果如下
result =

2.3333333

• 求二元积分函数：int2d（）

 

函数格式：[I,err]=int2d(X,Y,f [,params])

其中：

I：积分误差估计值；

X：一个3×n的矩阵，包含n个三角形顶点的X坐标信息；

Y:一个3×n的矩阵，包含n个三角形顶点的x坐标信息；

f：被积函数；

params：几分参数为实数向量[tol，iclose，maxtri，mevals，iflag]，

例子：

X=[0,0;1,1;1,0];//X，Y为两个三角形的坐标
Y=[0,0;0,1;1,1];
deff('z=f(x,y)','z=cos(x+y)')
[I,e]=int2d(X,Y,f)
// computes the integrand over the square [0 1]x[0 1]