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

Sub ИФН1()
Const n = 6, e = 0.0001, x = 1.5
Dim i As Integer, x0 As Double, x1 As Double, t As Double, dy((n + 1) * (n + 2) / 2 - 1) As Double, P As Double, c As Double, z As Double, k As Integer, j As Integer, P1 As Double
x0 = Cells(1, 1)
x1 = Cells(2, 1)
c = 1
z = 1
t = (x - x0) / (x1 - x0)
dy(0) = Cells(1, 2)
P = dy(i)
For i = 1 To n
k = i * (i + 1) / 2
dy(k) = Cells(i + 1, 2)
For j = 1 To i
dy(k + j) = ((dy(k + j - 1) - dy(k + j - i - 1)))
Next j
c = c * (t - (i - 1))
z = z * i
P1 = c / z * dy(k + i)
If Abs(P1) < e Then Exit For
P = P + P1
Next i
Debug.Print P, i - 1, Abs(P1)
End Sub