(35S) calculate Eigenvalues of Symmetric Matrices using Jacobi’s Method

[dmh]

This 35S program is a port of the HP 41C program written by Valentin Albillo in 1980.

Valentin's paper and original 41C program: https://albillo.hpcalc.org/programs/HP%2...Jacobi.pdf

Demonstration of 35S program: https://youtu.be/vuQlEPFerCE

Code:

E001    LBL E
E002    CF 10
E003    INPUT N
E004    10
E005    STO I
E006    RCL N
E007    1 E3
E008    /
E009    1
E010    +
E011    STO P
E012    RCL P
E013    INTG
E014    STO R
E015    RCL P
E016    STO Q
E017    INTG
E018    STO C
E019    EQN [R,C]
E020    PSE
E021    1
E022    STO +C
E023    RCL Q
E024    INPUT X
E025    STO (I)
E026    ISG I
E027    X<>Y
E028    ISG Q
E029    GTO E019
E030    ISG P
E031    GTO E012
E032    CF 0
E033    2
E034    STO Q
E035    1
E036    STO P
E037    CF 0
E038    RCL Q
E039    XEQ E195
E040    X<>Y
E041    RND
E042    X=0?
E043    GTO E164
E044    SF 0
E045    LAST X
E046    REGZ
E047    STO J
E048    RDN
E049    ST - (J)
E050    STO O
E051    ST+ O
E052    CHS
E053    STO T
E054    RCL P
E055    RCL P
E056    XEQ E195
E057    STO V
E058    RDN
E059    STO Z
E060    RCL Q
E061    RCL Q
E062    XEQ E195
E063    STO W
E064    RDN
E065    RCL Z
E066    X<>Y
E067    -
E068    STO U
E069    2
E070    /
E071    STO Y
E072    SGN
E073    x=0?
E074    1
E075    RCL T
E076    RCL Y
E077    X^2
E078    RCL T
E079    X^2
E080    +
E081    SQRT
E082    /
E083    *
E084    1
E085    REGY
E086    X^2
E087    -
E088    SQRT
E089    1
E090    +
E091    REGX
E092    +
E093    SQRT
E094    /
E095    STO T
E096    STO * O
E097    X^2
E098    ST * U
E099    1
E100    STO I
E101    X<>Y
E102    -
E103    SQRT
E104    ST * O
E105    X<> U
E106    RCL O
E107    +
E108    RCL V
E109    STO J
E110    RDN
E111    ST - (J)
E112    RCL W
E113    STO J
E114    RDN
E115    ST + (J)
E116    RCL Q
E117    RCL I
E118    X=Y?
E119    GTO E158
E120    RCL P
E121    X=Y?
E122    GTO E158
E123    XEQ E195
E124    STO S
E125    RDN
E126    STO V
E127    RCL U
E128    STO W
E129    *
E130    STO Y
E131    RCL Q
E132    RCL I
E133    XEQ E195
E134    STO Z
E135    RDN
E136    STO * W
E137    RCL T
E138    STO * V
E139    *
E140    RCL Y
E141    X<>Y
E142    -
E143    RCL S
E144    STO J
E145    RDN
E146    STO (J)
E147    RCL V
E148    RCL W
E149    +
E150    RCL Z
E151    STO J
E152    RDN
E153    STO (J)
E154    RCL I
E155    RCL N
E156    X<=Y?
E157    GTO E164
E158    1
E159    STO + I
E160    RCL N
E161    RCL I
E162    X<=y?
E163    GTO E116
E164    RCL Q
E165    RCL P
E166    1
E167    +
E168    X<>Y?
E169    GTO E036
E170    X<>Y
E171    RCL N
E172    X<>Y?
E173    ISG Q
E174    X<>Y
E175    X<>Y?
E176    GTO E035
E177    FS? 0
E178    GTO E032
E179    RCL N
E180    1 E3
E181    /
E182    1
E183    +
E184    STO T
E185    RCL T
E186    INTG
E187    ENTER
E188    XEQ E195
E189    X<>Y
E190    STO X
E191    VIEW X
E192    ISG T
E193    GTO E185
E194    RTN
E195    X>Y?
E196    X<>Y
E197    RCL N
E198    ENTER
E199    +
E200    X<>Y
E201    -
E202    -1
E203    LAST X
E204    +
E205    *
E206    2
E207    /
E208    +
E209    9
E210    +
E211    STO J
E212    RCL (J)
E213    X<>Y
E214    RTN