31 
0232 fac=111111;
0233 qel=N.elec(:,1)+N.elec(:,2)*fac;
0234 qa=AA(:,1);qm=MM(:,1);
0235 <span class="keyword">if</span> size(AA,2)&gt;1, qa=qa+AA(:,2)*fac; <span class="keyword">end</span>
0236 <span class="keyword">if</span> size(MM,2)&gt;1, qm=qm+MM(:,2)*fac; <span class="keyword">end</span>
0237 [jj,ii]=ismember(qa,qel);
0238 <span class="comment">% [aa,bb]=meshgrid(qa,qel);</span>
0239 <span class="comment">% [ii,jj]=find((aa-bb)==0);</span>
0240 N.a(jj)=ii;
0241 <span class="keyword">if</span> ~ise