clear;
clc;

f = @(x) (x^3 - x - 1);
fp = @(x) (imag(f(x + 1i*10^-8)) / 10^-8);

xl = 1;
xu = 2;
acc = 0.001;

x = (xl + xu) / 2;
itr = 0;

while abs(f(x)) > acc
    itr = itr + 1;
    x = x - f(x)/fp(x);
end

disp("Result: " + x);
disp("Iteration: " + itr);