+- HP Forums (https://www.hpmuseum.org/forum)
+-- Forum: HP Software Libraries (/forum-10.html)
+--- Forum: General Software Library (/forum-13.html)
+--- Thread: (35S) calculate Eigenvalues of Symmetric Matrices using Jacobi’s Method (/thread-20173.html)
(35S) calculate Eigenvalues of Symmetric Matrices using Jacobi’s Method - dmh - 07-14-2023 11:45 AM
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%20Program%20VA420%20-%20HP-41C%20Eigenvalues%20of%20Symmetric%20Matrices%20-%20Jacobi.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