Método de Verlet

Para o método de Euler implícito havíamos utilizado a derivada a esquerda:

$$f^{\prime }\left(t\right)\approx \frac{f\left(t\right)-f(t-\mathrm{\Delta }t)}{\mathrm{\Delta }t}$$

Então se $y\left(t\right)=f^{\prime }\left(t\right)$ a segunda derivada é $y^{\prime }\left(t\right)=f^{″}\left(t\right)$, pela definição, da derivada a direita:

$$y^{\prime }\left(t\right)=\underset{\mathrm{\Delta }t\to 0}{\lim }\frac{y(t+\mathrm{\Delta }t)-y\left(t\right)}{\mathrm{\Delta }t}=\underset{\mathrm{\Delta }t\to 0}{\lim }\frac{f^{\prime }(t+\mathrm{\Delta }t)-f^{\prime }\left(t\right)}{\mathrm{\Delta }t}$$

Logo utilizando as aproximações:

$$\begin{array}{rl}{f}^{″}\left(t\right)& \approx \frac{1}{\mathrm{\Delta }t}(f^{\prime }(t+\mathrm{\Delta }t)-f^{\prime }\left(t\right))\\ & \approx \frac{1}{\mathrm{\Delta }t}(\frac{f(t+\mathrm{\Delta }t)-f(t+\mathrm{\Delta }t-\mathrm{\Delta }t)}{\mathrm{\Delta }t}-\frac{f\left(t\right)-f(t-\mathrm{\Delta }t)}{\mathrm{\Delta }t})\\ & \approx \frac{f(t+\mathrm{\Delta }t)+f(t-\mathrm{\Delta }t)-2f\left(t\right)}{{\left(\mathrm{\Delta }t\right)}^2}\end{array}$$

Isolando então $f(t+\mathrm{\Delta }t)$:

$$\begin{array}{rl}f(t+\mathrm{\Delta }t)& ={\left(\mathrm{\Delta }t\right)}^2{f}^{″}\left(t\right)-f(t-\mathrm{\Delta }t)+2f\left(t\right)\end{array}$$

Temos o método de Verlet. Podemos notar que precisamos conhecer $f\left(t\right)$ em dois tempos anteriores. Podemos utilizar outro algoritmo para o primeiro passo. Se $f\left(t\right)=x\left(t\right)$ é posição, então $f^{″}\left(t\right)=\ddot{x}\left(t\right)=a(x,t)$ logo podemos reescrever:

$$\begin{array}{rl}x(t+\mathrm{\Delta }t)& =a(x,t){\left(\mathrm{\Delta }t\right)}^2-x(t-\mathrm{\Delta }t)+2x\left(t\right)\end{array}$$

Para calcular a energia, podemos obter a velocidade então utilizando a derivada centrada: $$v\left(t\right)=\dot{x}\left(t\right)=\frac{x(t+\mathrm{\Delta }t)-x(t-\mathrm{\Delta }t)}{2\mathrm{\Delta }t}$$

Principais materiais utilizados

1. The Second Derivative (Matthew Boelkins, David Austin & Steven Schlicker; LibreTexts)