Matlab The Bisection Method

MATLAB语言

function y=f(x)

y=f(x); %函数f(t)的表达式

i=0; %二分次数记数

a=a; %求根区间左端

b=b; %求根区间右端

fa=f(a); %计算f(a)的值

fb=f(b); %计算f(b)的值

c=(a+b)/2; %计算区间中点

fc=f(c); %计算区间中点f(c)

while abs(fc)>=ε; %推断f(c)是否为零点

if fa*fc>=0; %推断左側区间是否有根

fa=fc;

a=c;

else fb=fc;

b=c;

end

c=(a+b)/2;

fc=f(c);

i=i+1;

end

fprintf(' %s%.6f %s%d','c,'迭代次数i=',i) %计算结果输出