Ton polynôme P0 ne convient pas. Sa dérivée en 0 n'est pas nulle.
J'explique le système de congruences :
Les conditions en 0 peuvent se traduire par le fait que

est congru à

modulo

, et les conditions en 1 par le fait que

est congru à

modulo
^2)
.
Le théorème des restes chinois, avec l'égalité de Bezout
^2(2X+1)-X^2(2X-3))
nous dit que les solutions sont les polynômes

congrus à un certain polynôme que je te laisse calculer modulo
^2)
.
Ce sont des calculs un peu lourds. Il est moins fatigant de chercher l'unique polynôme de degré 3 satisfaisant les conditions (les autres lui sont congrus modulo
^2)
. Les conditions en 0 imposent au polynôme d'être de la forme

et les conditions en 1 permettent de déterminer

et

.
On peut aussi signaler la méthode d'interpolation de Hermite par les différences divisées :

qui donne
\times X^2 + \mathbf{3}\times X^2(X-1))
.