matlab二分法求方程x^3-x-1=0由介值定理，若一个函数在某个区间的两端点异号，则该函数在该区间内有根存在。求根的一种方法—二分法就是基于此原理而来。所谓二分法就是将已给区间两等

Question

matlab二分法求方程x^3-x-1=0由介值定理，若一个函数在某个区间的两端点异号，则该函数在该区间内有根存在。求根的一种方法—二分法就是基于此原理而来。所谓二分法就是将已给区间两等分，取中点的函数值作为近似值，若求得的近似值不满足精度要求，确定二分后新的有根区间，然后检查新的近似值是否达到精度要求，依此类推。用二分法求方程x3-x-1=0在区间[1，1.5]内的一个实根，要求两次近似值之间的误差不超过0.001。 哪位亲知道这道题的解法请告诉我，谢谢啦

Answer 1

这是源代码: 在matlab中保存为:bisection.m functionrtn=bisection(fx,xa,xb,n,delta) %BisectionMethod %Thefirstparameterfxisaexternalfunctionwithrespecttoviablex. %xaistheleftpointoftheinitialinterval %xbistherightpointoftheinitialinterval %nisthenumberofiterations. x=xa;fa=eval(fx); x=xb;fb=eval(fx); disp('[nxaxbxcfc]'); fori=1:n xc=(xa+xb)/2;x=xc;fc=eval(fx); X=[i,xa,xb,xc,fc]; disp(X), iffc*faf='x^3-x-1'; >>bisection(f,1,1.5,20,10^(-3)) [nxaxbxcfc] 1.00001.00001.50001.2500-0.2969 2.00001.25001.50001.37500.2246 3.00001.25001.37501.3125-0.0515 4.00001.31251.37501.34380.0826 5.00001.31251.34381.32810.0146 6.00001.31251.32811.3203-0.0187 7.00001.32031.32811.3242-0.0021 8.00001.32421.32811.32620.0062 9.00001.32421.32621.32520.0020 从结果可以看出, 这个解为:1.3262