Несобственный интеграл (численные методы)

Function f(ByVal m As Double, ByVal d As Double, ByVal x As Double) As Double
        f = 1 / (2 * Sqrt(PI)) * Exp(-(x - m + d) ^ 2 / (2 * 0.5 ^ 2))
        Return f
    End Function
 Sub Integral(ByVal m As Double, ByVal d As Double, ByVal E As Double, ByRef t As String, ByRef s As Double)
        Dim a, b, s1, h, x As Double
        Dim n, c As Integer
        a = 0 'нижний пределел интегрирования
        b = PI / 2 'верхний предел интегрирования
        t = "" 'строка для вывода хода решения
        n = 2 'кол-во разбиений
        h = (b - a) / n ' шаг
        s = (f(m, d, a) + 4 * f(m, d, (a + b) / 2) + f(m, d, b)) * (h / 3)
        t = (Format(h, "0.0000") & Space(5) & Format(n, "000") & Space(5) & Format(s1, "0.00000000000") _
    & Space(5) & Format(s, "0.00000000000")) & vbNewLine
        Do
            n = 2 * n
            h = (b - a) / n
            s1 = s
            c = 4
            x = a
            s = f(m, d, a) + f(m, d, b)
            For i = 1 To n - 1
                x = x + h
                s = s + c * f(m, d, x)
                c = 6 - c
            Next
            s = s * (h / 3)
            t = t & (Format(h, "0.00000") & Space(5) & Format(n, "000") & Space(5) & Format(s1, "0.00000000000") _
    & Space(5) & Format(s, "0.00000000000")) & vbNewLine
        Loop Until Abs((s - s1) / 15) < E ' условие окончания цикла DoLoop по оценке погрешности методом двойного просчета 
    End Sub