GENERIC INTERFACE RootApproximation(R, RRt, C, CRt);

Arithmetic for Modula-3, see doc for details

Abstract: Root finders

1/28/96 Harry George Initial version

FROM Arithmetic IMPORT Error;

TYPE
  RealPolynomial2 = ARRAY [0 .. 2] OF R.T;
  ComplexPolynomial2 = ARRAY [0 .. 2] OF C.T;
  RootArray2 = ARRAY [0 .. 1] OF C.T;

PROCEDURE RealQuadratic (READONLY x: RealPolynomial2;  (* coefs*)
  ): RootArray2;

Given a*x^2+b*x+c=0, solve for x.

PROCEDURE ComplexQuadratic (READONLY x: ComplexPolynomial2;  (* coefs*)
  ): RootArray2 RAISES {Error};

Given a*x^2+b*x+c=0, solve for x.

PROCEDURE RealNewtonMaehli (x: RRt.T; ): REF CRt.RootArray RAISES {Error};
PROCEDURE ComplexNewtonMaehli (x: CRt.T; ): REF CRt.RootArray
  RAISES {Error};

determine all roots of a polynomial with a variant of the Newton-Maehli method, that is the roots are not removed successively but the roots are refined simultanous

END RootApproximation.