时光流逝
(1)被积函数的原函数不能用初等函数表示。
)2)被积函数很难用公式表示，用图形或表格给出。
应该建立定积分的近似计算方法、数值积分方法。
梯形法：
z=trapz(x，y ) ) )。
x=0:0.5:1； y=exp(-x.^2)； z=trapz(x，y ) z=0.7314
x=0:0.05:1； y=exp(-x.^2)； z=trapz(x，y ) z=0.7467
抛物线法： Simpson(XXDhh )公式
自适应xxdhh公式：
z=quad(f、a、b、tol ) )。
a为积分下限，b为积分上限，tol为计算精度，默认值为0.001。
调用quad函数时，首先创建描述被积函数的函数文件或语句函数。
以下称为内联函数和函数句柄
g=inline('exp )-x.^2) )； z=quad(g，- 1，1 ) z=1.4936 f=@(x ) x ) exp(-x.^2)； z=quad(f，- 1，1 ) z=1.4936
分别使用左端点矩形法、右端点的矩形法、梯形法
n=100； x=linspace (0，1，n )； y=4.//(1x.^2)； JXL=sum(y ) 1:(n-1 ) )/njxl=3.1202jxr=sum ) y ) 2:n )./njxr=3.1002 tx=trapz(x ) x，y ) tx=3.
有两个常用的求数值积分的函数：
z=quadl(f，a，b，tol )自适应复合Lobatto数值积分法
z=quadgk(f、a、c、tol )自适应复合Gauss-Kronrod数值积分法适用于高精度和振荡数值积分以及广义数值积分
g=inline('4./) 1x.^2) )； quad(g，0，1 ) ans=3.1416 quadl(g ) g，0，1 ) ans=3.1416
g=@(x )4./) 1x.^2)； quadGK(g，0，1 ) ans=3.1416
f1=@(x )1./(x.^2)； Z1=quad GK (f1，1，inf ) z1=1 f2=@(x ) x ) exp(-x.^2)； Z2=quadGK(F2，-inf，inf ) Z2=1.7725F3=@(x )1./(x.* sqrt )1-log(x ).^2)； z3=quad GK (F3，1，exp )1) z3=1.5708