Помогите разработать программу нахождения значения определённого интеграла - VB

 Function f(ByVal x As Double) As Double
        'подынтегральная функция
        f = -(x ^ 3) + 4
    End Function
    'a,b - границы участка интегрирования, eps - точность, h- начао=льный шаг интегрирования
    Sub SM(ByVal a As Double, ByVal b As Double, ByVal eps As Double, ByVal h As Double)
        'n-число разбиений отрезка интегрирования
        Dim n, c As Integer
        Dim s, s1, x As Double
        s = (f(a) + f(b) + 4 * f((a + b) / 2)) * Abs(a - b) / 3
        Do
            s1 = s
            s = f(a) + f(b)
            x = a
            n = (b - a) / h
            c = 4
            For i = 1 To n - 1
                x = x + h
                s = s + c * f(x)
                c = 6 - c
            Next
            s = s * h / 3
            'вывод промежуточных данных
            ListBox1.Items.Add(CStr(n))
            ListBox2.Items.Add(CStr(s))
            ListBox3.Items.Add(CStr(h))
            h = h / 2
        Loop Until (Abs(s - s1) / 15 < eps) 'проверка точности оп правилу Рунге
        TextBox5.Text = CStr(s) 'Вывод точного значения
    End Sub