Os métodos de passo múltiplo são ferramentas poderosas para a resolução numérica de Problemas de Valor Inicial (PVI). Diferentemente dos métodos de passo simples (como Euler e Runge-Kutta), que utilizam apenas a informação no ponto atual para avançar, os métodos de passo múltiplo aproveitam informações de pontos anteriores para obter uma maior eficiência computacional e precisão. Hoje, vamos explorar o Método de Adams-Bashforth de quarta ordem, um dos esquemas mais utilizados na prática.
Fundamentação Teórica
Consideremos o problema de valor inicial padrão:
$$
\begin{cases}
y^\prime (x) = f(x,y(x)), \,\,\, x \in [a,b] \\
y(x_0) = y_0,
\end{cases}
$$
onde $a = x_0 < x_1 < \cdots < x_n = b$ são pontos igualmente espaçados com passo $h = x_{i+1} - x_i$. Integrando a equação diferencial no intervalo $[x_i, x_{i+1}]$, obtemos:
$$
y(x_{i+1}) = y(x_i) + \int_{x_i}^{x_{i+1}} f(x, y(x)) \, dx. \quad (1)
$$
A ideia central dos métodos de Adams-Bashforth é aproximar a integral em (1) substituindo a função $f(x, y(x))$ por um polinômio interpolador que passa por pontos anteriores. Para o método de quarta ordem, utilizamos um polinômio de grau 3 que interpola $f$ nos nós $x_i, x_{i-1}, x_{i-2}$ e $x_{i-3}$. Integrando esse polinômio no intervalo $[x_i, x_{i+1}]$, chegamos à seguinte fórmula:
$$
y_{i+1} = y_i + \frac{h}{24} \left[ 55 f(x_i, y_i) – 59 f(x_{i-1}, y_{i-1}) + 37 f(x_{i-2}, y_{i-2}) – 9 f(x_{i-3}, y_{i-3}) \right], \quad i \geq 3.
$$
Esta é a fórmula do Método de Adams-Bashforth de quarta ordem. Observe que ela é explícita, pois $y_{i+1}$ aparece apenas no lado esquerdo da equação. A ordem do método refere-se à ordem do erro global, que é proporcional a $h^4$, tornando-o uma excelente escolha para problemas que exigem alta precisão.
Erro Local
O erro de truncamento local (ETL) do método de Adams-Bashforth de quarta ordem é dado por:
$$
\text{ETL} = \frac{251}{720} h^5 y^{(5)}(\eta_i),
$$
onde $\eta_i \in (x_{i-3}, x_{i+1})$. Este erro local é de ordem $h^5$, o que garante que o erro global (acumulado) seja de ordem $h^4$.
Inicialização do Método
Um ponto crucial é que, para $i = 3$, a fórmula exige os valores de $y_1, y_2$ e $y_3$, além da condição inicial $y_0$. Como o método não é autoiniciante, precisamos obter esses valores iniciais utilizando um método de passo simples de alta ordem, como o Método de Runge-Kutta de quarta ordem (RK-4). Essa combinação é comum e eficiente em aplicações práticas.
Implementação Computacional
Para facilitar a aplicação do método, desenvolvemos um solver completo em JavaScript que implementa o Adams-Bashforth de quarta ordem com inicialização via RK-4. O código está bem comentado e é de fácil adaptação para diferentes problemas.
Você pode baixar ou visualizar o arquivo completo com o solver clicando no link abaixo:
📄 Acesse o Solver Adams-Bashforth (adamsbashforth.html)
Exemplo Numérico
Vamos testar o método com o seguinte problema clássico:
$$
\begin{cases}
y^\prime(x) = y(x), \,\,\, x \in [0,1] \\
y(0) = 1,
\end{cases}
$$
cuja solução exata é $y(x) = e^x$. Utilizando o solver com $N = 10$ subintervalos ($h = 0.1$), obtemos uma aproximação extremamente precisa, com erros da ordem de $10^{-6}$. Os resultados demonstram a superioridade do método em relação a esquemas de ordem inferior.
Vantagens e Aplicações
- Eficiência: Como é um método explícito, cada passo requer apenas uma avaliação de $f(x,y)$, sendo computacionalmente mais barato que métodos implícitos.
- Alta Precisão: A quarta ordem garante resultados muito bons com passos moderados.
- Ampla Utilização: É empregado em simulações de dinâmica de fluidos, astrofísica, engenharia e qualquer área que exija a solução de EDOs com alta eficiência.
No entanto, é importante lembrar que, por ser explícito, pode sofrer com instabilidades para problemas rígidos. Nesses casos, métodos implícitos como Adams-Moulton são mais adequados.
Conclusão
O Método de Adams-Bashforth de quarta ordem é uma excelente opção para a resolução numérica de PVI quando se busca um equilíbrio entre precisão e custo computacional. Sua formulação simples e seu desempenho o tornam uma ferramenta indispensável no kit de qualquer cientista computacional.
Você pode praticar e visualizar seus exemplos utilizando o arquivo completo com o solver clicando no link abaixo:
📄 Acesse o Solver Adams-Bashforth (adamsbashforth.html)