rootsl.m

来自「Polynomial Root Finder is a reliable and」· M 代码 · 共 12 行

function [z,e] = rootsl(p)

% function [z,e] = rootsl(p)
%
% computes the zeros z of the polynomial p (descending 
% order). It yields results with higher precision in much
% less time (especially for large degrees) than the 
% built-in function roots. e is an estimate for the maximum
% relative error of all roots. maximum degree is 10000. 
%
% author: m.lang & b.frenzel, 15.6.1993