使用二分法求方程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
执行文件, 结果如图:
标签:10,subs,xk,fx,二分法,ff,sin 来源: https://blog.csdn.net/weixin_62155702/article/details/121609978