function x=nabisect(fname,a,b,e)
%用途：二分法解非線性f(x)=0；
%格式：x=nabisect(fname,a,b,e)
%     fname為用函數(shù)句柄或內(nèi)嵌函數(shù)表達(dá)的f(x)
%     a, b為區(qū)間端點，e為精度(默認(rèn)值1e-4)，x返回解。
%     程序要求函數(shù)在兩端點值必須異號。

if nargin<4,e=1e-4;end;
fa=feval(fname,a);fb=feval(fname,b);
if fa*fb>0, error('函數(shù)在兩端點值必須異號');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