gausslegendre.m

来自「数值分析常用的一些公式」· M 代码 · 共 12 行

function s=GaussLegendre(a,b)
%用4点Gauss求积公式求数值积分,其中a与b为积分区间
x=[0.8611363116 -0.8611363116 0.3399810436 -0.3399810436];
u=[0.3478548451 0.3478548451 0.6521451549 0.6521451549];
t=0.0;
for i=1:4
    y=x(i)*(b-a)*0.5+(a+b)*0.5;
    t=t+u(i)*ff(y);
end
s=t*(b-a)*0.5;
return