【matlab】数值积分公式的程序实现

（一）专题实验（Newton-Cotes积分公式）

1、编写[a,b]上梯形积分公式、Simpson积分公式。

2、利用自己编写的程序计算定积分，计算一下数值解和精确解之间差的绝对值。

梯形积分：

function T=TX_int(f,a,b)

T=(b-a)/2*(f(a)+f(b));

>> TX_int(@(x)cos(x),0,pi/4)

ans =

0.6704

function T=TX_int(f,a,b)

T=(b-a)/2*(f(a)+f(b));

a=abs(sqrt(2)/2-T);

disp(a)

>> TX_int(@(x)cos(x),0,pi/4)

   0.0367

ans =

   0.6704

Simpson积分;

function T=TX_int(f,a,b)

T=(b-a)/6*(f(a)+4*f((a+b)/2)+f(b));

>> TX_int(@(x)cos(x),0,pi/4)

ans =

   0.7072

function T=TX_int(f,a,b)

T=(b-a)/6*(f(a)+4*f((a+b)/2)+f(b));

d=abs(sqrt(2)/2-T);

disp(d);

>> TX_int(@(x)cos(x),0,pi/4)

   9.5166e-05

ans =

0.7072