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

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

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

[ Hits: 10.591 ]

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

Download Euler_Modificado.m

0 0

Denuncie Favoritos Indicar

Resolver o PVI pelo método de Euler Modificado usando o Octave.

Esconder código-fonte

function [vetx,vety] = Euler_Modificado(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  ');
for i=1:m
    yi = y;
    x = x+h/2;    
    y = y+(h/2)*Fxy;
    Fxy = eval(funcao);
    y = yi + h*Fxy; 
   
    disp([i  x    y    Fxy]);
    vetx(i+1) = x;
    vety(i+1) = y;
end
end

Scripts recomendados

Cálculo da raiz de uma função pelo Método da Bisseção - Octave

Octave - Sistemas lineares por Gauss-Seidel

Script em Bath "DOS" para Upload e Download de arquivos usando o SSH

Octave - Calcular raiz pelo método de Newton

Octave - Decomposição LU

Comentários

Nenhum comentário foi encontrado.