How to count number of points in a volumetric grid (4D histogram)?

Hi, I have 3D point data in XYZ and i wish to creat a 4D histogram, i.e. a 3D matrix containing the point count in a volumetric lattice. now, i know how to do it using a ton of for loops but the question is what is th good implementation for this ? btw i found this : https://www.mathworks.com/matlabcentral/fileexchange/23897-n-dimensional-histogram

However it uses histc (which i saw is not recommended) and it gives me a strange result (for some unkown reason it adds a bin in one of the dimensions) Thanks, Itzik

[Itzik Ben Shabat]

I changed line 108 in the link to

   [~,~, loc(:,d)] = histcounts (Xd, ed);

I think it worked

[Mikk Laanes]

Could you share how you made the referenced code to work in 3D? It appears to only accept M-by-N matrices?