function [total]= histdiffnorm(a,b);
A = imread(a);
B = imread(b);
M = imhist(A);
N = imhist(B);
total = 0;
for i = 1:256
    total = total+sqrt((M(i,:)-N(i,:))*(M(i,:)-N(i,:)));
end
return;
