er shapes, the base kinematic object will be used.</span>
            dir   = p2.coord - p1.coord;
            dist  = norm(dir);
            separ = dist - p1.radius - p2.radius;

            <span class="comment">% Check if particle is within verlet distance</span>
            <span class="keyword">if</span> (addVL &amp;&amp; separ &lt; this.verlet_dist)
                p1.verlet_p(end+1) = p2;
            <span class="keyword">end</span>

            <span class="comment">% Check if interaction exists</s