function DiffVar = f10(Ng,d,grad)
Ngtilde=fix((Ng-1)./grad)+1-d;
DiffVar = grad.^2.*d.*(Ngtilde-d);
end