Description

Eigenvalue example.

Contributor: Erwin Kalvelagen, October 2008.

Small Model of Type : GAMS

Category : GAMS Test library

Main file : eigval01.gms

$title Test eigenvalue utility (EIGVAL01,SEQ=409)

$ontext

Eigenvalue example.

Contributor: Erwin Kalvelagen, October 2008.


octave:1> a=[9 1 1; 1 9 1; 1 1 9]
a =

  9  1  1
  1  9  1
  1  1  9

octave:2> eig(a)
ans =

   8
   8
  11


$offtext

set i /i1*i3/;
alias (i,j);

table a(i,j)
      i1   i2   i3
i1     9    1    1
i2     1    9    1
i3     1    1    9

;

parameter e(i) 'eigenvalues';

$libInclude linalg eigenvalue i a e
abort$execerror 'Trouble with calculating Eigenvalue';

display a,e;