Sobre o Método de Runge-Kutta de 4ª Ordem (RK-4)

O método de Runge-Kutta de 4ª ordem (RK-4) é uma das técnicas mais utilizadas para resolver equações diferenciais ordinárias numericamente. Ele fornece excelente precisão sem exigir derivadas de ordem superior, sendo amplamente aplicado em problemas de engenharia e ciências aplicadas.

1️⃣ Versão Clássica do RK-4

Esta é a forma mais conhecida do método, também chamada de forma padrão:

function r = rk4(t, y, h) // RK-4 clássico
    k1 = h * gg(t, y);
    k2 = h * gg(t + h/2, y + k1/2);  
    k3 = h * gg(t + h/2, y + k2/2);
    k4 = h * gg(t + h, y + k3);
    r = y + (k1 + 2*(k2 + k3) + k4)/6;
endfunction
    

2️⃣ Versão RK-4 com Regra 3/8

Outra formulação possível é a chamada 3/8-rule, que usa diferentes coeficientes de ponderação:

function r = rk4f(t, y, h) // RK-4 - 3/8-rule
    k1 = h * gg(t, y);
    k2 = h * gg(t + h/3, y + k1/3);  
    k3 = h * gg(t + 2*h/3, y + (k2 - k1/3));
    k4 = h * gg(t + h, y + (k3 - k2 + k1));
    r = y + (k1 + 3*(k2 + k3) + k4)/8;
endfunction
    

Embora menos comum, essa versão também é de quarta ordem e pode apresentar melhor estabilidade em certos problemas.

3️⃣ Método de Ralston (RK-4)

O Método de Ralston é uma variação otimizada do RK-4, ajustando os coeficientes para minimizar o erro de truncamento local:

function r = rk4g(t, y, h) // RK-4 - Ralston's method
    c3 = (14 - 3*sqrt(5)) / 16;
    a21 = (-2889 + 1428*sqrt(5)) / 1024;
    a22 = (3785 - 1620*sqrt(5)) / 1024;
    a31 = (-3365 + 2094*sqrt(5)) / 6040;
    a32 = (-975 - 3046*sqrt(5)) / 2552;
    a33 = (467040 + 203968*sqrt(5)) / 240845;
    b1 = (263 + 24*sqrt(5)) / 1812;
    b2 = (125 - 1000*sqrt(5)) / 3828;
    b3 = (3426304 + 1661952*sqrt(5)) / 5924787;
    b4 = (30 - 4*sqrt(5)) / 123;

    k1 = h * gg(t, y);
    k2 = h * gg(t + h*2/5, y + k1*2/5);  
    k3 = h * gg(t + c3*h, y + a22*k2 + a21*k1);
    k4 = h * gg(t + h, y + a33*k3 + a32*k2 + a31*k1);
    r = y + (b1*k1 + b2*k2 + b3*k3 + b4*k4);
endfunction
    

Essa versão é mais complexa, mas pode oferecer maior precisão em problemas específicos, especialmente quando o comportamento da função é altamente não linear.

📊 Explicação Visual do Método RK-4

O método RK-4 calcula quatro estimativas intermediárias da derivada — k₁, k₂, k₃ e k₄ — e combina essas informações para obter uma aproximação precisa de y(t+h). O fluxo de cálculo pode ser visualizado da seguinte forma:

y(t+h) = y(t) + (k₁ + 2k₂ + 2k₃ + k₄)/6

📘 Conclusão

O método RK-4 possui diversas formulações (teoricamente podem ser 'infinitas') equivalentes em ordem e precisão, mas com diferentes coeficientes e estratégias de avaliação. Todas seguem o mesmo princípio: combinar múltiplas estimativas da derivada para obter uma aproximação mais precisa da solução da equação diferencial.

Essas três versões ilustram bem a flexibilidade e o poder do método de Runge-Kutta, sendo excelentes exemplos para estudo e implementação prática em Scilab ou outras linguagens numéricas.

---

📌 Notas Finais

1. Minibiografia de Carl Runge

Carl David Tolmé Runge (1856–1927) foi um matemático e físico alemão, conhecido por suas contribuições em análise numérica e espectroscopia. Estudou em Berlim com Karl Weierstrass e Ernst Kummer, tornando-se professor em Hannover e depois em Göttingen. Em parceria com Martin Kutta, desenvolveu o método de Runge-Kutta para integração numérica de equações diferenciais. Também é lembrado pelo Fenômeno de Runge e pelo vetor de Laplace-Runge-Lenz.

2. Minibiografia de Martin Kutta

Martin Wilhelm Kutta (1867–1944) foi um matemático alemão, nascido na Alta Silésia (atual Polônia). Estudou em Breslau e Munique, sendo assistente de Walther von Dyck. Lecionou em Jena, Aachen e Stuttgart, onde permaneceu até sua aposentadoria. Em 1901, desenvolveu junto com Runge o método de Runge-Kutta. Também contribuiu para a aerodinâmica com o Teorema de Kutta-Joukowski e a condição de Kutta.

3. Nota de Transparência

Esta postagem foi elaborada com o auxílio do Copilot, um sistema de inteligência artificial da Microsoft, utilizado aqui como apoio didático para enriquecer o conteúdo da disciplina de Cálculo Numérico.

---

Extra:
% Método de Runge-Kutta de 4ª Ordem (RK-4) % --- RK-4 Clássico --- \[ \begin{aligned} k_1 &= h \cdot g(t, y), \\ k_2 &= h \cdot g\left(t + \frac{h}{2},\, y + \frac{k_1}{2}\right), \\ k_3 &= h \cdot g\left(t + \frac{h}{2},\, y + \frac{k_2}{2}\right), \\ k_4 &= h \cdot g(t + h,\, y + k_3), \\ y_{n+1} &= y_n + \frac{k_1 + 2(k_2 + k_3) + k_4}{6}. \end{aligned} \] % --- RK-4 com Regra 3/8 --- \[ \begin{aligned} k_1 &= h \cdot g(t, y), \\ k_2 &= h \cdot g\left(t + \frac{h}{3},\, y + \frac{k_1}{3}\right), \\ k_3 &= h \cdot g\left(t + \frac{2h}{3},\, y + \left(k_2 - \frac{k_1}{3}\right)\right), \\ k_4 &= h \cdot g\left(t + h,\, y + (k_3 - k_2 + k_1)\right), \\ y_{n+1} &= y_n + \frac{k_1 + 3(k_2 + k_3) + k_4}{8}. \end{aligned} \] % --- RK-4 - Método de Ralston --- \[ \begin{aligned} c_3 &= \frac{14 - 3\sqrt{5}}{16}, \\ a_{21} &= \frac{-2889 + 1428\sqrt{5}}{1024}, \quad a_{22} = \frac{3785 - 1620\sqrt{5}}{1024}, \\ a_{31} &= \frac{-3365 + 2094\sqrt{5}}{6040}, \quad a_{32} = \frac{-975 - 3046\sqrt{5}}{2552}, \quad a_{33} = \frac{467040 + 203968\sqrt{5}}{240845}, \\ b_1 &= \frac{263 + 24\sqrt{5}}{1812}, \quad b_2 = \frac{125 - 1000\sqrt{5}}{3828}, \\ b_3 &= \frac{3426304 + 1661952\sqrt{5}}{5924787}, \quad b_4 = \frac{30 - 4\sqrt{5}}{123}, \\ k_1 &= h \cdot g(t, y), \\ k_2 &= h \cdot g\left(t + \frac{2h}{5},\, y + \frac{2k_1}{5}\right), \\ k_3 &= h \cdot g\left(t + c_3h,\, y + a_{22}k_2 + a_{21}k_1\right), \\ k_4 &= h \cdot g\left(t + h,\, y + a_{33}k_3 + a_{32}k_2 + a_{31}k_1\right), \\ y_{n+1} &= y_n + (b_1k_1 + b_2k_2 + b_3k_3 + b_4k_4). \end{aligned} \]

---

Nota de aula – Cálculo Numérico
Prof. Dr. Francisco J. A. de Aquino