Radau and Lobatto Quadrature

For practical computations, the following weighting formxxla was derived by Villadsen and Michelsen (1978) 

Similarly, when the boundary point at x = 0 is added to the N interior interpolation points, the interior points must be chosen as roots of the following Nth degree polynomial 

5 LINEAR BOUNDARY VALUE PROBLEM—DIRICHLET BOUNDARY CONDITION 

The diffusion-reaction problem for slab catalyst particles is a classic problem used to illustrate the orthogonal collocation method. We consider this problem next. 

---