next up previous
Next: Exercise 1.2: One dimensional Up: Ordinary differential equations: a Previous: 4th order Runge-Kutta

Equations of motion (2nd order ODEs)

We know that the motion of an object is determined by Newton's equations. In the one-dimensional case, we can define the instantaneous position $y(t)$, velocity $v(t)$ and acceleration $a(t)$ of an object using the language of differential calculus:
\begin{displaymath}
v(t)=\frac{dy}{dt},
a(t)=\frac{dv}{dt}.
\end{displaymath} (19)

The motion of the particle is defined by:

\begin{displaymath}
\frac{d^2y}{dt^2}=\frac{F(t)}{m}
\end{displaymath}

This is a second order differential equation that can written as two first order differential equations:
    $\displaystyle \frac{dv}{dt}=\frac{F(t)}{m};$ (20)
    $\displaystyle \frac{dy}{dt}=v(t).$ (21)

To solve it we can apply any of the methods described in the previous sections. If we pick Euler's, we obtain:
    $\displaystyle v_{n+1}=v_n+\frac{F(t)}{m}\Delta t = v_n+a_n\Delta t,$ (22)
    $\displaystyle y_{n+1}=y_n+v_n\Delta t,$ (23)

where $a_n=F(t)/m$.



Subsections
next up previous
Next: Exercise 1.2: One dimensional Up: Ordinary differential equations: a Previous: 4th order Runge-Kutta
Adrian E. Feiguin 2004-06-01