21 Oct 16:33 2014

### diagonalize matrix

K M <erik0709 <at> live.dk>

2014-10-21 14:33:41 GMT

2014-10-21 14:33:41 GMT

Hello Everyone!

I started to use wxMaxima today :) and now trying to diagonalize a 3x3 matrix looks like:

matrix( [0,0,a], [0,b,-a], [a,-a,c]).

(a,b,c are all symbolic).

with an example:

assume(a>0,b>0); M : matrix([a,a+b,a+b],[a+b,a,a+b],[a+b,a+b,a]); load("eigen"); [myeigval,myeigvec]:similaritytransform(ev(M,hermitianmatrix));

; however end up with errors

(%o3) "C:/Maxima-5.31.2/share/maxima/5.31.2/share/matrix/eigen.mac"

length: argument cannot be a symbol; found M

#0: eigenvectors(mat=M)(eigen.mac line 107)

#1: uniteigenvectors(m=M)(eigen.mac line 154)

#2: similaritytransform(mat=M)(eigen.mac line 169)

-- an error. To debug this try: debugmode(true);

Could you help me with this?

Thanks!

/Erik

I started to use wxMaxima today :) and now trying to diagonalize a 3x3 matrix looks like:

matrix( [0,0,a], [0,b,-a], [a,-a,c]).

(a,b,c are all symbolic).

with an example:

assume(a>0,b>0); M : matrix([a,a+b,a+b],[a+b,a,a+b],[a+b,a+b,a]); load("eigen"); [myeigval,myeigvec]:similaritytransform(ev(M,hermitianmatrix));

; however end up with errors

(%o3) "C:/Maxima-5.31.2/share/maxima/5.31.2/share/matrix/eigen.mac"

length: argument cannot be a symbol; found M

#0: eigenvectors(mat=M)(eigen.mac line 107)

#1: uniteigenvectors(m=M)(eigen.mac line 154)

#2: similaritytransform(mat=M)(eigen.mac line 169)

-- an error. To debug this try: debugmode(true);

Could you help me with this?

Thanks!

/Erik

------------------------------------------------------------------------------ Comprehensive Server Monitoring with Site24x7. Monitor 10 servers for $9/Month. Get alerted through email, SMS, voice calls or mobile push notifications. Take corrective actions from your mobile device. http://p.sf.net/sfu/Zoho

_______________________________________________ Maxima-discuss mailing list Maxima-discuss <at> lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/maxima-discuss