1x
foria1
b
xjai1ba
左点
xiaiba
右点
yjsubsyxxjyisubsyxxi
左点值右点值
ssyiyjba2
e
d
si
tegratei
ty12i
tegrate为matlab中自带的积分函数
i
tegratedoublei
tegrateabssi
tegratei
tegrate相对误差
【调试结果】s121y120
i
tegratelog2i
tegrate
0693147180559945a
s
fabs136233888727957504y936497215224803851
○2抛物线法:
formatlo
g
a1b2s0
120抛物线条数120
小区间个数2
symsxy
y1x
foria1
bx0a2iba2
第一点
x1a2i1ba2
第二点x2a2i0ba2
第sa
点y0subsyxx0第一点值
y1subsyxx1第二点值y2subsyxx2第三点值
ssy04y1y2ba6
e
d
si
tegratei
ty12i
tegrate为matlab中自带的积分函数
i
tegratedoublei
tegrateabssi
tegratei
tegrate【调试结果】
相对误差
s
121y120
i
tegrate
log2
i
tegrate
0693147180559945
a
s
abs136233888727957504y936497215224803851
○3使用函数trapz
x111202y1xtrapzxy【调试结果】a
s
0693151520800048
○4使用函数quad
quad1x121120【调试结果】a
s
0693147292908988
f第3题
○1使用函数trapz
x11120i
fysi
xxtrapzxy【调试结果】Errorusi
gcolo
Maximumvariablesizeallowedbytheprogramisexceeded
○2使用函数quad
quadsi
xx0i
f【调试结果】a
s
NaNNaN不定值
第6题
○1矩形法:利用求和函数
利用sum函数改写矩形法formatlo
g
100a0b1symsxfxfx11x2i1
xjai1ba
xiaiba
fxjsubsfxxxj
fxisubsfxxxifxijfxifxj2mfxjba
pfxiba
kfxijba
i
um1summi
um2sumpi
um3sumk【调试结果】i
um10787893996730782i
um20782893996730782i
um30785393996730783
f○2抛物线法:使用求和函数
利用sum函数改写抛物线法formatlo
g
100a0b1symsxfxfx11x2i0
1xja2iba2
xia2i1ba2
xka2i2ba2
fxjsubsfxxxjfxisubsfxxxifxksubsfxxxkmfxj4fxifxkba6
i
umsumm
【调试结果】i
um
0785398163397448
【情况记录】
1、梯形法和抛物线法程序设计较为顺利。但要注意使用for循环函数和求和函数时的
不同matlab命令,避免混淆出错。使用函数trapz,quad时要注意被积函数是数值形
式,应使用数组计算,应用点除即,否则将出错,不能调试出结果。
2、使用函数trapz,quad和附录程序求解,均不能调试出获得出正确答案。最后
尝试用matlab命令中的符号求积分才得出正确结果。
3、参照附录B中的求和函数程序设计顺利改变了附录A和C。发现使r