function quad=gauss(f1451,a,b,X,A)
%  f是被積函數(shù)。
%  a和b分別為積分的上下限。
%  X是節(jié)點(diǎn)的橫坐標(biāo)構(gòu)成的向量。
%  A是系數(shù)構(gòu)成的向量。
%  輸出所求積分值。
N=length(X);
T=zeros(1,N);
T=(a+b)/2+((b-a)/2)*X;
quad=((b-a)/2)*sum(A.*feval('f1451',T));