Использование массива в функции

 'процедура осуществляющая решение уравнения методом Гаусса
    Sub Gauss(ByVal P As Integer, ByVal G(,) As Double)
        Dim t, c() As Double
        Dim N1, N2 As Integer
        Dim S As Double
        ReDim c(P + 1)
        ReDim cc(P) 'массив для коэффициентов аппроксимации
        N1 = P + 1
        N2 = P + 2
        For i = 1 To N1 - 1
            t = G(i, i)
            For k = 1 To N2
                G(i, k) = G(i, k) / t
            Next
            For n = i + 1 To N1
                t = G(n, i)
                For j = 1 To N2
                    G(n, j) = G(n, j) - (t * G(i, j))
                Next
            Next
        Next
        c(N1) = G(N1, N1 + 1) / G(N1, N1)
        For i = N1 - 1 To 1 Step -1
            S = 0
            For k = i + 1 To N1
                S = S + G(i, k) * c(k)
            Next
            c(i) = G(i, N1 + 1) - S
        Next
        For i = 0 To P
            cc(i) = Round(c(i + 1)) 'массив с полученными коэффицентами аппроксимации
        Next
    End Sub
    'процедура, решающая задачу аппроксимации
    Sub Approx(ByVal P As Integer, ByVal nu As Integer, ByVal x() As Double, ByVal y() As Double, ByRef d As Double)
        Dim G(P + 1, P + 2) As Double
        For i = 1 To P + 1
            G(i, P + 2) = Summa2(x, y, i - 1, nu)
            For j = 1 To P + 1
                G(i, j) = Summa(x, j + i - 2, nu)
                If i = j Then
                Else
                    G(j, i) = G(i, j)
                End If
            Next
        Next
        Gauss(P, G)
        'вычесление погрешности аппроксимаци
        For j = 0 To nu
            d = d + (cc(0) + cc(1) * x(j) ^ 2 - y(j)) ^ 2
        Next
        d = Round(Sqrt((nu + 1) / d))
    End Sub
    'функция для определения площади криволинейной трапеции' 
    Function q(ByVal x As Single) As Single
        If Simps.RadioButton1.Checked = True Then
            Return Exp(x ^ 2) 'Уравнение для проверки
        Else
            Return(49.1024 + 13.6838 * x - 1.6139 * x ^ 2 + 0.0555 * x ^ 3 - 0.000577 * x ^ 4) ^ 2  'уравнение полученное в результате аппроксимации
        End If
    End Function

'было
Return(49.1024 + 13.6838 * x - 1.6139 * x ^ 2 + 0.0555 * x ^ 3 - 0.000577 * x ^ 4) ^ 2
'стало
Return(cc(0) + cc(1) * x - cc(2) * x ^ 2 + cc(3) * x ^ 3 - cc(4) * x ^ 4) ^ 2