% This little program computes the value of V (Harman's
% notation) given also in formula 7.4 in the text.
bh=lding./(hj*ones(1,nf));
Vtemp=n*sum(sum(bh.^4))-sum(sum(bh.^2).^2);