-----maxima-bounces at math.utexas.edu wrote: -----
> I noticed that none of the included eigenvalue routines (that I am
> aware of at least) work in a reasonable way with bfloat and irrational
> numbers.
For real selfadjoint matrices, we have double and bigfloat versions of
the Jacobi method.
(%i1) B: matrix([sqrt(2),1],[1,1]);
(%o1) matrix([sqrt(2),1],[1,1])
(%i2) eigens_by_jacobi(B, 'bigfloatfield);
....(deleted)
Barton