 <span class="keyword">if</span> (vparticle &gt; vmax)
                        vmax = vparticle;
                    <span class="keyword">end</span>
                    <span class="keyword">if</span> (p1.radius &gt; rmax)
                        rmax = p1.radius;
                    <span class="keyword">end</span>
                <span class="keyword">end</span>

                <span class="comment">% Reset verlet lists</span>
                p1.verlet_p = Particle.empty;
                p1.verlet_w = W