guass1.m

来自「Gauss-Legendre 采用五点 Gauss-Legendre 求积公式计」· M 代码 · 共 7 行

function s=guass1(a,b,n)
h=(b-a)/n;
s=0.0;
for m=0:(1*n/2-1)
s=s+h*(guassf(a+h*((1-1/sqrt(3))+2*m))+guassf(a+h*((1+1/sqrt(3))+2*m)));
end