使用二分法求方程x^5+sin(x)-1=0在[0,1]内的近似根(误差<10^-5)。

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

ezplot('x^5+sin(x)-1',[0,10,0,10])

执行命令得到如图

使用二分法求方程x^5+sin(x)-1=0在[0,1]内的近似根(误差<10^-5)。

编写文件(新建脚本)

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

执行文件, 结果如图:

使用二分法求方程x^5+sin(x)-1=0在[0,1]内的近似根(误差<10^-5)。

 使用二分法求方程x^5+sin(x)-1=0在[0,1]内的近似根(误差<10^-5)。

上一篇:轴角、四元数转换


下一篇:scala高阶函数之匹配功能