les(i);

                <span class="comment">% Update maximum velocity</span>
                <span class="keyword">if</span> (this.auto_freq)
                    vparticle = vecnorm(p1.veloc_trl);
                    <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 c