[Axiom-developer] [numerical linear algebra]

[anonymous]

++added:

Why would you expect to be able to recover the characteristic polynomial?  
There is always round off error for finite precision arithmetic. For integer 
approximations, it would be worse.
The command {\tt solve: (Polynomial Fraction Integer, PositiveInteger)->List 
Equation Polynomial Integer}
solves the equation over the integers, so it is {\it not} accurate. For example:

\begin{axiom}
solve(x+11/10,3)
ev:= solve(rhs(eigen.1),1.0*10^(-50))
cp:= reduce(*, [rhs(x)-lhs(x) for x in ev])
\end{axiom}

--
forwarded from http://page.axiom-developer.org/zope/mathaction/address@hidden