 the verlet distance</span>
                        p1.verlet_w(end+1) = w;

                        <span class="comment">% Check if separation is greater than cutoff distance</span>
                        <span class="comment">% Assumption: cutoff ratio applies to particle radius</span>
                        <span class="keyword">if</span> (int.kinemat.separ &gt;= this.cutoff * p1.radius)
                            <span class="comment">% Remove interaction references from particle</span>
            