Построить интерполяционный полином Ньютона для функции

Построим интерполяционный полином Ньютона второй степени
Составим таблицу, содержащую значения x , y и конечные разности:
x y ∆y
   ∆y2
1,3 3,6693 0,8124 0,1798
1,5 4,4817 0,9922 0,2198
1,7 5,4739 1,212
1,9 6,6859
Соответствующий полином имеет вид
Px=a0+a1x-1,3+a2x-1,3(x-1,5)
Используя конечные разности из этой таблицы и, учитывая, что в нашем случае h=1, вычислим коэффициенты этого полинома:
a0=y0=3,6693
a1=∆y0h=0,81240,2=4,062
a2=∆2y02h2=0,17982∙0,22=2,2475
Искомый полином будет
Px=3,6693+4,062x-1,3+2,2475x-1,3(x-1,5)
После преобразований получим
Px=2,2475x2-2,231x+2,7713
P1,4=2,2475∙1,42-2,231∙1,4+2,7713=4,053