function [x, itr] = newtonraphson(f, f1, a, b, acc)
x = (a + b) / 2;

itr = 0;
disp(["itr", "x", "y"])
while true
    itr = itr + 1;
    
    x = x - f(x)/f1(x);
    y = f(x);
    disp([itr, x, y])

    if (abs(y) < acc)
        break;
    end
end