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

end