Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Verlet integration - first iteration

  1. Jun 20, 2012 #1
    I'm trying to make a simple simulation of a projectile motion with air drag. I have already implemented it using Euler and fourth order Runge Kutta methods. I am however unsure about Verlet integration. The equation goes as follows:

    [tex] r(t + Δt) = 2r(t) - r(t-Δt) + a(t)Δt^2[/tex]

    I don't really know what value should I use during first iteration for


    Could someone please help me?
  2. jcsd
  3. Jun 20, 2012 #2

    D H

    User Avatar
    Staff Emeritus
    Science Advisor

    Verlet is one of many "multistep" integration techniques that need past data. As you noted, that past data doesn't exist on the first step. Some other technique, ideally one with comparable error characteristics, needs to be used to bootstrap a multistep integrator. For verlet integration, just use a quadratic estimation for that very first step,
    [tex]r(t+\Delta t) = r(t) + v(t)\Delta t + \frac 1 2 a(t)\Delta t^2[/tex]
    This has a global error of order 2, which is the same as verlet.

    From the second step onward you can use verlet integration.
  4. Jun 20, 2012 #3
    Thank you very much, that really helped!
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook