(35S) calculate Eigenvalues of Symmetric Matrices using Jacobi’s Method
07-14-2023, 11:45 AM
Post: #1
 dmh Senior Member Posts: 353 Joined: Nov 2020
(35S) calculate Eigenvalues of Symmetric Matrices using Jacobi’s Method
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