Matematicheskie_metody / Мат. методы. 2 семестр(1 модуль) / Метод Эйткена

Sub иипэ()
Const n = 6, e = 0.0001, x = 1.5, i0 = n \ 2 + 1
Dim i As Single, s As Single, k As Integer, p(-1 To (n + 1) * (n + 2) / 2 - 1) As Double, dx(n) As Double, xz(n) As Double, j As Integer
For i = 0 To n
s = (i / 2 - i \ 2) * 4 - 1
k = i * (i + 1) / 2
p(k) = Cells(i0 + s * (i + 1) \ 2, 2)
xz(i) = Cells(i0 + s * (i + 1) \ 2, 1)
dx(i) = x - xz(i)
For j = 1 To i
p(k + j) = (dx(j - 1) * p(k + j - 1) - dx(i) * p((j + 1) * j / 2 - 1)) / (xz(i) - xz(j - 1))
Next j
If Abs(p(k + i) - p(k - 1)) < e Then Exit For
Next i
Debug.Print p(k - 1), i - 1
End Sub