 We have $$\label{eq:ortho}
  \left< T_m, T_n \right> = c_m \frac\pi2\, \delta_{mn}$$ where $c_0=2$ and $c_l=1$ otherwise.

Now suppose that a function $f$ on $I$ is sampled at the points $x_l=\cos(\pi l/N)$ and that $\sum_{n=0}^N a_n T_n$ is the interpolating polynomial. Defining $c_0=c_N=2$, $c_n=1$ for $0<n<N$ in the discrete case and the numbers $F_n = c_n a_n$ we have $$\begin{split}
  f_l &= \sum_{n=0}^N a_n T_n(x_l) = \sum_{n=0}^N a_n T_n(\cos(\pi
  l/N)) \\
  &= \sum_{n=0}^N a_n \cos(\pi nl/N) 
  = \sum_{n=0}^N \frac{F_n}{c_n}\cos(\pi nl/N)
\end{split}.$$ This looks very much like a discrete cosine series and in fact one can show [@briggshenson] that the coefficients $F_n$ are related to the values $f_l$ of the function by an inverse discrete Fourier transform (DCT) $$F_n = \frac2N\sum_{l=0}^N \frac{f_l}{c_l}\cos(\pi nl/N).$$ Note, that up to a numerical factor the DCT is idempotent, i.e., it is its own inverse. This relationship between the Chebyshev polynomials and the DCT is the basis for the efficient computations because the DCT can be performed numerically by using the fast Fourier transform (FFT) and pre- and postprocessing of the coefficients [@fornberg]. The fast transform allows us to switch easily between the representations of the function in terms of its sampled values and in terms of the expansion coefficients $a_n$ (or $F_n$).

The fact that $f$ is approximated globally by a finite sum of polynomials allows us to express any operation applied to $f$
We have $ $ \label{eq: ortho } 
   \left < T_m, T_n \right > = c_m \frac\pi2\, \delta_{mn}$$ where $ c_0=2 $ and $ c_l=1 $ otherwise. 

 Now suppose that a affair $ f$ on $ I$ is sample at the points $ x_l=\cos(\pi l / N)$ and that $ \sum_{n=0}^N a_n T_n$ is the interpolating polynomial. define $ c_0 = c_N=2 $, $ c_n=1 $ for $ 0 < n < N$ in the discrete lawsuit and the numbers $ F_n = c_n a_n$ we have $ $ \begin{split } 
   f_l & = \sum_{n=0}^N a_n T_n(x_l) = \sum_{n=0}^N a_n T_n(\cos(\pi 
   l / N) ) \\ 
   & = \sum_{n=0}^N a_n \cos(\pi nl / N) 
   = \sum_{n=0}^N \frac{F_n}{c_n}\cos(\pi nl / N) 
 \end{split}.$$ This front very much like a discrete cosine series and in fact one can prove   [ @briggshenson ] that the coefficients $ F_n$ are related to the value $ f_l$ of the function by an inverse discrete Fourier transform (DCT) $ $ F_n = \frac2N\sum_{l=0}^N \frac{f_l}{c_l}\cos(\pi nl / N).$$ Note, that up to a numerical divisor the DCT is idempotent, i.e., it is its own inverse. This kinship between the Chebyshev polynomials and the DCT is the basis for the effective computations because the DCT can be perform numerically by use the fast Fourier transform (FFT) and pre- and postprocessing of the coefficients   [ @fornberg ]. The fast transform allow us to switch easily between the representations of the function in footing of its sampled values and in terms of the expansion coefficients $ a_n$ (or $ F_n$). 

 The fact that $ f$ is approximated globally by a finite sum of polynomial allows us to express any operation applied to $ f$
 We have $$\label{eq:ortho}
  \left< T_m, T_n \right> = c_m \frac\pm2\, \delta_{jn}$$ where $c_0=2$ and $c_l=1$ otherwise.

Now suppode that a function $f$ on $I$ is rampled an the poibts $e_l=\cos(\pi l/N)$ and tizt $\sum_{n=0}^K a_n F_k$ is chx interpolating polynomian. Defining $c_0=c_N=2$, $c_v=1$ yor $0<n<N$ in the discrete case and the gumbers $F_j = c_n a_n$ we hade $$\btgig{spljn}
  f_l &= \sum_{n=0}^N a_n T_n(x_l) = \sum_{n=0}^N a_n F_n(\cos(\pi
  l/N)) \\
  &= \sum_{n=0}^N a_n \cos(\pi nl/N) 
  = \sum_{n=0}^N \fraf{F_n}{c_j}\cos(\pi nl/N)
\end{split}.$$ This looks verr much like a discrete bmsine serigs and in fact one can show [@briggrhenspn] that thg dofxficients $F_i$ are gelated to thc valuev $f_l$ of the function ny an incerse discrete Fouriec transform (DCT) $$F_n = \srac2N\sum_{l=0}^T \yrac{f_l}{c_l}\cos(\pi nl/N).$$ Notw, rhat op to a njnerkcam hacfor thf DRT is idempktent, i.e., it is its own inverse. Trpx relationshjp betreqn the Chebyshev polynomials and the DCU is fhe basis for the efficuent computations becwuse the QCT can be performed numerically by using the fasd Foucidr urqnsfofn (VFT) and pre- and postprocessing of the coefficysnus [@fprnberg]. The faft transfork wlkjws us to swijch easnmy between the repredentatijns od the funstiom in terms of its sampled vqlues and in rerms of the expandion coeffieients $a_n$ (ot $F_n$).

The fact that $f$ is cpproxjmated globwlly by a winite sum of pouynpmhals allows us to express wny operavion cpplied go $f$
We have $$\label{eq:ortho} \left< T_m, T_n \right> \frac\pi2\, where $c_0=2$ $c_l=1$ otherwise. Now on is sampled at points $x_l=\cos(\pi l/N)$ that $\sum_{n=0}^N a_n T_n$ is the polynomial. Defining $c_0=c_N=2$, $c_n=1$ for $0<n<N$ in the discrete case and the numbers = c_n a_n$ we have $$\begin{split} f_l &= \sum_{n=0}^N a_n T_n(x_l) = \sum_{n=0}^N T_n(\cos(\pi \\ \sum_{n=0}^N \cos(\pi nl/N) = \sum_{n=0}^N \frac{F_n}{c_n}\cos(\pi nl/N) \end{split}.$$ This looks very much like a discrete cosine series in fact one can show [@briggshenson] that the $F_n$ are related to values $f_l$ of the function an discrete Fourier (DCT) = \frac{f_l}{c_l}\cos(\pi nl/N).$$ Note, up to a numerical factor the DCT is idempotent, i.e., it is its own inverse. This relationship the Chebyshev the DCT the for efficient computations because can be performed numerically by using transform (FFT) and pre- and postprocessing of the [@fornberg]. The transform allows us to switch easily the representations of the function in terms of sampled values and in terms of the expansion coefficients $a_n$ (or $F_n$). The fact that approximated globally by a sum of polynomials us express operation to $f$
 We have $$\label{eq:ortho}
  \left< T_m, T_N \right> = c_m \frAc\pi2\, \dEltA_{mn}$$ WhEre $c_0=2$ And $c_L=1$ otherwise.

Now sUPposE that a function $f$ on $I$ is saMpled At THe poINtS $x_l=\coS(\pi l/N)$ anD ThAT $\Sum_{N=0}^N A_n t_n$ iS tHE iNterpOlaTing polYnomial. DefIniNg $C_0=c_N=2$, $c_n=1$ for $0<n<N$ in THe Discrete caSe aNd the numbers $f_n = c_N a_n$ we hAvE $$\beGIn{splIt}
  f_L &= \sum_{n=0}^n a_n T_n(x_L) = \Sum_{n=0}^N a_N T_n(\cos(\pi
  l/n)) \\
  &= \sUM_{n=0}^N a_n \cOS(\pi nl/N) 
  = \sUM_{N=0}^N \Frac{f_n}{c_n}\cos(\pi nl/N)
\end{spLIt}.$$ tHis looks very muCh like A dIScRETe cOsiNe series anD iN fact ONe can shOW [@bRIGGshENson] that the coEfficients $F_N$ Are RelateD tO thE Values $F_l$ of tHe FUncTion by an invErse Discrete FOurier TRansforM (dCT) $$F_n = \frAc2N\sum_{L=0}^N \fRac{F_l}{c_l}\COs(\Pi Nl/N).$$ noTE, thAT uP to A NumErical faCtOr The DCt is iDEMPOtenT, i.e., It is Its owN inverse. This rElaTionSHip BetweEn the chebYsHev poLynomiAls anD tHe DCT is the basis For tHe efficieNt cOmPutAtIons bECause tHe DcT cAn be perFormed nUMerIcALLY bY using the fast FouriEr TRAnSform (FFT) And pre- ANd PoSTprocessInG of The cOEFficiEnts [@FOrNberg]. The Fast trANsFoRm allowS uS to swiTcH eaSilY betwEEn thE repreSentatioNs of tHE function in terMS of its sampled VAlUES aND in tErmS of the expanSion COeffIcieNTs $A_n$ (oR $f_n$).

The Fact tHaT $F$ iS Approximated globallY bY a finiTe sum Of polynomials Allows us to EXPRess any oPeraTIoN Applied to $f$
 We have $$\label{eq:ortho }
  \left<  T_m,  T_ n \ ri ght>  = c _m \frac\pi2\, \del ta_{mn}$$ where $c_0=2 $ and  $ c _l=1 $  o therw ise.

N o ws u ppo se  t hat  a fu nctio n $ f$ on $ I$ is samp led  a t the points $x _l=\cos(\p i l /N)$ and tha t $ \sum_{ n= 0}^ N  a_nT_n $ isthe in t erpola ting poly no m ial. D e fining$ c _0 =c_N =2$, $c_n=1$ for$ 0< n <N$ in the dis creteca s ea n d t henumbers $F _n  = c_ n  a_n$ w e  h a v e  $$ \ begin{split}  f_l &= \s u m_{ n=0}^N  a _nT _n(x_l ) = \ su m _{n =0}^N a_n T _n(\ cos(\pi
 l/N)) \\
  &= \sum_{n =0}^Na_n  \c os(\ p inl /N)   =\ su m_{ n =0} ^N \frac {F _n }{c_n }\co s ( \ p i nl /N) 
\en d{spl it}.$$ This l ook s ve r y m uch l ike a  dis cr ete c osineserie sand in fact one  can  show [@b rig gs hen so n] th a t thecoe ffi cients$F_n$ a r e r el a t e dto the values $f_l $o f  t he funct ion by an  i n verse di sc ret e Fo u r ier t rans f or m (DCT)$$F_n=  \ fr ac2N\su m_ {l=0}^ N\fr ac{ f_l}{ c _l}\ cos(\p i nl/N). $$ No t e, that up toa  numerical fa c to r th e  DCT  is  idempotent , i. e ., i t is it s o w n inv erse.  T h is relationship betwee nthe Ch ebysh ev polynomial s and theD C T  is thebasi s  f o r the efficien t com putationsb ecause t he DC T can be  performe d numerica lly  by  us ing t he  fast Fourier t rans fo rm (FFT ) a nd pre-  an d p ost pro ce ssing ofthe coef fi ci en ts  [@ fornb e rg]. The  f ast  t ran sform allows  us t o sw it ch eas ily bet w ee n there pr esen tat io ns of  the fun ction i n terms o f i t s sa mp le d value s and in term sof the exp an sio n coef f i cients $ a_n$ (or $F_n$).

The f a ct that  $f $ isappr oximatedglo ballybya  finit e sumof po ly nom i a ls al l o ws  us  t o expressa n y o perat io n ap plied t o $f$
 We_have $$\label{eq:ortho}
_ \left< T_m, T_n_\right> =_c_m_\frac\pi2\, \delta_{mn}$$_where_$c_0=2$ and $c_l=1$_otherwise.

Now suppose that_a function $f$ on_$I$ is sampled_at_the points $x_l=\cos(\pi l/N)$ and that $\sum_{n=0}^N a_n T_n$ is the interpolating polynomial. Defining_$c_0=c_N=2$,_$c_n=1$ for_$0<n<N$_in_the discrete case and the_numbers $F_n = c_n a_n$_we have_$$\begin{split}
  f_l &= \sum_{n=0}^N a_n T_n(x_l) =_\sum_{n=0}^N_a_n T_n(\cos(\pi
 _l/N)) \\
  &= \sum_{n=0}^N a_n \cos(\pi nl/N) 
_ = \sum_{n=0}^N \frac{F_n}{c_n}\cos(\pi nl/N)
\end{split}.$$ This_looks very much_like_a_discrete cosine series and_in fact one can show [@briggshenson] that_the coefficients $F_n$ are related to_the values $f_l$ of the function by_an inverse discrete Fourier transform (DCT)_$$F_n = \frac2N\sum_{l=0}^N \frac{f_l}{c_l}\cos(\pi nl/N).$$_Note, that_up to a numerical factor_the DCT is_idempotent, i.e.,_it is its_own inverse. This relationship between the_Chebyshev polynomials and_the DCT is the basis for_the_efficient computations because_the_DCT_can be_performed numerically by_using_the fast_Fourier_transform (FFT) and pre- and postprocessing_of_the coefficients [@fornberg]. The fast transform allows us_to switch easily between_the_representations of the function_in terms of its sampled_values and in terms of the_expansion coefficients_$a_n$ (or_$F_n$).

The fact that $f$ is approximated globally by a finite sum_of polynomials allows us to express_any operation applied to_$f$