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.251 ]

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

Crivo de Eratóstenes Simples em Fortran 95

CobWeb Plot em Octave

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

Algoritmo para calcular a tabuada

Algorítimo de maior e menor


  

Comentários

Nenhum comentário foi encontrado.


Contribuir com comentário




Patrocínio

Site hospedado pelo provedor RedeHost.
Linux banner

Destaques

Artigos

Dicas

Tópicos

Top 10 do mês

Scripts