Octave - Método de Euler para solução de EDO
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.
Descrição
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.
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