chebyshev polynomials code snippets in matlab