You found a bug, I'm sorry to say. I'm pretty sure that the bug is due
to a faulty computation of the multiplicities of the roots of the
characteristic polynomial:
(%i26) solve(g33664^2*(g33664^2+1/16)+(g33664^2+1/16)/16,g33664);
(%o26) [g33664=-%i/4,g33664=%i/4]
(%i27) multiplicities;
(%o27) [1,1]
I think (%o27) should be [2,2]. I can't think of a workaround; I'll
file a bug report.
Barton
-----maxima-bounces at math.utexas.edu wrote: -----
>I'm?trying?to?find?the?eigenvalues?of?a?simple?matrix
>
>b:matrix([0,-1/4,0,0],[1/4,0,0,0],[0,0,0,-1/4],[0,0,1/4,0]);
>
>by
>
>eigenvalues(b);
>
>But?Maxima?give?me
>
>part?fell?off?end.
>#0:
>eigenvalues(mat=matrix([0,-1/4,0,0],[1/4,0,0,0],[0,0,0,-1/4],[0,0,1/4,0]))
>(eigen.mac?line?94)
>--?an?error.??To?debug?this?try?debugmode(true);
>
>What?wrong?with?this??Can?anyone?help?
>
>Noparit?Jinuntuya
>Department?of?Physics,?Kasetsart?University
>Thailand.
>
>
>
>
>
>_______________________________________________
>Maxima?mailing?list
>Maxima at math.utexas.edu
>http://www.math.utexas.edu/mailman/listinfo/maxima