function [KN] = normalise(K);

KN = repmat(1./sqrt(diag(K)),1,size(K,1)).*K.*repmat((1./sqrt(diag(K)))',size(K,1),1);
