Processing math: 100%

 

 

 

The Verlet method

If we now perform a Taylor expansion x(t+h)=x(t)+hx(1)(t)+h22x(2)(t)+O(h3).

In our case the second derivative is known via Newton's second law, namely x(2)(t)=a(x,t). If we add to the above equation the corresponding Taylor expansion for x(th), we obtain, using the discretized expressions x(ti±h)=xi±1andxi=x(ti),
we arrive at xi+1=2xixi1+h2x(2)i+O(h4).