先画出图形,确定大概几个根
ezplot('x^5+sin(x)-1',[0,10,0,10])
执行命令得到如图
编写文件(新建脚本)
function y =exam2_241(m,n,er)
syms x xk
a=m;b=n;k=0;
ff=x^5+sin(x)-1;
while b-a>er
xk=(a+b)/2;
fx=subs(ff,x,xk);
fa=subs(ff,x,a);
k=k+1;
if fx==0
y(k)=xk;
break;
elseif fa*fx<0
b=xk;
else
a=xk;
end
y(k)=xk;
end
plot(y,'.-');
grid on
执行文件, 结果如图: