MATLAB实现Steffensen迭代法.doc

syms x;f=input(输入函数形式 fx=);a=input(输入初始近似值 a=);tic;x0=subs(f,x,a)+a;x1=subs(f,x,x0)+x0;x2=a-(x0-a)2)/(x1-2*x0+a);k=1;x2=double(x2);fprintf(k=%d x=%fn,k,x2);while abs(x2-a)0.0001a=x2;x0=subs(f,x,a)+a;x1=subs(f,x,x0)+x0;x2=a-(x0-a)2)/(x1-2*x0+a);k=k+1;x2=double(x2);fprintf(k=%d x=%fn,k,x2);endtoc;