+ S * fmax;
            <span class="keyword">end</span>

            <span class="comment">% Residual overlap</span>
            this.residue = v0 * sqrt(m/this.stiff) * (1-e2);
        <span class="keyword">end</span>

        <span class="comment">%------------------------------------------------------------------</span>
        <span class="keyword">function</span> this = evalForce(this,int)
            <span class="comment">% Needed properties</span>
            dir  = int.kinemat.dir_n;
            ov