一、定积分

定积分的求法；

y = quad(@Fun,a,b)       Fun为函数，a为下限，b为上限

y = quad(@Fun,a,b，ε)    ε为误差，默认为10^(-6)

用quad求积分必须要定义函数，调用函数时，要用@

function [ y ] = Fun(x)
 y = x.^2;

end

>> y = quad(@Fun,-1,1)

y =

    0.6667

还有一种 quad1求积分，方法与quad相同

quad:采用递推自适应的Simpson法来计算

quad1:采用递推自适应的Lobatto法来计算