Octave - Método de Gauss-Legendre

Descrição

Esconder código-fonte

function [Integral,CondErro] = Gauss_Legendre(funcao,a,b,n)


%entrada: funcao,a,b,n -> funcao,limite inferior, limite superior, numero de pontos

%saida: Integral,CondErro -> integral e condicao de erro



Integral = 0;

[Avet,Tvet,CondErro] = PesAbsGL(n);

if CondErro ~= 0

    return;

end

e1 = (b-a)/2;

e2 = (a+b)/2;

if (mod(n,2) == 0)

    c1 = 1;

    c2 = 0.5;

else

    c1 =0;

    c2 =1;

end

for i=1:n

    k = fix(i-0.5*(n+1)+sign(i-0.5*(n+c1))*c2);

    t = sign(k)*Tvet(abs(k));

    x = e1*t+e2;

    y = eval(funcao); 

    c = Avet(abs(k));

    Integral = Integral+y*c;

    %escreva i,t,x,y,c

end

Integral = e1*Integral;

end