function y = horner(c,t)
% Evaluation of polynomial in monomial basis representation.
% Note: c(1) is the leading coefficient, c(end) the constant term
% (This follows the convention of MATLAB's ployval() function
y = c(1)*ones(size(t)); for j=2:length(c), y = t.*y+c(j); end