Octave - Método de Runge-Kutta

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

[ Hits: 16.267 ]

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

Matriz de Hilbert e resolução de sistemas lineares

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

Octave - Sistemas lineares por Jacobi

Octave - Calcular raiz pelo método Regula-Falsi

Octave - Método de Gauss-Legendre


  

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