先画出图形,确定大概几个根
ezplot('exp(x)-2*x^2-1',[-1,10,0,10])
执行命令得到如图
得知有三个根,编写文件(新建脚本)
function [y,k]=exam2_5(er,n,xa,xb)
syms x xk
x0=xa;x1=xb;
ff=exp(x)-2*x^2-1;
y(1)=xa;y(2)=xb;
k=2;
while abs(x1-x0)>er&k<n
fx1=subs(ff,x,x1);
fx0=subs(ff,x,x0);
x2=x1-fx1*(x1-x0)/(fx1-fx0);
k=k+1;
y(k)=x2;
x0=x1;
x1=x2;
end
执行文件:第1个根
[result,n]=exam2_5(1e-5,100,0,0.2)
结果如图
第2个根
[result,n]=exam2_5(1e-5,100,0.6,0.8)
结果如图
第3个根
[result,n]=exam2_5(1e-5,100,2.8,3)
结果如图