Octave - Método de Gauss-Legendre

Publicado por Daniel Moreira dos Santos (última atualização em 13/07/2010)

[ Hits: 8.127 ]

Homepage: http://www.danielmoreira.wordpress.com

Download Gauss_Legendre.m




Calcular integral pelo método de Gauss-Legendre usando o Octave.

  



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

Scripts recomendados

Octave - Sistemas lineares por Jacobi

CobWeb Plot em Octave

Matriz de Hilbert e resolução de sistemas lineares

Octave - Calcular raiz pelo método Regula-Falsi

Octave - Calcular raiz pelo método da Secante


  

Comentários

Nenhum comentário foi encontrado.


Contribuir com comentário




Patrocínio

Site hospedado pelo provedor RedeHost.
Linux banner
Linux banner
Linux banner

Destaques

Artigos

Dicas

Tópicos

Top 10 do mês

Scripts