Function DivDif2(x0, x1)

DivDif2 =( f(x1) -f(x0) )/(x1 -x0)

end


Function DivDif3(x0, x1, x2)

DivDif3 =( DivDif2(x1, x2) -DivDif2(x0, x1) )/(x2 -x0)

END