1 | test_eigen = ->
|
2 | run_test [
|
3 | "eigen(A)",
|
4 | "Stop: eigen: argument is not a square matrix",
|
5 |
|
6 | "eigenval(A)",
|
7 | "eigenval(A)",
|
8 |
|
9 | "eigenvec(A)",
|
10 | "eigenvec(A)",
|
11 |
|
12 | "eigen([1,2])",
|
13 | "Stop: eigen: argument is not a square matrix",
|
14 |
|
15 | "eigen([[1,2],[1,2]])",
|
16 | "Stop: eigen: matrix is not symmetrical",
|
17 |
|
18 | "eigenval([[1,1,1,1],[1,2,3,4],[1,3,6,10],[1,4,10,20]])",
|
19 | "[[0.038016,0.0,0.0,0.0],[0.0,0.453835,0.0,0.0],[0.0,0.0,2.20345,0.0],[0.0,0.0,0.0,26.3047]]",
|
20 |
|
21 | "eigenvec([[1,1,1,1],[1,2,3,4],[1,3,6,10],[1,4,10,20]])",
|
22 | "[[0.308686,-0.72309,0.594551,-0.168412],[0.787275,-0.163234,-0.532107,0.265358],[0.530366,0.640332,0.391832,-0.393897],[0.0601867,0.201173,0.458082,0.863752]]",
|
23 |
|
24 | "eigen(hilbert(50))",
|
25 | "",
|
26 |
|
27 | "1+trace(hilbert(50))-trace(dot(transpose(Q),D,Q))",
|
28 | "1",
|
29 |
|
30 | "D=quote(D)",
|
31 | "",
|
32 |
|
33 | "Q=quote(Q)",
|
34 | "",
|
35 |
|
36 | "A=hilbert(3)",
|
37 | "",
|
38 |
|
39 | "eigen(A)",
|
40 | "",
|
41 |
|
42 | "D-eigenval(A)",
|
43 | "[[0,0,0],[0,0,0],[0,0,0]]",
|
44 |
|
45 | "Q-eigenvec(A)",
|
46 | "[[0,0,0],[0,0,0],[0,0,0]]",
|
47 |
|
48 | "A=quote(A)",
|
49 | "",
|
50 |
|
51 | "D=quote(D)",
|
52 | "",
|
53 |
|
54 | "Q=quote(Q)",
|
55 | "",
|
56 | ]
|