Das Interpolationsverfahren nach Newton macht keine Vorgaben zur Reihenfolge der Stützpunkte. Wenn Du also die Terme mit \(b\) vermeiden willst, so stelle den Punkt \((1|\, b)\) nach hinten. Die übrigen drei Punkte liegen auf einer Gerade, was dazu führt, dass das \(a_2\) zu \(0\) wird; das vereinfachte das ganze weiter. Die Rechnung sollte dann so aussehen:
$$ \begin{array}{cccc|ccc} &&& x_i & y_i \\ \hline &&& 0 & \colorbox{#FFFF00}{2} \\ &&3 & 3 & 5 & \colorbox{#FFFF00}{$\frac{5-2}{3}$} \\ &4 &1 & 4 & 6 & \frac{6-5}1=1 & \colorbox{#FFFF00}{0} \\ 1 &-2 &-3 & 1 & b & \frac{b-6}{-3} & \frac{6-b-3}{3 \cdot (-2)}=\frac{b-3}{6} & \colorbox{#FFFF00}{$\frac{b-3}{6}$}\end{array}$$
und das Polynom ist dann:
$$\begin{aligned} y(x) &=2 + 1(x-0) + 0 + \frac{b-3}{6}(x-0)(x-3)(x-4) \\ &= \frac{b-3}{6}x^3 - \frac{7(b-3)}{6}x^2 + (2(b-3)+1)x + 2 \end{aligned}$$
anbei drei Graphen mit \(b=\{3, \, 6, \, 9\}\):
~plot~ {0|2};{3|5};{4|6};[[-2|6|-10|15]];x+2;x^3/2-7x^2/2+7x+2;{1|3};{1|6};{1|9};x^3-7x^2+13x+2 ~plot~
Gruß Werner