         int = findobj(p1.interacts,'elem2',p2);
                        
                        % Compute separation between elements
                        int.kinemat = int.kinemat.setRelPos(p1,p2);
                        
                        % Check if particle is within verlet distance
                        % Assumption: particle do not go directly from
                        % interaction to outside the verlet distance
                        p1.verlet_p(end+1) = p2;
                        