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;