二分法求解非线性方程的近似解:M-文件 funn1+fun
程序:
function [x0 fval]=funn1(fun,a,b,D)
if nargin<3
error('a and b should redefinated');
else if nargin==3
D=1e-6;
end
end
fa=feval(fun,a);
fb=feval(fun,b);
if fa*fb>0
error('this equation has no answer');
else if fa==0
m=a;
fm=fa;
else if fb==0
m=b;
fm=fb;
else
while abs(b-a)>D
m=(a+b)/2;
1 / 3