Re: eignenvalue for nonsymmetric matrix

[Mike Miller]

> >   octave:1> stiff = [1,2,3,4;5,6,7,8;9,10,11,12;13,14,15,16]

By the way, this gives the same result:

reshape([1:16],4,4)'

and it's extensible:

N=4; stiff=reshape([1:N^2],N,N)'

Just change N.  You can then do things like this:

N=30; stiff=reshape([1:N^2],N,N)'; eig(stiff)

Mike



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------