Hi all,
Why can't I get the eigenvalues by the following lines?
Code:
LIB "linalg.lib";
ring r=0,x,lp;
> matrix m[2][2]=1,2,3,4;
> print(m);
1,2,
3,4
> eigenvals(m);
[1]:
_[1]=x2-5x-2
[2]:
1
It seems that the output is the characteristic polynomial instead of the eigenvalues. Is it because it can't be solved in the field of rational numbers?! Let's try the complex numbers as the field of coefficients.
Code:
ring r=complex,x,lp;
// ** redefining r **
> matrix m[2][2]=1,2,3,4;
> print(m);
1,2,
3,4
> eigenvals(m);
? not implemented
? error occurred in or before STDIN line 9: ` return(system("eigenvals",jet(M,0)));`
Login
Register
Forum FAQ
Forum Search
