</em></span>. To
        apply numerical integration one first has to design the right hand side of
        the equation <span class="emphasis"><em>w' = f(w) = (p,-q+&#947; p)</em></span>:
      </p>
<p>
</p>
<pre class="programlisting"><span class="comment">/* The type of container used to hold the state vector */</span>
<span class="keyword">typedef</span> <span class="identifier">std</span><span class="special">::</span><span class="identifier">vector</span><span class="special">&lt;</span> <span class="