二分法求解单变量非线性方程及其应用与实现(2)
2013-07-15 01:12
导读：0 1 2 3 4 5 61.0 1.25 1.3125 1.32031.5 1.375 1.3438 1.32811.25 1.375 1.3125 1.3438 1.3281 1.3203 1.3242－ + － + + － － 6. 基本二分法的matlab实现与C语言实现 6.1 %二分法的算法及
0
1
2
3
4
5
61.0
1.25
1.3125
1.32031.5
1.375
1.3438
1.32811.25
1.375
1.3125
1.3438
1.3281
1.3203
1.3242－
+
－
+
+
－
－
6. 基本二分法的matlab实现与C语言实现
6.1 %二分法的算法及MATLAB实现
function [c, err, yc] = bisect(f, a, b, delta)
% f 是所要求解的函数
% a 和 b 分别是有根区间的左右限
% delta 是允许的误差界
% c 为所求的近似解
% yc 为函数 f 在 c 上的值
% err 是 c 的误差估计
if nargin < 4
delta = 1e -5;
end
ya = feval (’f’, a);
yb = feval (’f’, b);
if yb == 0, c = b, return
end
if ya * yb > 0
disp(’(a, b)不是有根区间’);
return
