function p = newton(f,df,p0,TOL,Nmax)
 %write your function here

end