= Solution
For the intended collocation weights, interpolation makes the associated <quadrature rule> exact for every polynomial of degree at most one, so the method has order at least two. It has order at least three precisely when the node polynomial is orthogonal to constants:
$$
0=\int_0^1(t-\alpha)(t-1)\,dt
=\frac\alpha2-\frac16.
$$
Thus
$$
\boxed{p=3\text{ when }\alpha=\frac13,
\qquad p=2\text{ for every other }\alpha\ne1.}
$$
At $\alpha=1/3$ this is the two-stage <Radau IIA method>. Since one node is fixed at the endpoint, no value of $\alpha$ makes the two-node quadrature exact through degree three, so order four cannot occur.
For completeness, the tableau exactly as printed has order zero when $\alpha\ne1/2$, since $b^T\mathbf1\ne1$. At $\alpha=1/2$ the two signs coincide because the second weight vanishes, and the resulting method has order two.
Back to article page
