.(8)令err=|b-a|.(9)判断err<delta是否成立,若成立,输出结果。若不成立,令k=k+1,执行(4)。流程图:startInputf,a,b,delta,max1,,k=k+1Yf(a)*f(b)>0NoutputYk>max1Ndx=yb*(b-a)/(b-a)c=b-xyc=feval(f,c)Yyc=0NYNyc*yb>0b=cyb=yca=cya=ycerr=|b-a|Nerr<deltaYoutputend当输入[c,err,yc]=regula(@f,0.15,0.16,1e-10,100)输出c=0.157539310266284err=-1.257285475730896e-008yc=2.775557561562891e-0173.(a)算法:(1)输入f,a,b,N,令k=1.(2)判断k>N+1是否成立,若成立,输出并执行步骤(3);若不成立,令p(k)=feval(f,b*(k-1)/(a*N)),k=k+1,再执行步骤(2)。(3)令m=p(1),k=2.(4)判断k>N+1是否成立,若成立,输出并执行步骤(6);若不成立,执行步骤(5)。(5)判断m>p(k)是否成立,若成立,m=p(k),执行步骤(4);若不成立,执行步骤(4)。(6)令n=p(1),k=2.(7)判断k>N+1是否成立,若成立,输出;若不成立,执行步骤(8)。(8)判断n>p(k)是否成立,若成立,n=p(k),输出;若不成立,执行步骤(7)。流程图:startInputf,a,b,Nk=1k=k+1Yk>N+1Noutputp(k)=feval(f,b*(k-1)/(a*N))m=p(1),k=2k=k+1Yk>N+1NNoutputm>p(k)Ym=p(k)n=p(1),k=2k=k+1Yk>N+1NNn>p(k)outputYendn=p(k)(b)
全文预览
