(15C) Solving a system of nonlinear equations

[Marcel]

Hi!
I upload two documents for the HP 15C for solving a system of three nonlinear equations, with the program take in the Advanced Functions Handbook, pp.122-127.

The Newton's Method is used and explained in one of the document with the example I use to solve with the HP 15C.

With a modern calculator, the procedure is more simple but this is more fun with the 15C !

Marcel.

[Namir]

Nice code and article! Thanks for sharing!!!

Namir

[Marcel]

Hi Namir!

Thank you very much for your comment.

Marcel.

[Thomas Klemm]

Here's the listing of the program that can be used with this HP-15C simulator:

Code:

   001 { 42 21 11 } f LBL A
   002 {        3 } 3
   003 {       36 } ENTER
   004 { 42 23 13 } f DIM C
   005 {        1 } 1
   006 { 42 23 12 } f DIM B
   007 {    32 12 } GSB B
   008 { 45 16 11 } RCL MATRIX A
   009 { 45 16 12 } RCL MATRIX B
   010 { 45 16 13 } RCL MATRIX C
   011 { 42 26 14 } f RESULT D
   012 {       10 } /
   013 { 42 26 11 } f RESULT A
   014 {       30 } -
   015 {    43 36 } g LSTx
   016 { 42 16  8 } f MATRIX 8
   017 { 45 16 12 } RCL MATRIX B
   018 { 42 16  8 } f MATRIX 8
   019 {    43 32 } g RTN
   020 { 42 21 12 } f LBL B
   021 { 42 16  1 } f MATRIX 1
   022 {        2 } 2
   023 {    44 25 } STO I
   024 {  45 11 u } u RCL A
   025 { 42 21  1 } f LBL 1
   026 {    44 24 } STO (i)
   027 {        1 } 1
   028 { 44 40 25 } STO + I
   029 {  45 11 u } u RCL A
   030 {    22  1 } GTO 1
   031 {    44 24 } STO (i)
   032 { 42 16  1 } f MATRIX 1
   033 {        3 } 3
   034 { 45 20  2 } RCL * 2
   035 {    45  3 } RCL 3
   036 {    45  4 } RCL 4
   037 {       20 } *
   038 {       24 } COS
   039 {       30 } -
   040 {       48 } .
   041 {        5 } 5
   042 {       30 } -
   043 {  44 12 u } u STO B
   044 {    45  2 } RCL 2
   045 {    43 11 } g x^2
   046 {    45  3 } RCL 3
   047 {       48 } .
   048 {        1 } 1
   049 {       40 } +
   050 {    43 11 } g x^2
   051 {        8 } 8
   052 {        1 } 1
   053 {       20 } *
   054 {       30 } -
   055 {    45  4 } RCL 4
   056 {       23 } SIN
   057 {       40 } +
   058 {        1 } 1
   059 {       48 } .
   060 {        0 } 0
   061 {        6 } 6
   062 {       40 } +
   063 {  44 12 u } u STO B
   064 {    45  2 } RCL 2
   065 {    45  3 } RCL 3
   066 {       20 } *
   067 {       16 } CHS
   068 {       12 } eX
   069 {        2 } 2
   070 {        0 } 0
   071 { 45 20  4 } RCL * 4
   072 {       40 } +
   073 {    43 26 } g PI
   074 {        1 } 1
   075 {        0 } 0
   076 {       20 } *
   077 {        3 } 3
   078 {       30 } -
   079 {        3 } 3
   080 {       10 } /
   081 {       40 } +
   082 {  44 12 u } u STO B
   083 { 42 21  2 } f LBL 2
   084 { 42 16  1 } f MATRIX 1
   085 {        3 } 3
   086 {  44 13 u } u STO C
   087 {    45  4 } RCL 4
   088 {    45  3 } RCL 3
   089 {    45  4 } RCL 4
   090 {       20 } *
   091 {       23 } SIN
   092 {       20 } *
   093 {  44 13 u } u STO C
   094 {    45  3 } RCL 3
   095 {    45  3 } RCL 3
   096 {    45  4 } RCL 4
   097 {       20 } *
   098 {       23 } SIN
   099 {       20 } *
   100 {  44 13 u } u STO C
   101 {        2 } 2
   102 { 45 20  2 } RCL * 2
   103 {  44 13 u } u STO C
   104 {        1 } 1
   105 {        6 } 6
   106 {        2 } 2
   107 {       16 } CHS
   108 {    45  3 } RCL 3
   109 {       48 } .
   110 {        1 } 1
   111 {       40 } +
   112 {       20 } *
   113 {  44 13 u } u STO C
   114 {    45  4 } RCL 4
   115 {       24 } COS
   116 {  44 13 u } u STO C
   117 {    45  2 } RCL 2
   118 {    45  3 } RCL 3
   119 {       20 } *
   120 {       16 } CHS
   121 {       12 } eX
   122 { 45 20  3 } RCL * 3
   123 {       16 } CHS
   124 {  44 13 u } u STO C
   125 {    45  2 } RCL 2
   126 {    45  3 } RCL 3
   127 {       20 } *
   128 {       16 } CHS
   129 {       12 } eX
   130 { 45 20  2 } RCL * 2
   131 {       16 } CHS
   132 {  44 13 u } u STO C
   133 {        2 } 2
   134 {        0 } 0
   135 {  44 13 u } u STO C
   136 { 42 21  3 } f LBL 3
   137 {    43 32 } g RTN