Let $x_1, x_2 \in \mathbb R$ and let $(a_0, a_1, \dots, a_n), (b_0, b_1, \dots, b_n)$ be $(n+1)$-tuples of real numbers. Show that there is a unique polynomial of degree at most $2n+1$ such that $$q^{[k]}(x_1)=a_k,$$$$q^{[k]}(x_2)=b_k$$ for $k=0, \dots, n$.
Any hints on how to get started with this exercise?
As a hint, I would use the following fact:
In your case, the functionals $\cdot^{[0]}(x_1), \ldots, \cdot^{[n]}(x_1), \cdot^{[0]}(x_2), \ldots, \cdot^{[n]}(x_2)$ on $P_{2n+1}(\mathbb{R})$ have count exactly equal to the dimension of the dual space. Therefore, they form a basis of the dual space (equivalent to the desired conclusion) if and only if they span the dual space.
This reduces the original problem to proving that if $p^{[0]}(x_1) = \cdots = p^{[n]}(x_1) = p^{[0]}(x_2) = \cdots = p^{[n]}(x_2) = 0$ for $p \in P_{2n+1}(\mathbb{R})$, then $p=0$.