function A = bin_numbers_read(filename)
fid = fopen(filename,'r');
if fid < 0
   fprintf('error opening file\n');
   return;
end
dims = fread(fid,2,'uint16');
values = fread(fid,inf,'double');
A = reshape(values,dims(1),dims(2));
fclose(fid);