mpty(BB),
0242     qb=BB(:,1);
0243     <span class="keyword">if</span> size(BB,2)&gt;1, qb=qb+BB(:,2)*fac; <span class="keyword">end</span>
0244     [jj,ii]=ismember(qb,qel);
0245 <span class="comment">%     [aa,bb]=meshgrid(qb,qel);</span>
0246 <span class="comment">%     [ii,jj]=find((aa-bb)==0);</span>
0247     N.b(jj)=ii;
0248 <span class="keyword">end</span>
0249 [jj,ii]=ismember(qm,qel);
0250 <span class="comment">% [aa,bb]=meshgrid(qm,qel);</span>
0251 <span class="comment">% [ii,jj]=find((aa-bb)==0