function x=nabisect(fname,a,b,e) %fname为内嵌函数
表达式;a,b为区间端点;e为输入定义的精度
if nargin<4,e=1e-4;end; % nargin为输入变量的个数,
若未输入精度变量,则给定默认精度
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
程序如下:
Matlab
Matlab
应
应
用
用
举
举
例
例
(
(
二
二
分
分
法
法
)
)
评论0
最新资源