数值分析程序纪要(MATLAB)
2018-03-29 本文已影响0人
叫酿么
一、求方程的根
(1)二分法
数值分析程序纪要(MATLAB)
function xc=bisect(f,a,b,tol)
if sign(f(a))*sign(f(b)) ＞＝0
error('f(a)f(b)<0 not satisfied')
end
fa=f(a);
fb=f(b);
while (b-a)/2＞tol
c=(a+b)/2;
fc=f(c);
if fc==0
break
end
if sign(fc)*sign(fa)<0
b=c;fb=fc;
else
a=c;fa=fc;
end
end
xc=(a+b)/2;
定义函数f=@(x) 表达式，
xc=bisect(f，a，b，tool)
＊注记:
