this is for holding javascript data
Ghlen Livid edited First_step_in_Velocity_Verlet__.tex
almost 8 years ago
Commit id: a7d9626c8745813f4bacf0ca7d2abf33e5c7d7fc
deletions | additions
diff --git a/First_step_in_Velocity_Verlet__.tex b/First_step_in_Velocity_Verlet__.tex
index 615cc98..b00c635 100644
--- a/First_step_in_Velocity_Verlet__.tex
+++ b/First_step_in_Velocity_Verlet__.tex
...
$$\vec{v}(\Delta t) = \vec{v}(0) + \frac{1}{2}\vec{a}(0)\Delta t + \frac{1}{\Delta t} (\vec{x}(2 \Delta t) - \vec{x}(\Delta t) - \vec{v}(\Delta t) \Delta t)$$
$$2 \vec{v}(\Delta t) = \vec{v}(0) + \frac{1}{2}\vec{a}(0)\Delta t + \frac{1}{\Delta t} (\vec{x}(2 \Delta t) - \vec{x}(\Delta t))$$
Now
$$ \frac{1}{\Delta t} (\vec{x}(2 \Delta t) - \vec{x}(\Delta t)) $$
is a forward finite-difference approximation of $\vec{v}(\Delta t)$, i.e.
$$ \frac{1}{\Delta t} (\vec{x}(2 \Delta t) - \vec{x}(\Delta t)) = \vec{\hat{v}}(\Delta t) + O(\Delta t^2) $$
Again, Euler's iteration has $O(\Delta t^2)$ error.
$$ \vec{v}(0) + \frac{1}{2}\vec{a}(0)\Delta t = \vec{\hat{v}}(\Delta t) + O(\Delta t^2) $$
So we get
$$\vec{v}(\Delta t) = \vec{\hat{v}}(\Delta t) + O(\Delta t^2) $$
where $\vec{\hat{v}}$ is a precise value.