Skrypt realizujący tablicowanie funkcji w języku Visual Basic w określonym przedziale.

Cześć,
Mam do przygotowanie skrypt w Visual Basicu realizujący tablicowanie funkcji w określonym przedziale -0,2<x<0,2 z wykorzystaniem jej wartości dokładnych i przybliżonych.
Wzór funkcji:
(1+5x)^1/2=1+1/2(5x)-1/(2^22!)(5x)^2+13/(2^33!)(5x)^3-135/(2^44!)(5x)^4+1357/(2^5*5!)...
Próbowałem coś napisać, jednak dość słabo mi idzie.
Poniżej załączam to, co udało mi się napisać..
Jest ktoś w stanie pomóc?
Z góry dziękuję.

Option Explicit
Dim a As Double, b As Double
Dim n As Integer, i As Integer, j As Integer
Dim eps As Double
Dim x As Double, s As Double, w As Double
Dim f As Double
Dim dx As Double
Dim q As Double
Sub szereg()
a = CDbl(Cells(1, 1).Value)
b = CDbl(Cells(1, 2).Value)
n = CInt(Cells(1, 3).Value)
eps = CDbl(Cells(1, 4).Value)
dx = (b - a) / CDbl(n)
For i = 0 To n Step 1
x = a + CDbl(i) * dx
f = (1 + 5 * x) ^ (1 / 2)
s = 1
w = 1 / 2 * (5 * x)
j = 2
q = 2

Do
s = s + w
w = (-1) * w * q * x / (2* j)
j = j + 1
q = 1 + (3 * (j - 2))

Loop While Abs(w) >= eps
Cells(4 + i, 1).Value = i + 1
Cells(4 + i, 2).Value = x
Cells(4 + i, 3).Value = f
Cells(4 + i, 4).Value = s
Cells(4 + i, 5).Value = (f - s) / (f * 100)
Next i
End Sub