gausslegendre.m

来自「《MATLAB实用指南》系列丛书DE源代码」· M 代码 · 共 12 行

function s=GaussLegendre(a,b)
  % 用4点Gauss-Legendre求积公式求数值积分，其中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