使用二分法求方程x^5+sin（x）-1=0在[0，1]内的近似根（误差＜10^-5）。

先画出图形，确定大概几个根

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

执行文件, 结果如图: