贝塞尔函数积分用matlab_MATLAB基础学习之数值积分（一次函数定积分）

数值积分的实现：

(一)比较quad与quadl函数的差别

基于自适应辛普森方法：[I,n]=quad(filename,a,b,tol,trace)

基于自适应Gauss-Lobatto方法：[I,n]=quadl(filename,a,b,tol,trace)

其中，filename是被积函数名；a和b分别是定积分的下限和上限，积分限[a,b]必须是有限的，不能为无穷大(inf)；tol用来控制积分精度，默认时取tol=10-6；trace控制是否展现积分过程，若取非0展现积分过程，取0则不展现，默认时取trace=0；返回参数I即定积分的值，n为被积函数的调用次数。

例：分别用quad函数和quadl函数求定积分的近似值，并在相同的积分精度下，比较被积函数的调用次数。

在命令行输入如下指令：

>> format long

>> f=@(x)4./(1+x.^2);

>> [I,n]=quad(f,0,1,1e-8)

I =

3.141592653733437

n =