using python latest version

#-------------------------------------------------------------

def NewtonsMethod(F,FP,x0,epsilon) :

#-------------------------------------------------------------

provides missing code that implements the Newton-Raphson Method

algorithm used in the C source code shown below. Pass lambda functions

actual parameters for F() and FP() in lieu of call-back function pointers

that Python does not support! Try epsilon = 10^-3 or smaller.

//------------------------------------------------------------

// C function that returns the root of F() "close to" x0 using Newton-Raphson Method

//------------------------------------------------------------

double NewtonsMethod(double (*F)(double x),double (*FP)(double x),

double x0,double epsilon)

//------------------------------------------------------------

{

double x = x0

do

{

x = x - F(x)/FP(x);

} while ( fabs(F(x)) >= epsilon );

return( x );

}