gauss.asv

来自「课程设计：编写连续函数最佳平方逼近的算法程序（以函数形式）。并用此程序进行数值实」· ASV 代码 · 共 14 行

function quad=GLintegration(f)
% Gauss-Lengendre 求积公式
%输入f是被积函数的字符表达式

%输出－－quad是积分值
A=[0.9061798,0.5384693,0,-0.5384693,-0.9061798];% A是1XN的Gauss-Lengendre节点坐标
W=[0.2369239,0.4786286,0.5688889,0.4786286,0.2369239];% W是1XN的Gauss-Lengendre节点对应的权值
a=-1;
b=1;
N=length(A);
T=zeros(1,N);
T=((a+b)/2)+((b-a)/2)*A;
quad=((b-a)/2)*sum(W.*feval(f,T));