Julia analog to R's poly function, orthogonal polynomials by QR decomposition

[stevengj]

@Daniel_Johansson

@danielwe @stevengj do any of you have some reference (introductory if possible) for this kind of topics (i.e. obtaining orthogonal polynomials provided a collection of knots), or that talk about reparametrization of a function basis using design matrix in general?

Trefethen and Bau lecture 37 talks about the connection of orthogonal polynomials to Lanczos. The only difference between e.g. Legendre polynomials and your case (discrete collection of points) is the inner product.