15 5 decimal

&[DATE] &[TIME] - &[FILENAME]

PRIMER 4.7

x1 0.7 0.05 0.15 0.05 0.05 0071 mat :

x1 0.7 0.05 0.15 0.05 0.05 0071 mat

S 1 - 0001 - 0001 - 2 - 5 - 0.5 - 00013124073 mat :

N.s 3 :

N.c 7 :

N.r 3 :

d N.s 1 - (N.c * N.r + :

d 17

F1 1 :

x2 000 0.21 0.79 0071 mat :

Δν S 1 col sum S 2 col sum S 3 col sum 31 mat :

Δν 01 0.5 - 31 mat

x3 100000071 mat :

ε 0.1 0.1 0.1 31 mat :

F2 10 :

F3 10 :

X0 ε F2 F3 x3 stack :

ε F2 F3 x3 f kbilans F1 x1 * F2 x2 * + F3 x3 * - S ε * + : ubilans 1 el F1 F2 + F3 - Δν ε * + : sastav3 x3 1 el x3 2 el x3 3 el x3 7 el 0.1235 - stack : f kbilans ubilans sastav3 stack : 41 line :

ε F2 F3 x3 f 9.4 - 0.05 - 0.05 1.4 7.95 0.5 0.6 1.05 100 0.1235 - 121 mat

XX# X# 1 el X# 2 el X# 3 el X# 4 el X# 5 el X# 6 el X# 7 el X# 8 el X# 9 el X# 10 el X# 11 el X# 12 el 121 mat :

ε F2 F3 x3 f 9.4 - 0.05 - 0.05 1.4 7.95 0.5 0.6 1.05 100 0.1235 - 121 mat

SOLUTION - do not know how to get it???

X f1 i 13 range ε i el X i el : for F2 X 4 el : F3 X 5 el : i 612 range x3 i 5 - el X i el : for kbilans F1 x1 * F2 x2 * + F3 x3 * - S ε * + : ubilans F1 F2 + F3 - Δν ε * + : sastav3 x3 1 el x3 2 el x3 3 el x3 7 el 0.1235 - stack : f kbilans ubilans sastav3 stack : 81 line :

X0 f1 9.4 - 0.05 - 0.05 1.4 7.95 0.5 0.6 1.05 100 0.1235 - 121 mat

res 0.7 0.05 0.15 11.98047 12.95547 000 0.06491 0.73441 0.07719 0.1235 121 mat :

res f1 000 4.08577 105 -^* - 5.54227 105 -^* - 3.27293 105 -^* - 5.45 107 -^* - 00000121 mat

X f2 i 13 range ε i el X i el : for F2 X 4 el : F3 X 5 el : i 612 range x3 i 5 - el X i el : for kbilans F1 x1 * F2 x2 * + F3 x3 * - S ε * + : ubilans F1 F2 + F3 - Δν ε * + : sastav3 x3 1 el x3 2 el x3 3 el x3 7 el 0.1235 - stack : i 17 range f i el kbilans i el : for f 8 el ubilans : i 912 range f i el sastav3 i 8 - el : for f 111 line :

X0 f2 9.4 - 0.05 - 0.05 1.4 7.95 0.5 0.6 1.05 100 0.1235 - 121 mat

res f2 000 4.08577 105 -^* - 5.54227 105 -^* - 3.27293 105 -^* - 5.45 107 -^* - 00000121 mat

X f3 i 13 range ε i el X i el : for F2 X 4 el : F3 X 5 el : i 612 range x3 i 5 - el X i el : for i 17 range kbilans i el F1 x1 * (i el F2 x2 * (i el + F3 x3 * (i el - S ε * (i el + : for ubilans F1 F2 + F3 - Δν ε * + : sastav3 x3 1 el x3 2 el x3 3 el x3 7 el 0.1235 - stack : i 17 range f i el kbilans i el : for f 8 el ubilans : i 912 range f i el sastav3 i 8 - el : for f 111 line :

X0 f3 9.4 - 0.05 - 0.05 1.4 7.95 0.5 0.6 1.05 100 0.1235 - 121 mat

res f3 000 4.08577 105 -^* - 5.54227 105 -^* - 3.27293 105 -^* - 5.45 107 -^* - 00000121 mat

DIFFERENT JACOBIANS - FUNCTIONS ARE NONLINEAR BUT JACOBIAN CONSTANT VALUES!!!!

XX# f3 XX# Jacob 1 - 0000000000001 - 0000000000001 - 0000000002 - 5 - 12 / - 0000000000000000000001310000000002400000000000112 / - 11 - 00000000000010000000000001000000000000100000000000000011212 mat

XX# f1 XX# Jacob #

XX# f2 XX# Jacob 00000000000000000000000000000000000000000000000000000000000000000000000000000000000001 0.5 - 11 - 00000000000010000000000001000000000000100000000000000011212 mat

XX# f1 XX# X0 roots #

XX# f2 XX# X0 roots #

X# f1 res 1010 -^100 Broyden 0.700000000002705 0.0499999999998 0.150000000001129 11.980465587103 12.9554655870955 2.773 1012 -^* - 1.897 1012 -^* 2.388 1012 -^* - 0.064906796875273 0.734405703125275 0.077187500000903 0.1235 121 mat

res 0.7 0.05 0.15 11.98047 12.95547 000 0.06491 0.73441 0.07719 0.1235 121 mat

XX# f3 XX# X0 roots #

X# f2 X0 105 -^100 NewtonRaphson #

X# f3 X0 105 -^100 Broyden #

X# f1 X0 HRE #

X# f1 X0 HRE.B #

X# f3 res HRE.B 0.7 0.05 0.15 11.98047 12.95547 000 0.06491 0.73441 0.07719 0.1235 121 mat

X# f2 res 1010 -^100 Broyden 0.700000000002705 0.0499999999998 0.150000000001129 11.980465587103 12.9554655870955 2.773 1012 -^* - 1.897 1012 -^* 2.388 1012 -^* - 0.064906796875273 0.734405703125275 0.077187500000903 0.1235 121 mat

res 0.7 0.05 0.15 11.98047 12.95547 000 0.06491 0.73441 0.07719 0.1235 121 mat

X# f3 res 1010 -^100 Broyden 0.700000000002586 0.049999999999748 0.15000000000014 11.9804655871016 12.9554655870955 3.696 1012 -^* - 2.145 1012 -^* 1.941 1012 -^* - 0.064906796875276 0.734405703125276 0.077187500000797 0.1235 121 mat

res 0.7 0.05 0.15 11.98047 12.95547 000 0.06491 0.73441 0.07719 0.1235 121 mat

MUCH BETTER INITIAL CONDITIONS!

X01 0.5 0.05 0.15 1213000 0.07 0.73 0.07 0.13 121 mat :

X# f1 X01 HRE.B #

X# f1 X01 1010 -^100 Broyden #

X# f2 X01 HRE.B #

X# f2 X01 1010 -^100 Broyden #

X# f3 X01 HRE.B #

X# f3 X01 1010 -^100 Broyden #

XX f1 XX 4 el 00021100 / 79100 / 0071 mat * XX 5 el XX 6 el XX 7 el XX 8 el XX 9 el XX 10 el XX 11 el XX 12 el 71 mat * - 710 XX 1 el * - 10 / 120 XX 2 el * - 20 / 320 XX 3 el * - 20 / 11022 XX 1 el * 5 XX 2 el * + (* XX 3 el + (* - 20 / 120 / 3 XX 2 el * XX 1 el + XX 3 el + 22 XX 2 el * XX 1 el + (* 71 mat + 21 XX 2 el + XX 4 el + XX 5 el - (* XX 3 el - 2 / XX 6 el XX 7 el XX 8 el 247 - 2000 XX 12 el * + 2000 / 41 mat stack

XX f2 XX 4 el 00021100 / 79100 / 0071 mat * XX 5 el XX 6 el XX 7 el XX 8 el XX 9 el XX 10 el XX 11 el XX 12 el 71 mat * - 710 XX 1 el * - 10 / 120 XX 2 el * - 20 / 320 XX 3 el * - 20 / 11022 XX 1 el * 5 XX 2 el * + (* XX 3 el + (* - 20 / 120 / 3 XX 2 el * XX 1 el + XX 3 el + 22 XX 2 el * XX 1 el + (* 71 mat + 1 el XX 4 el 00021100 / 79100 / 0071 mat * XX 5 el XX 6 el XX 7 el XX 8 el XX 9 el XX 10 el XX 11 el XX 12 el 71 mat * - 710 XX 1 el * - 10 / 120 XX 2 el * - 20 / 320 XX 3 el * - 20 / 11022 XX 1 el * 5 XX 2 el * + (* XX 3 el + (* - 20 / 120 / 3 XX 2 el * XX 1 el + XX 3 el + 22 XX 2 el * XX 1 el + (* 71 mat + 2 el XX 4 el 00021100 / 79100 / 0071 mat * XX 5 el XX 6 el XX 7 el XX 8 el XX 9 el XX 10 el XX 11 el XX 12 el 71 mat * - 710 XX 1 el * - 10 / 120 XX 2 el * - 20 / 320 XX 3 el * - 20 / 11022 XX 1 el * 5 XX 2 el * + (* XX 3 el + (* - 20 / 120 / 3 XX 2 el * XX 1 el + XX 3 el + 22 XX 2 el * XX 1 el + (* 71 mat + 3 el XX 4 el 00021100 / 79100 / 0071 mat * XX 5 el XX 6 el XX 7 el XX 8 el XX 9 el XX 10 el XX 11 el XX 12 el 71 mat * - 710 XX 1 el * - 10 / 120 XX 2 el * - 20 / 320 XX 3 el * - 20 / 11022 XX 1 el * 5 XX 2 el * + (* XX 3 el + (* - 20 / 120 / 3 XX 2 el * XX 1 el + XX 3 el + 22 XX 2 el * XX 1 el + (* 71 mat + 4 el XX 4 el 00021100 / 79100 / 0071 mat * XX 5 el XX 6 el XX 7 el XX 8 el XX 9 el XX 10 el XX 11 el XX 12 el 71 mat * - 710 XX 1 el * - 10 / 120 XX 2 el * - 20 / 320 XX 3 el * - 20 / 11022 XX 1 el * 5 XX 2 el * + (* XX 3 el + (* - 20 / 120 / 3 XX 2 el * XX 1 el + XX 3 el + 22 XX 2 el * XX 1 el + (* 71 mat + 5 el XX 4 el 00021100 / 79100 / 0071 mat * XX 5 el XX 6 el XX 7 el XX 8 el XX 9 el XX 10 el XX 11 el XX 12 el 71 mat * - 710 XX 1 el * - 10 / 120 XX 2 el * - 20 / 320 XX 3 el * - 20 / 11022 XX 1 el * 5 XX 2 el * + (* XX 3 el + (* - 20 / 120 / 3 XX 2 el * XX 1 el + XX 3 el + 22 XX 2 el * XX 1 el + (* 71 mat + 6 el XX 4 el 00021100 / 79100 / 0071 mat * XX 5 el XX 6 el XX 7 el XX 8 el XX 9 el XX 10 el XX 11 el XX 12 el 71 mat * - 710 XX 1 el * - 10 / 120 XX 2 el * - 20 / 320 XX 3 el * - 20 / 11022 XX 1 el * 5 XX 2 el * + (* XX 3 el + (* - 20 / 120 / 3 XX 2 el * XX 1 el + XX 3 el + 22 XX 2 el * XX 1 el + (* 71 mat + 7 el 21 XX 2 el + XX 4 el + XX 5 el - (* XX 3 el - 2 / XX 6 el XX 7 el XX 8 el 247 - 2000 XX 12 el * + 2000 / 121 mat

XX f3 710 XX 4 el 00021100 / 79100 / 0071 mat * 1 el - XX 5 el XX 6 el XX 7 el XX 8 el XX 9 el XX 10 el XX 11 el XX 12 el 71 mat * 1 el + XX 1 el + (* - 10 / 120 XX 4 el 00021100 / 79100 / 0071 mat * 2 el - XX 5 el XX 6 el XX 7 el XX 8 el XX 9 el XX 10 el XX 11 el XX 12 el 71 mat * 2 el + XX 2 el + (* - 20 / 320 XX 4 el 00021100 / 79100 / 0071 mat * 3 el - XX 5 el XX 6 el XX 7 el XX 8 el XX 9 el XX 10 el XX 11 el XX 12 el 71 mat * 3 el + XX 3 el + (* - 20 / 11022 XX 1 el * 5 XX 2 el * + XX 4 el 00021100 / 79100 / 0071 mat * 4 el - XX 5 el XX 6 el XX 7 el XX 8 el XX 9 el XX 10 el XX 11 el XX 12 el 71 mat * 4 el + (* XX 3 el + (* - 20 / 120 XX 4 el 00021100 / 79100 / 0071 mat * 5 el XX 5 el XX 6 el XX 7 el XX 8 el XX 9 el XX 10 el XX 11 el XX 12 el 71 mat * 5 el - (* + 20 / 3 XX 2 el * XX 4 el 00021100 / 79100 / 0071 mat * 6 el + XX 5 el XX 6 el XX 7 el XX 8 el XX 9 el XX 10 el XX 11 el XX 12 el 71 mat * 6 el - XX 1 el + XX 3 el + 22 XX 2 el * XX 1 el + (* XX 4 el 00021100 / 79100 / 0071 mat * 7 el + XX 5 el XX 6 el XX 7 el XX 8 el XX 9 el XX 10 el XX 11 el XX 12 el 71 mat * 7 el - 21 XX 2 el + XX 4 el + XX 5 el - (* XX 3 el - 2 / XX 6 el XX 7 el XX 8 el 247 - 2000 XX 12 el * + 2000 / 121 mat