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.257 ]

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 - Calcular raiz pelo método de Newton

CobWeb Plot em Octave

Crivo de Eratóstenes Simples em Fortran 95

Octave - Sistemas lineares por Jacobi

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


  

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