nabisect.m

Matlab数值解算法实现代码

function x=nabisect(fname,a,b,e)
%用途：二分法解非线性f(x)=0；
%格式：x=nabisect(fname,a,b,e)
%     fname为用函数句柄或内嵌函数表达的f(x)
%     a, b为区间端点，e为精度(默认值1e-4)，x返回解。
%     程序要求函数在两端点值必须异号。

if nargin<4,e=1e-4;end;
fa=feval(fname,a);fb=feval(fname,b);
if fa*fb>0, error('函数在两端点值必须异号');end
x=(a+b)/2
while (b-a)>(2*e),
   fx=feval(fname,x);
   if fa*fx<0,b=x;fb=fx;else a=x;fa=fx;end
   x=(a+b)/2
end