% function that describes Cheby Type I polynomials in monomial basis function M = M_Cheb_A(N) if N == 0 M = 1; elseif N == 1 M = [1,0;1,0]; else %% Chebyshev's recurrence M = zeros(N+1,N+1); M(1,1)= 1; M(2,2) = 1; for k = 2:N M(k+1,:) = [0,M(k,1:N)]*2 - M(k-1,:); end %% %% end