function [Hx_y,Hy_x]=cp(a)

y_marg=sum(a);
x_marg=sum(a');
[r,c] = size(a);

Hx_y=0;
for i=1:r;    %  rows 
   for j=1:c;   %   cols
      if( a(i,j)/y_marg(1,j)==0 )
         %do nothing
      else
         Hx_y = Hx_y + (-a(i,j)*log2(a(i,j)/y_marg(1,i)));
      end
  end
end

Hy_x=0;
for i=1:r;    %  rows 
   for j=1:c;   %   cols
      if( a(i,j)/y_marg(1,j)==0 )
         %do nothing
      else
         Hy_x = Hy_x + (-a(i,j)*log2(a(i,j)/x_marg(1,i)));
      end
  end
end