polyrt0.m

来自「一套科学运算语言 MATLAB 的程序设计有关的函数代码!!!」· M 代码 · 共 17 行

%POLYRT0 用 Newton 搜索算法求解多项式方程一个实根，这是一个
%   M脚本文件，我们当然可以容易地将它改写为一个 M 函数。

%Designed by Prof D Xue (c) 2000
B=polyder(A); ikey=0;
while (ikey==0)
   C=polyval(B,x0); V=polyval(A,x0);
   try
      x1=x0-V/C; 
      if (norm(x1-x0)<eps), ikey=1; break; 
      else, 
         err=norm(x1-x0); disp([x1, C, err]); x0=x1;  
      end
   catch, error(lasterr); end
end
disp(['The solution is ', num2str(x1), ', Error=', num2str(err)])