Octave - Método de Euler para solução de EDO

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

[ Hits: 21.784 ]

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

Download Euler.m




Em matemática e ciência computacional, o método de Euler, cujo nome relaciona-se com Leonhard Euler, é um procedimento numérico de primeira ordem para solucionar equações diferenciais ordinárias com um valor inicial dado. É o tipo mais básico de método explícito para integração numérica para equações diferenciais ordinárias.

  



Esconder código-fonte

function [vetx,vety] = Euler(funcao,a,b,m,y0);


%parametros de entrada: funcao,a,b,m,y0 -> funcao,limites inferior, superior, num. de

%intervalos e valor inicial

%parametros de saida: vetx, vety -> abcissas e solucao do PVI



h = (b-a)/m;

x=a;

y = y0;

Fxy = eval(funcao);

vetx(1) = x;

vety(1) = y;

disp('     i      x     y     Fxy  ');

disp([0  x    y    Fxy]);

for i=1:m

    x = a+i*h;

    y = y+h*Fxy;

    Fxy = eval(funcao);  

    disp([i  x    y    Fxy]);

    vetx(i+1) = x;

    vety(i+1) = y;

end

end



    

Scripts recomendados

Octave - Sistemas lineares por Gauss-Seidel

Octave - Método de Runge-Kutta

Octave - Método de Euler-Modificado para solução de EDO

Octave - Sistemas lineares por Jacobi

Crivo de Eratóstenes Simples em Fortran 95


  

Comentários

Nenhum comentário foi encontrado.


Contribuir com comentário




Patrocínio

Site hospedado pelo provedor RedeHost.
Linux banner

Destaques

Artigos

Dicas

Tópicos

Top 10 do mês

Scripts