数值分析的个简单算法实现.docVIP

数值分析的几个简单算法实现??Matlab6.5 ? ? 2006-11-24 07:42:12 大 中 小 %埃特金算法Aitken.m function x=Aitken(fname,x0,N,e) k=1;e=1e-5;x2=x0+2*e; while abs(x2-x0)ekN x1=feval(x0);x2=feval(x1); x2=x2-(x2-x1)^2/(x2-2*x1+x0); k=k+1; x2=x0;end; disp(x2); if k==N,warning(迭代失败);end; end ? ? ? %迭代法diedai.m文件 function x=diedai(fname,x0,e,N) if nargin5;N500;end; if nargin4;e=1e-4;end; x=x0;x0=x+2*e;e=1e-5;k=1; while abs(x0-x)ekN; x=feval(fname,x0); k=k+1; x0=x; disp(x); end if k==N,warning(迭代失败);end ? ? %二分法erfenfa.m function x=erfenfa(fname,a,b,e) y0=feval(fname,a);e=1e-5; while (b-a)e