function are not permitted in this context, a am using matlab R2012

function[d] = hcompare_EMD(h1,h2)

d = sum(abs(cumsum(h1) - cumsum(h2))); end

img = imread('C:\Users\Hp\Desktop\col\001.jpg');

function[h] = histImage(img)

% This function calculates normalized histogram of image.

% Normalized histogram is histogram h, that has at

% each i place in it, value:

% (number of picture pixels with gray level i-1) /

% (total num of pixels in picture).

sum = 0;

[y,x] = size(img); % getting sizes of image

h = zeros(1,256); % creating output histogram array

for i = 1:1: y % runing on rows

    for j = 1:1: x % running on colomns

        % gray level is addtess to cell in output histogram array

   % we add there 1 (bacause of current pixel (y,x) has this gray level

h(img(i,j)) = h(img(i,j)) + 1;

        % pay attention to fact, that we use here pixel value as index!
    end
end

h = h./(y*x);

end % to calculate the distance between the histograms of the two images %% %(histArray, histPattern):

function[dmap] = patDistMAp(histArray, histPattern)

[y,x,z] = size(histArray);

dmap = zeros(y,x); % output array

for i = 1:1: y % runing on rows

    for j = 1:1: x   % running on colomns

        hist = histArray(i,j,:);

% for k = 1:1:256 % making array 1x256 from 1x1x256 % h1(k) = hist(1,1,k); % there is a permute function, % end % but we will use it next time)

        h1 = permute(squeeze(hist),[2,1]);

        % Using temp variable, as MATLAB7 wants it:

        temp = hcompare_EMD(histPattern,h1);

        dmap(i,j) = temp;

    end

end

end