File: eigen.pcm

package info (click to toggle)
extrema 4.3.6-1
  • links: PTS
  • area: main
  • in suites: lenny
  • size: 19,212 kB
  • ctags: 6,452
  • sloc: cpp: 86,428; sh: 8,229; makefile: 814
file content (17 lines) | stat: -rw-r--r-- 239 bytes parent folder | download | duplicates (3) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
m=[[2;-1;0;0];[-1;2;-1;0];[0;-1;2;-1];[0;0;-1;2]]
e=eigen(m)
='matrix m'
=m
=' '
='matrix e'
=e
!
! check
!
n=vlen(m)[1]
do j = [1:n]
! these should be all zero (or close to zero)
 =' '
 ='j = '//rchar(j)
 =m<>e[*,j]-e[j,n+1]*e[*,j]
enddo