span class="comment">% D: full overlap of domain and element (right)</span>
0042   Svp = 1-(xp-xv)/h;  
0043   dSvp=-1/h;
0044 <span class="keyword">elseif</span> ( h-lp)&lt;(xp-xv) &amp;&amp; (xp-xv)&lt;=( h+lp)                                  <span class="comment">% E: partial overlap of domain and element</span>
0045   Svp = (h+lp-(xp-xv))^2/(4*h*lp); 
0046   dSvp=-(h+lp-(xp-xv))/(2*h*lp);
0047 <span class="keyword">else</span>                                                                        <span