求解非线性方程的实根matlab

1.求解方程P(k+1)=cos(sin(P(k)))的解
代码可供参考
clear;
clc;

p0=0.3;
tol=0.00000001;
max=2000;
P(1)=p0;
for k=1:max
P(k+1)=cos(sin(P(k)));
err(k+1)=abs(P(k+1)-P(k));

if err(k+1)<tol
    break;
end
p=P(k+1);

end
fprintf(“此非线性方程的解为；%f”,p);

2.求round((log(b-a)-log(tol))/log(2))的解

clc;
tol=0.00000001;
a=-1;
b=2;
ya=feval(@f,a);
yb=feval(@f,b);
max=round((log(b-a)-log(tol))/log(2));

if yayb>0
fprintf(“试位法不适用”);
else
for k=1:max
c=(a+b)/2;
yc=feval(@f,c);
if yc==0
b=c;
a=c;
elseif ycyb>0
b=c;
else
a=c;
end
err(k)=b-a;
if err(k)<tol
break
end
end
fprintf(“结果为；%f”,(a+b)/2);
end
function [y]=f(x)
y=sin(x);
end