function estimate=NewtonsMethodFunction(f,fp,x0,maxabserror)

x1=x0-f(x0)/fp(x0);
while abs(x0-x1)>maxabserror
    x0=x1;
    x1=x0-f(x0)/fp(x0);
end
estimate=x1;