wolffd@0: function x = logb(y, base)
wolffd@0: % logb - logarithm to base b
wolffd@0: % function x = logb(y, base)
wolffd@0: 
wolffd@0: x = log(y) ./ log(base);