--- a
+++ b/Supporting Functions/CRI_ShannonEntropy.m
@@ -0,0 +1,12 @@
+function [ H, prob ] = CRI_ShannonEntropy(x, bin_min, bin_max, binWidth)
+%UNTITLED Summary of this function goes here
+%   Detailed explanation goes here
+[counts,binCenters] = hist(x,[bin_min:binWidth:bin_max]);
+binWidth = diff(binCenters);
+binWidth = [binWidth(end),binWidth]; % Replicate last bin width for first, which is indeterminate.
+nz = counts>0; % Index to non-zero bins
+prob = counts(nz)/sum(counts(nz));
+H = -sum(prob.*log2(prob./binWidth(nz)));
+
+end
+