% specdist - spectral distance using I-S divergence
D=sum(X./V - 1 + log(V./X),1);

