gauss10.m

来自「学matlab还有比这个更好的吗?我千兴万苦才找到的,这是一个学习matlab的」· M 代码 · 共 12 行

function g = gauss10(fun,a,b)
%GAUSS10(fun,a,b) 	利用10参数Gauss求积法近似计算a<x<b 区间上的积分
%  				fun		被积函数。字符串、内联函数名、M函数文件名都可。
%======================================================
x = [0.1488743390;0.4333953941;0.6974095683;...
     0.8650633667;0.9739065285];
w = [0.2955242247;0.2692667193;0.2190863625;...
     0.1494513492;0.0666713443];
t = .5*(b+a)+.5*(b-a)*[-flipud(x);x];
W = [flipud(w);w];
g = sum(W.*feval(fun,t))*(b-a)/2;