function x = logb(y, base)
% logb - logarithm to base b
% function x = logb(y, base)

x = log(y) ./ log(base);