Помогите запилить программу по методу Симпсона

Module Module1
    Public a, b, E, T, h, x11, x12, x21, x22, n1, I11, I12, I21, I22, C As Double
    Sub Main()
        a = CDbl(Console.ReadLine())
        b = CDbl(Console.ReadLine())
        E = CDbl(Console.ReadLine())
        T = CDbl(Console.ReadLine())
        Dim I1, I2, n As Double
        While E >= Math.Abs(I1 - I2) / 15
            n = 3
            h = (b - a) / n
            I11 = 0
            x11 = a
            x12 = a + h
            For i = 1 To n Step 1
                x11 = x11 + h
                I11 = I11 + F(x11)
                x12 = x12 + h
                I12 = I12 + F(x12)
            Next i
            I1 = h * (F(a) + F(b) + 2 * I1) / 2
            n1 = 2 * n
            h = (b - a) / n1
            I2 = 0
            x21 = a
            x22 = a + h
            For i = 1 To n1 Step 1
                x21 = x21 + h
                x22 = x22 + h
                I21 = I21 + F(x21)
                I22 = I22 + F(x22)
            Next i
            I2 = h * (F(a) + F(b) + 2 * I21 + 4 * I22) / 6
        End While
        Console.WriteLine(I1)
        Console.ReadKey()
    End Sub
    Function F(ByVal x) As Single
        F = x * Math.Exp(-T * x ^ 2) * Math.Sin(T ^ 3 * X ^ 2)
    End Function
End Module