function f = fun_test_newton(x)

f = x.^6-x-1;