samer@4: function k=log2pie
samer@4: 	persistent K;
samer@4: 	if isempty(K) K=1+log(2*pi); end
samer@4: 	k=K;
samer@4: end