function z = cheb(x, order)
if (nargin<2)
    order=14;
end

x = reshape(x,length(x),1);
z(:,1) = ones(size(x));
z(:,2) = x;
if (order<3)
    return
end
twox = 2.0*x;
for (n=3:order)
    z(:,n)=z(:,n-1).*twox-z(:,n-2);
end

Published with MATLAB® 7.6