Octave - Método de Runge-Kutta

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

[ Hits: 10.921 ]

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

Download RungeKutta.m




Em análise numérica, os métodos de Runge–Kutta formam uma família importante de metódos iterativos implícitos e explícitos para a resolução numérica (aproximação) de soluções de equações diferenciais ordinárias. Aqui, vamos resolver o PVI por Runge-Kutta de ordem 4.

  



Esconder código-fonte

function [Vetx,Vety] = RungeKutta(funcao,a,b,m,y0)


%parametros de entrada: funcao,a,b,m,y0 - > lim. inf., lim. sup., num.

%subintervalos e valor inicial

%parametros de saida: Vetx,Vety -> abcissas e solucao do PVI



h = (b-a)/m;

xt = a;

yt = y0;

Vetx(1) = xt;

Vety(1)=yt;

disp ('   i   xt    yt');

disp([ 0   xt   yt]);

for i=1:m

    x = xt;

    y = yt;

    k1 = eval(funcao);

    x = xt + h/2;

    y = yt + h/2 * k1;

    k2 = eval(funcao);

    y = yt + h/2 * k2;

    k3 = eval(funcao);

    x = xt + h;

    y = yt + h * k3;

    k4 = eval(funcao);

    xt = a+i*h;

    yt = yt+h/6*(k1+2*(k2+k3)+k4);

    disp([ i   xt   yt]);

end

end


Scripts recomendados

Octave - Calcular raiz pelo método Regula-Falsi

Função Fatorial Recursiva em Fortran 95

CobWeb Plot em Octave

Octave - Calcular raiz pelo método de Newton

Octave - Método de Euler-Modificado 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