      <span class="comment">% Assumption: particle do not go directly from</span>
                        <span class="comment">% interaction to outside the verlet distance</span>
                        p1.verlet_p(end+1) = p2;

                        <span class="comment">% Check if separation is greater than cutoff distance</span>
                        <span class="comment">% Assumption: cutoff ratio applies to maximum radius</span>
                        <span class="keyword">if</span> (int.kinemat.