function [locsn,nein] = RD_gen_cells2D(num_cells,err)
%

x = [1:num_cells];
y = x;
[xx,yy] = meshgrid(x,y);
locs = [xx(:),yy(:)];
locsn = locs + err*(0.5-rand(num_cells^2,2));
nein = RD_neighbors(locsn);
