function SumAvg = f6(Ng,d,grad)
Ngtilde=1+fix((Ng-1)./grad);%integer part
SumAvg=(Ngtilde-1).*grad.*d./d;%although f6 is independent of d
%we must used *.d./d to get unit vecor of the size of d
end