(一)专题实验(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