[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Axiom-developer] [numerical linear algebra]

From: anonymous
Subject: [Axiom-developer] [numerical linear algebra]
Date: Fri, 11 Feb 2005 08:05:55 -0600


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:

ev:= solve(rhs(eigen.1),1.0*10^(-50))
cp:= reduce(*, [rhs(x)-lhs(x) for x in ev])

forwarded from

reply via email to

[Prev in Thread] Current Thread [Next in Thread]