clc
clear
x1=input('Enter the value of x1:');
x2=input('Enter the value of x2:');
for i=1:1:5
t(i)=i;
a(i)=x1;
b(i)=x2;
f1(i)=a(i)^3-9*a(i)+1;
f2(i)=b(i)^3-9*b(i)+1;
x0(i)=a(i)-(f1(i)*(b(i)-a(i)))/(f2(i)-f1(i));
f3(i)=x0(i)^3-9*x0(i)+1;
if i>1
k(i)=((x0(i)-x0(i-1))/x0(i))*100;
m(i)=abs(k(i));
if m(i)<=0.001
break;
end
end
if f1(i)*f3(i)>0
x1=x0(i);
end
if f1(i)*f3(i)<0
x2=x0(i);
end
end
result=[t' a' b' f1' f2' x0' f3' m'];
disp('Iteration x1 x2 f1 fb x0 f3 Error');
disp(result)
root=x0(i);
fprintf('root is %0.4f',root);