clc,clear

syms x
%% 第一小题
phi = x^3-1;
f = x^3-x-1;
a = 1; % 区间左端点
b = 2; % 区间右端点
x = 1.5; % 初始值

%% 第二小题
% phi = 2*log(x)+log(3);
% f = 3*x^2-exp(x);
% a = 3; % 区间左端点
% b = 4; % 区间右端点
% x = 3.5; % 初始值


n = 1; % 迭代次数
tol = 1e-4; % 精度


fprintf('第 %d 步\t迭代值为：%f, 函数值为：%f\n',n,x,subs(f,x));

hold on;
plot(a:0.01:b,subs(f,a:0.01:b),'b','Linewidth',2);
plot(x,subs(f,x),'ro','Linewidth',2);


while abs( subs(f,x) )>tol
     
    n = n+1;
    
    x = x - (subs(phi,x)-x)^2/(subs(phi,subs(phi,x)) - 2*subs(phi,x) + x);
    
    
    fprintf('第 %d 步\t迭代值为：%f, 函数值为：%f\n',n,x,subs(f,x));
    plot(x,subs(f,x),'ro','Linewidth',2);
    pause(0.5);
     
end