第五章数值积分-席.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 1. 建立被积函数文件： function y=f(x) y=sqrt(1 + exp(x)) 2.调用： I = Newt_itg(@f, 0, 2, 5) 例： 四、Gauss－Legendre公式 采用不等距节点（Legendre多项式的根）可以大幅度提高积分精度，在相同的节点数下，其精度是Newton－Cotes公式的两倍。 Legendre多项式： function pn = Legen_pw(n) pbb=[1]; if n==0, pn=pbb; break; end pb=[1 0]; if n==1, pn=pb; break; end for i=2:n; pn= ((2*i-1)*[pb,0] - (i-1)*[0, 0, pbb] )/i; pbb=pb; pb=pn; end 调用： p=Legen_pw(n) 其幂级数形式为： （系数可由Legen_pw(n)确定） n 阶Legendre多项式的 n 个根全部位于-1与1之间，n＝7时的根如下