gauss.m

来自「余向阳博士Matlab教程的源代码。 包括：1.非线性方程的解法；2.线性方程」· M 代码 · 共 11 行

function quad=gauss(f1451,a,b,X,A)
%  f是被积函数。
%  a和b分别为积分的上下限。
%  X是节点的横坐标构成的向量。
%  A是系数构成的向量。
%  输出所求积分值。
N=length(X);
T=zeros(1,N);
T=(a+b)/2+((b-a)/2)*X;
quad=((b-a)/2)*sum(A.*feval('f1451',T));