Matlab编写程序用二分法求解非线性方程的根

Matlab编写程序用二分法求解非线性方程的根

用二分法求方程x*x-x-1=0的正根，要求精确到小数点后四位。（matlab）
l1 计算公式
f(ak)*f(bk)<0;
bk-ak=1/2k-1*(b-a);
a1<=a2<=… <=ak<=…,b1<=b2<=…<=bk<=…。
l2 算法分析
设f(x)∈C[a,b],且f(a)f(b)<0，在[a,b]间寻找实根，记[a1,b1]=[a,b]，取x1=(a1+b1)/2,若f(x1)=0,则x1是f(x)=0的根，f(x1)f(a1)>0,则a2=x1,b2=b1,否则a2=a1,b2=x1。得到[a2,b2]满足：f(a2)f(b2)<0,b2-a2=(b1-a1)/2=(b-a)/2,a2>=a1,b2<=b1。以[a2,b2]取代[a1,b1]，继续以上过程，直到精度达到要求为止。
l3 源程序
function f1=fun(x);
f1=x-cos(x);
function [x,k]=erfen(a,b,s) %a,b为根区间，s为精度
a=0;b=1;s=1e-4;k=0;
while abs(a-b)>s
x=(a+b)/2;
if fun(a)*fun(x)<0
b=x;
else
a=x;
end
k=k+1;
end
x=(a+b)/2 %x为方程的解
k % k为计算次数
实验结果讨论和分析
本题使用二分法得到的x=0.7391，满足基本要求，题目要求精确到小数点后四位，告诉了本题二分法得应达到得精确度；计算次数为14，二分法收敛性很好，收敛速度不快。