\centerline{\bf ORTHOGONAL POLYNOMIALS} \bigskip Let $w:(a,b)\to E$ be continuous, nonnegative, and integrable, with at most isolated zeros on $(a,b)$. Then $$\la f,g\ra=\int_a^b f(x)g(x)w(x)\,dx$$ defines an inner product on the vector space $V=\left\{f\in C(a,b)\mid \int_a^b [f(x)]^2 w(x)dx<\infty\right\}$. Given distinct points $x_1$, $\ldots$, $x_M\in[a,b]$ and positive constants $w_1$, $\ldots$, $w_M$, $$\la f,g\ra= \sum_{i=1}^M f(x_i)g(x_i)w_i$$ defines a (discrete) inner product on any vector space of functions $V$ with the property that $f\in V$, $f(x_1)=\cdots= f(x_M)=0 \Rightarrow f\equiv0$. In the following, $\la f,g\ra$ can be either a continuous (integral) or a discrete (sum) inner product. \noindent{\bf Definition}. $\{Q_i\}_{i=0}^m$ ($m+1-1\cr (a,b)=(-\infty,\infty)& [a,b)=[0,\infty)\cr H_{k+1}(x)=2x\,H_k(x)-2k\,H_{k-1}(x)& L_{k+1}^\alpha(x)={-1\over k+1} (x-2k-\alpha-1)L_k^\alpha(x) \cr 1& - {k+\alpha\over k+1}L_{k-1}^\alpha(x) \cr 2x\cr 4x^2-2\cr 8x^3-12x\cr 16x^4-48x^2+12\cr}}$$ \hrule\medskip Given $P_S(f)=R(x)=\sum_{k=0}^n \beta_k Q_k(x)$, where $\beta_k= \la f,Q_k\ra/\la Q_k,Q_k\ra$, and $\{Q_i\}_{i=0}^n$ are orthogonal, $R(x)$ can be efficiently evaluated using the three term recurrence relation from above: $$\vbox{\tt \hbox{$d_{n+2}\: d_{n+1}\: 0$;}\hbox{for $k\:n$ step $-1$ until $0$ do} \hbox{\qquad$d_k\:\beta_k+A_k\bigl(x-B_k\bigr)d_{k+1} - C_{k+1}d_{k+2}$; \qquad( $R(x)=d_0Q_0(x)$ )}}$$ \vfil\eject