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: 9.215 ]

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

Download Euler_Modificado.m




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

Função Fatorial Recursiva em Fortran 95

Octave - Calcular raiz pelo método Regula-Falsi

Crivo de Eratóstenes Simples em Fortran 95

Teste de hardware e outros

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
Linux banner
Linux banner

Destaques

Artigos

Dicas

Tópicos

Top 10 do mês

Scripts