function Entropy = f9(Ng,d,grad)
Ngtilde=1+fix((Ng-1)./grad);%integer part
Entropy=log(Ngtilde).*d./d;%because f9 is independent of d we must 
%use *.d./d to get a unit vector of the size of d 
end