function y = logchev(n, x) m = size(x, 2); if x(1) < 1 y(1) = logchev(n, x(1).*exp(1)) - 1; elseif x(1) > exp(1) y(1) = logchev(n, x(1)./exp(1)) + 1; else y(1) = chebyshev(n, @log, 1, exp(1), x(1)); end if m > 1 y(2:m) = logchev(n, x(2:m)); end