使用quadl函数可以进行Lobatto求积。该函数的语法格式如下。
q=quadl(fun,a,b)：求函数fun从a到b的积分近似，误差容限为1e-6。fun为M文件函数或匿名函数的句柄。函数y=fun(x)应该接受向量变量x并返回结果向量y。
q=quadl(fun,a,b,tol)：使用绝对误差容限tol代替默认误差容限1.0e-6。当tol的值更大时，需要的计算次数更少，计算更快，但计算精度差一些。
q=quadl(fun,a,b,tol,trace)：在递归过程中使用非0的trace参数。
[q,fcnt]=quadl(…)：返回函数计算次数。
【例23】用quadl函数求下面积分的数值积分
\[\int_{0}^{2} {\frac{1}{{x}^{3}-2x-5}dx}\]
首先编写函数的M文件myfun.m。
code.matlab
function y = myfun(x)
y = 1./(x.^3-2*x-5);然后在命令窗口中输入
code.matlab
>> Q = quadl(@myfun,0,2)
Q=
-0.4605