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

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

[ Hits: 4.487 ]

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

Download Euler_Melhorado.m




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




Esconder código-fonte

function [vetx,vety] = Euler_Melhorado(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

    xa = x;

    ya = y;

    x = x+h;

    y = y+h*Fxy;

    ter = eval(funcao);

    y = ya + (h/2)*(Fxy+ter);

    x= a+i*h;

    Fxy = eval(funcao);    

    disp([i  x    y    Fxy]);

    vetx(i+1) = x;

    vety(i+1) = y;

end 

end



    

Scripts recomendados

Octave - Método de Gauss-Legendre

Octave - Método de Runge-Kutta

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

Algorítimo de maior e menor

Crivo de Eratóstenes Simples em Fortran 95

Comentários

Nenhum comentário foi encontrado.


Contribuir com comentário