guass1.m

Gauss-Legendre 采用五点 Gauss-Legendre 求积公式计算定积分

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