gauss.m

matlab程序

function I = G_quad(fun,a,b,N)
h = (b-a)/N;I=0;
for k = 1:N
    t = [-sqrt(3/5) 0 sqrt(3/5)];
    A = [5/9 8/9 5/9];
    F = feval(fun,h/2*t+a+(k-1/2)*h);
    I = I + sum(A.*F);
end
I = h/2 * I;