y = @(x) x^3-(0.165*x^2)+(3.993*10^-4);
x0 = input("Enter first guess: " );
x1 = input("Enter second guess: " );
max_iterations = input("Enter max iterations: " );
tolerance = input("Enter tolerance: " );

iterations = 0;

while iterations < max_iterations
    iterations = iterations + 1;

    x2=x1-(y(x1)*(x1-x0))/(y(x1)-y(x0))

    relative_approximate_error =abs(x2-x1)/x2;
    if relative_approximate_error < tolerance
    break;
    endif

x0=x1;
x1=x2; 
    
endwhile
fprintf("No of iterations: %d \n", iterations);
fprintf("Root: %d \n",x2);