..</span>
0033         N.elec(N.m(fb),1:2)+N.elec(N.n(fb),1:2))/4;
0034 <span class="keyword">end</span>
0035 midx=unique(sort(midpoint(:,1)));
0036 midy=unique(sort(midpoint(:,2)));
0037 dmx=min(diff(midx));
0038 midx=min(midx):dmx:max(midx);
0039 dmy=min(diff(midy));
0040 midy=min(midy):dmy:max(midy);
0041 [aa,bb]=meshgrid(midpoint(:,1),midx);
0042 [ix,jx]=find((aa-bb)==0);
0043 [aa,bb]=meshgrid(midpoint(:,2),midy);
0044 [iy,jy]=find((aa-bb)==0);
0045 datums=ones(length(midx),length(midy),length(konfs))*N