' преобразуем исходную функцию к виду: 
'   f(x) = (x - 3)*(x^2 + 4) + 3
' формула для вычисления следующего приближения x:
'   [xi] = 3 - 3 / ([xi-1]^2 + 4), где
'   [xi] - следующее приближение для x
'   [xi-1] - приближение, полученное на предидущем шаге


' код:

Sub iteration_method()

Const epsilon! =  0 . 001     ' точность вычислений
Const MAXITER% =  1000      ' максимальное кол-во итераций
Const X0! =  2 !            ' нулевое приближение

Dim x!, x2!
Dim count%

x2 = X0

Debug.Print "Iter#", "xi", "delta"

Do Until Abs(x2 - x) < epsilon Or count >  1000 
    x = x2
    x2 =  3  -  3  / (x * x +  4 )
    count = count +  1 
    Debug.Print Format(count, "000"), x2, Format(Abs(x2 - x), "0.0000")
Loop
if Abs(x2 - x) > epsilon Then _
Debug.Print "(!) требуемая точность не достигнута. Ищите другую формулу для приближений."

End Sub