Work at SourceForge, help us to make it a better place! We have an immediate need for a Support Technician in our San Francisco or Denver office.
% Function not found: FZ_ROOTS
If your polynomials do not have complex coefficients, and if it's OK to use a different algorithm, please try using the IMSL_ZEROPOLY() from the current CVS version. It has the same calling sequence as FZ_ROOTS(). The following wrapper might enable your code to work without changes:
function fz_roots, c, double=d
return, imsl_zeropoly(c, double=d)
GDL implements IMSL_ZEROPOLY() using GSL's gsl_poly_complex_solve() which uses the "balanced-QR reduction of the companion matrix" (as IMSL_ZEROPOLY() does in IDL).