It's not clear where the bin limits actually are, but:
% the data
x = [1.33 1.89 9.27 9.46 9.20 7.43 6.08 5.57 6.70];
y = [8.89 0.77 1.49 9.36 8.69 1.61 1.34 4.60 2.77];
% edge breakpoints
tiling = [10 10]; % [x y]
xedges = linspace(min(x),max(x),tiling(1));
yedges = linspace(min(y),max(y),tiling(2));
% the member counts
counts = histcounts2(x,y,xedges,yedges)
% visualize the same thing
histogram2(x,y,xedges,yedges);
