二分法,不动点迭代法,艾特肯加速迭代法,牛顿切线法的matlab程序及举例.docVIP

§ 2. 1. 1二分法的MATLAB主程序 function [k, x, wuca, yx]=erfen (a，b，abtol) a(l)=a; b(l)=b; ya二fun(a(l)) ; yb二fun(b(1)) ; %程序中调用的 fun. m为函数 if ya* yb〉0, dispC注意：ya*yb〉0,请重新调整区间端点 a和b.’)，return end maxl=-l+ceil ((log (b~a)- log (abtol))/ log(2)); % ceil是向+OO方向取整 for k=l: maxl+1 a;ya=fun(a); b;yb=fun(b); x=(a+b)/2; yx=fun (x); wuca=abs(b~a)/2; k=k_l; [k，a，b，x，wuca, ya, yb, yx] if yx==0 a=x; b=x; elseif yb氺yx0 b=x;yb=yx; else a=x; ya=yx; end if b~a abtol ， return, end end k=maxl; x; wuca; yx=fun(x); §2. 1.2不动点迭代法的MATLAB主程序 function [k, piancha, xdpiancha, xk, yk]=diedai2 (xO, tol, ddm ax) x (l)=xO; for i=