-----maxima-bounces at math.utexas.edu wrote: -----
>is there a Maxima function for the numerical
>computation of eigenystems?
The linearalgebra package has code for the Jacobi
method, but this method is restricted to real
symmetric matrices. I don't know of any other
numerical code for eigenvalues / vectors. Would you
like to contribute?
(%i1) load(linearalgebra)$
(%i2) m : genmatrix(lambda([i,j],1/(i+j-1)),5,5)$
(%i3) eigens_by_jacobi(m,floatfield);
(%o3) [[1.5670506910982309,0.2085342186110134,3.2879287721665683
*10^-6,0.011407491623419795,
3.0589804015119409*10^-4],
And
(%i4) eigens_by_jacobi(m,bigfloatfield);
Barton