kv

-

kernel handle returned by StartMaple

s

-

Maple object object

fn

-

Maple FUNCTION object

n

-

number of arguments

args

-

array of n hardware floats

Function MyNewton(ByVal kv As Long) As Double

    Dim i As Long

    Dim guess(2) As Double

    Dim tolerance As Double

    Dim f, fprime, x, initguess As Long

    ' define a function in x to solve, and find it's derivative

    f = EvalMapleStatement(kv, "(x) -> x^3+2*x^2-9*x+2;")

    MapleAssign kv, ToMapleName(kv, "f", True), f

    fprime = EvalMapleStatement(kv, "unapply(diff(f(x),x),x);")

    x = ToMapleName(kv, "x", True)

    guess(1) = 1.1

    tolerance = 0.0001

    ' find a root (abort after 500 iterations)

    For i = 0 To 500

        'test current guess

        If Abs(EvalhfMapleProc(kv, f, 1, guess(0))) <= tolerance Then

            i = 501 'break

        Else

            'refine guess

            guess(1) = guess(1) - EvalhfMapleProc(kv, f, 1, guess(0)) / _

                        EvalhfMapleProc(kv, fprime, 1, guess(0))

        End If

    Next i

    If i = 500 Then

        MyNewton = -999999  'couldn't find root

    Else

        MyNewton = guess(1)

    End If

End Function