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

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

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

Algorítimo de maior e menor

Octave - Decomposição LU

Octave - Método de Gauss-Legendre

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


  

Comentários

Nenhum comentário foi encontrado.


Contribuir com comentário




Patrocínio

Site hospedado pelo provedor HostGator.
Linux banner
Linux banner
Linux banner

Destaques

Artigos

Dicas

Viva o Android

Tópicos

Top 10 do mês

Scripts