Cannot find Multiple Roots using Bisection

2 次查看（过去 30 天）

James 2012-3-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/31523-cannot-find-multiple-roots-using-bisection

Function 1:

function    f = example1( x, oldroots)
Bi = 0.5;
             f = 1-x*cot(x)-Bi;;
          if  length( oldroots ) > 0 ,
               f = f / polyval( poly ( oldroots ), x );
          end

Function 2:

function  x  =  findallzeros( fhandle, xguess, maxroots)
  x   = [ ]                ;
    if nargin==1
        xguess =0;
        maxroots = 1;
    end
    if nargin==2
        maxroots  = 1;
    end
    nroots = 0;
    while nroots<maxroots
        if  abs(feval(fhandle,xguess,[ ]))<=1e-16,
            xroot=xguess;
        else
            [xroot,val,flag] = fzero(fhandle,xguess,[ ],x);
            if flag <= 0
                return
            end
        end
        x=[x,xroot];
        nroots=nroots+1;
        xguess=xguess+0.01;
    end