function SumVar = f7(Ng,d,grad)
Ngtilde=1+fix((Ng-1)./grad);%integer part
SumVar=grad.^2*(1/3*Ngtilde.^2-d.*Ngtilde+d.^2-1/3);
end