Module Module1

    'changed EXAMPLE FROM ALGLIB PACKAGE

    Public Sub function_cx_1_func(c As Double(), x As Double(,), ByRef func As Double, obj As Object)
        '
        '
        func = (c(0) * x(0, 0)) + (c(0) * x(0, 1)) + (c(0) * x(0, 2)) + (c(1) * x(0, 1)) + (c(1) * x(1, 1)) + (c(1) * x(2, 1))

    End Sub

    Public Sub Main()
        '
        ' 
        '
        '  Gradient is estimated using combination of numerical differences
        '  and secant updates. diffstep variable stores differentiation step
        '  (we have to tell algorithm what step to use).
        '
        Dim x(,) As Double = New Double(,) {{0.0, 1.0}, {3.5, 0.0}, {1.0, 2.0}}
        Dim y() As Double = New Double() {1, 4, 7}
        Dim c() As Double = New Double() {0.6, 0.3}
        Dim epsf As Double = 0
        Dim epsx As Double = 0.000001
        Dim maxits As Integer = 0
        Dim info As Integer
        Dim state As lsfitstate = New XAlglib.lsfitstate() ' initializer can be dropped, but compiler will issue warning
        Dim rep As lsfitreport = New XAlglib.lsfitreport() ' initializer can be dropped, but compiler will issue warning
        Dim diffstep As Double = 0.0001

        '
        '  Fitting without weights
        '
        XAlglib.lsfitcreatef(x, y, c, diffstep, state)


        xalglib.lsfitsetcond(state, epsf, epsx, maxits)
        XAlglib.lsfitfit(state, AddressOf function_cx_1_func, Nothing, Nothing)
        xalglib.lsfitresults(state, info, c, rep)
        System.Console.WriteLine("{0}", info) ' EXPECTED: 2
        System.Console.WriteLine("{0}", alglib.ap.format(c, 1)) ' EXPECTED: 2,1


       
        System.Console.ReadLine()
        Environment.Exit(0)
    End Sub
End Module