Tuesday, December 19, 2023

Hilbert Matrix 20x20 LDLT decomposition

 The following is LDLT decomposition of 20x20 Hilbert Matrix A.


It seems Firefox render this page correctly.

Hilbert Matrix H=[ 1 1 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10 1 11 1 12 1 13 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10 1 11 1 12 1 13 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1 21 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10 1 11 1 12 1 13 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1 21 1 22 1 4 1 5 1 6 1 7 1 8 1 9 1 10 1 11 1 12 1 13 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1 21 1 22 1 23 1 5 1 6 1 7 1 8 1 9 1 10 1 11 1 12 1 13 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1 21 1 22 1 23 1 24 1 6 1 7 1 8 1 9 1 10 1 11 1 12 1 13 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1 21 1 22 1 23 1 24 1 25 1 7 1 8 1 9 1 10 1 11 1 12 1 13 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1 21 1 22 1 23 1 24 1 25 1 26 1 8 1 9 1 10 1 11 1 12 1 13 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1 21 1 22 1 23 1 24 1 25 1 26 1 27 1 9 1 10 1 11 1 12 1 13 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1 21 1 22 1 23 1 24 1 25 1 26 1 27 1 28 1 10 1 11 1 12 1 13 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1 21 1 22 1 23 1 24 1 25 1 26 1 27 1 28 1 29 1 11 1 12 1 13 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1 21 1 22 1 23 1 24 1 25 1 26 1 27 1 28 1 29 1 30 1 12 1 13 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1 21 1 22 1 23 1 24 1 25 1 26 1 27 1 28 1 29 1 30 1 31 1 13 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1 21 1 22 1 23 1 24 1 25 1 26 1 27 1 28 1 29 1 30 1 31 1 32 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1 21 1 22 1 23 1 24 1 25 1 26 1 27 1 28 1 29 1 30 1 31 1 32 1 33 1 15 1 16 1 17 1 18 1 19 1 20 1 21 1 22 1 23 1 24 1 25 1 26 1 27 1 28 1 29 1 30 1 31 1 32 1 33 1 34 1 16 1 17 1 18 1 19 1 20 1 21 1 22 1 23 1 24 1 25 1 26 1 27 1 28 1 29 1 30 1 31 1 32 1 33 1 34 1 35 1 17 1 18 1 19 1 20 1 21 1 22 1 23 1 24 1 25 1 26 1 27 1 28 1 29 1 30 1 31 1 32 1 33 1 34 1 35 1 36 1 18 1 19 1 20 1 21 1 22 1 23 1 24 1 25 1 26 1 27 1 28 1 29 1 30 1 31 1 32 1 33 1 34 1 35 1 36 1 37 1 19 1 20 1 21 1 22 1 23 1 24 1 25 1 26 1 27 1 28 1 29 1 30 1 31 1 32 1 33 1 34 1 35 1 36 1 37 1 38 1 20 1 21 1 22 1 23 1 24 1 25 1 26 1 27 1 28 1 29 1 30 1 31 1 32 1 33 1 34 1 35 1 36 1 37 1 38 1 39 ]

L=[ 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 2 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 3 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 4 9 10 3 2 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 5 4 5 12 7 2 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 6 5 7 25 14 25 9 5 2 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 7 9 14 25 14 10 3 45 11 3 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 8 7 12 7 4 245 66 245 44 147 26 7 2 1 0 0 0 0 0 0 0 0 0 0 0 0 1 9 8 15 56 33 392 99 980 143 112 13 112 15 4 1 0 0 0 0 0 0 0 0 0 0 0 1 10 27 55 18 11 588 143 1134 143 756 65 63 5 162 17 9 2 1 0 0 0 0 0 0 0 0 0 0 1 11 5 11 225 143 600 143 1260 143 189 13 315 17 300 17 225 19 5 1 0 0 0 0 0 0 0 0 0 1 12 11 26 275 182 55 13 495 52 7623 442 847 34 9075 323 1815 76 605 42 11 2 1 0 0 0 0 0 0 0 0 1 13 36 91 132 91 55 13 4455 442 4356 221 10164 323 13068 323 5445 133 220 7 396 23 6 1 0 0 0 0 0 0 0 1 14 13 35 39 28 143 34 715 68 14157 646 61347 1615 122694 2261 16731 266 9295 161 1859 46 507 25 13 2 1 0 0 0 0 0 0 1 15 7 20 91 68 637 153 7007 646 77077 3230 143143 3230 22308 323 39039 437 13013 138 91091 1150 1274 25 637 27 7 1 0 0 0 0 0 1 16 45 136 175 136 15925 3876 28665 2584 33033 1292 65065 1292 1254825 14858 418275 3496 13013 92 63063 460 637 6 2275 36 1575 58 15 2 1 0 0 0 0 1 17 16 51 400 323 3920 969 3640 323 8736 323 416416 7429 743600 7429 66924 437 4576 23 224224 1035 2912 15 36400 261 2240 29 960 31 8 1 0 0 0 1 18 17 57 68 57 680 171 2380 209 12376 437 80444 1311 252824 2185 82654 437 165308 621 330616 1035 420784 1305 210392 783 161840 899 2890 31 1156 33 17 2 1 0 0 1 19 27 95 153 133 816 209 55080 4807 12852 437 723996 10925 286416 2185 495924 2185 165308 483 1487772 3335 360672 725 420784 899 327726 899 78030 341 1224 11 1377 35 9 1 0 1 20 19 70 171 154 969 253 2907 253 17442 575 40698 575 16796 115 214149 805 29838094 70035 29838094 50025 16275324 22475 1356277 1798 6572727 9889 165699 341 110466 385 18411 140 3249 74 19 2 1 ]

D=[ 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 180 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 2800 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 44100 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 698544 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 11099088 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 176679360 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 2815827300 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 44914183600 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 716830370256 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 11445589052352 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 182811491808400 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 2920656969720000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 46670906271240000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 745904795339462400 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 11922821963004219300 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 190600129650794094000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 3047248986392325330000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 48722219250572027160000 ]

Hilbert Matrix H= L D L T

No comments:

Post a Comment