5547:
4365:
5542:{\displaystyle {\begin{aligned}\mathbf {x} ^{(1)}&={\begin{bmatrix}0.000&-0.1875\\0.000&-0.1193\end{bmatrix}}{\begin{bmatrix}1.0\\1.0\end{bmatrix}}+{\begin{bmatrix}0.6875\\-0.7443\end{bmatrix}}={\begin{bmatrix}0.5000\\-0.8636\end{bmatrix}}.\\\mathbf {x} ^{(2)}&={\begin{bmatrix}0.000&-0.1875\\0.000&-0.1193\end{bmatrix}}{\begin{bmatrix}0.5000\\-0.8636\end{bmatrix}}+{\begin{bmatrix}0.6875\\-0.7443\end{bmatrix}}={\begin{bmatrix}0.8494\\-0.6413\end{bmatrix}}.\\\mathbf {x} ^{(3)}&={\begin{bmatrix}0.000&-0.1875\\0.000&-0.1193\end{bmatrix}}{\begin{bmatrix}0.8494\\-0.6413\\\end{bmatrix}}+{\begin{bmatrix}0.6875\\-0.7443\end{bmatrix}}={\begin{bmatrix}0.8077\\-0.6678\end{bmatrix}}.\\\mathbf {x} ^{(4)}&={\begin{bmatrix}0.000&-0.1875\\0.000&-0.1193\end{bmatrix}}{\begin{bmatrix}0.8077\\-0.6678\end{bmatrix}}+{\begin{bmatrix}0.6875\\-0.7443\end{bmatrix}}={\begin{bmatrix}0.8127\\-0.6646\end{bmatrix}}.\\\mathbf {x} ^{(5)}&={\begin{bmatrix}0.000&-0.1875\\0.000&-0.1193\end{bmatrix}}{\begin{bmatrix}0.8127\\-0.6646\end{bmatrix}}+{\begin{bmatrix}0.6875\\-0.7443\end{bmatrix}}={\begin{bmatrix}0.8121\\-0.6650\end{bmatrix}}.\\\mathbf {x} ^{(6)}&={\begin{bmatrix}0.000&-0.1875\\0.000&-0.1193\end{bmatrix}}{\begin{bmatrix}0.8121\\-0.6650\end{bmatrix}}+{\begin{bmatrix}0.6875\\-0.7443\end{bmatrix}}={\begin{bmatrix}0.8122\\-0.6650\end{bmatrix}}.\\\mathbf {x} ^{(7)}&={\begin{bmatrix}0.000&-0.1875\\0.000&-0.1193\end{bmatrix}}{\begin{bmatrix}0.8122\\-0.6650\end{bmatrix}}+{\begin{bmatrix}0.6875\\-0.7443\end{bmatrix}}={\begin{bmatrix}0.8122\\-0.6650\end{bmatrix}}.\end{aligned}}}
7204:
1494:
6834:
4184:
6657:
1158:
472:
7199:{\displaystyle {\begin{aligned}\mathbf {x} ^{(1)}&={\begin{bmatrix}0&-1.500\\0&1.071\\\end{bmatrix}}{\begin{bmatrix}1.1\\2.3\\\end{bmatrix}}+{\begin{bmatrix}5.500\\-2.071\\\end{bmatrix}}={\begin{bmatrix}2.050\\0.393\\\end{bmatrix}}.\\\mathbf {x} ^{(2)}&={\begin{bmatrix}0&-1.500\\0&1.071\\\end{bmatrix}}{\begin{bmatrix}2.050\\0.393\\\end{bmatrix}}+{\begin{bmatrix}5.500\\-2.071\\\end{bmatrix}}={\begin{bmatrix}4.911\\-1.651\end{bmatrix}}.\\\mathbf {x} ^{(3)}&=\cdots .\end{aligned}}}
3910:
6386:
125:
1489:{\displaystyle \mathbf {A} =\underbrace {\begin{bmatrix}a_{11}&0&\cdots &0\\a_{21}&a_{22}&\cdots &0\\\vdots &\vdots &\ddots &\vdots \\a_{n1}&a_{n2}&\cdots &a_{nn}\end{bmatrix}} _{\textstyle \mathbf {L} }+\underbrace {\begin{bmatrix}0&a_{12}&\cdots &a_{1n}\\0&0&\cdots &a_{2n}\\\vdots &\vdots &\ddots &\vdots \\0&0&\cdots &0\end{bmatrix}} _{\textstyle \mathbf {U} }.}
1669:
4179:{\displaystyle {\begin{aligned}\mathbf {T} &=-{\begin{bmatrix}0.0625&0.0000\\0.0398&-0.0909\end{bmatrix}}{\begin{bmatrix}0&3\\0&0\end{bmatrix}}={\begin{bmatrix}0.000&-0.1875\\0.000&-0.1194\end{bmatrix}},\\\mathbf {c} &={\begin{bmatrix}0.0625&0.0000\\0.0398&-0.0909\end{bmatrix}}{\begin{bmatrix}11\\13\end{bmatrix}}={\begin{bmatrix}0.6875\\-0.7439\end{bmatrix}}.\end{aligned}}}
6652:{\displaystyle {\begin{aligned}\mathbf {T} &=-{\begin{bmatrix}0.500&0.000\\-0.357&0.143\\\end{bmatrix}}{\begin{bmatrix}0&3\\0&0\\\end{bmatrix}}={\begin{bmatrix}0.000&-1.500\\0.000&1.071\\\end{bmatrix}},\\\mathbf {c} &={\begin{bmatrix}0.500&0.000\\-0.357&0.143\\\end{bmatrix}}{\begin{bmatrix}11\\13\\\end{bmatrix}}={\begin{bmatrix}5.500\\-2.071\\\end{bmatrix}}.\end{aligned}}}
467:{\displaystyle \mathbf {A} ={\begin{bmatrix}a_{11}&a_{12}&\cdots &a_{1n}\\a_{21}&a_{22}&\cdots &a_{2n}\\\vdots &\vdots &\ddots &\vdots \\a_{n1}&a_{n2}&\cdots &a_{nn}\end{bmatrix}},\qquad \mathbf {x} ={\begin{bmatrix}x_{1}\\x_{2}\\\vdots \\x_{n}\end{bmatrix}},\qquad \mathbf {b} ={\begin{bmatrix}b_{1}\\b_{2}\\\vdots \\b_{n}\end{bmatrix}}.}
7865:
8703:
1509:
8327:
9866:
7570:
8338:
3750:
6229:
2172:
3903:
6379:
3303:
5783:
7952:
1664:{\displaystyle {\begin{alignedat}{1}\mathbf {A} \mathbf {x} &=\mathbf {b} \\(\mathbf {L} +\mathbf {U} )\mathbf {x} &=\mathbf {b} \\\mathbf {L} \mathbf {x} +\mathbf {U} \mathbf {x} &=\mathbf {b} \\\mathbf {L} \mathbf {x} &=\mathbf {b} -\mathbf {U} \mathbf {x} \end{alignedat}}}
9602:
6045:
3565:
5626:
1823:
3637:
6119:
1922:
3779:
7860:{\displaystyle {\begin{array}{rrrrl}10x_{1}&-x_{2}&+2x_{3}&&=6,\\-x_{1}&+11x_{2}&-x_{3}&+3x_{4}&=25,\\2x_{1}&-x_{2}&+10x_{3}&-x_{4}&=-11,\\&3x_{2}&-x_{3}&+8x_{4}&=15.\end{array}}}
5881:
3401:
6258:
3200:
7305:
5683:
4370:
972:
2835:
8698:{\displaystyle {\begin{aligned}x_{1}&=3/5=0.6,\\x_{2}&=(3/5)/11+25/11=3/55+25/11=2.3272,\\x_{3}&=-(3/5)/5+(2.3272)/10-11/10=-3/25+0.23272-1.1=-0.9873,\\x_{4}&=-3(2.3272)/8+(-0.9873)/8+15/8=0.8789.\end{aligned}}}
5956:
3476:
4357:
6827:
6839:
8322:{\displaystyle {\begin{aligned}x_{1}&=x_{2}/10-x_{3}/5+3/5,\\x_{2}&=x_{1}/11+x_{3}/11-3x_{4}/11+25/11,\\x_{3}&=-x_{1}/5+x_{2}/10+x_{4}/10-11/10,\\x_{4}&=-3x_{2}/8+x_{3}/8+15/8.\end{aligned}}}
6391:
3915:
9861:{\displaystyle x_{i}^{(k+1)}={\frac {1}{a_{ii}}}\left(b_{i}-\sum _{j<i}a_{ij}x_{j}^{(k+1)}-\sum _{j>i}a_{ij}x_{j}^{(k)}\right),\quad {\begin{array}{l}i=1,2,\ldots ,n\\k=0,1,2,\ldots \end{array}}}
2880:
2427:-th iteration. This means that, unlike the Jacobi method, only one storage vector is required as elements can be overwritten as they are computed, which can be advantageous for very large problems.
5554:
2679:
2627:
1725:
1084:
5964:
5676:
3484:
3195:
8343:
7957:
114:
2720:
5790:
3310:
2230:
8707:
Using the approximations obtained, the iterative procedure is repeated until the desired accuracy has been reached. The following are the approximated solutions after four iterations.
665:
7233:
2358:
2317:
1893:
894:
833:
7535:
7435:
6763:
4290:
2915:
2755:
2393:
750:
597:
7920:
7362:
2973:
Since elements can be overwritten as they are computed in this algorithm, only one storage vector is needed, and vector indexing is omitted. The algorithm goes as follows:
7228:
6725:
6703:
6681:
6253:
6114:
6092:
6070:
4252:
4230:
4208:
3774:
3612:
3590:
2937:
2589:
2556:
2528:
2503:
1852:
1719:
1695:
1046:
1020:
994:
855:
792:
619:
562:
540:
518:
496:
715:
2473:
10167:"Über ein Verfahren, die Gleichungen, auf welche die Methode der kleinsten Quadrate führt, sowie lineäre Gleichungen überhaupt, durch successive Annäherung aufzulösen"
5886:
3406:
10152:
2963:
10394:
8821:
8793:
8765:
8737:
7947:
7562:
7489:
7462:
7389:
3745:{\displaystyle \mathbf {L} ={\begin{bmatrix}16&0\\7&-11\\\end{bmatrix}}\quad {\text{and}}\quad \mathbf {U} ={\begin{bmatrix}0&3\\0&0\end{bmatrix}}.}
2425:
2257:
1131:
6224:{\displaystyle \mathbf {L} ={\begin{bmatrix}2&0\\5&7\\\end{bmatrix}}\quad {\text{and}}\quad \mathbf {U} ={\begin{bmatrix}0&3\\0&0\\\end{bmatrix}}.}
2760:
881:
2167:{\displaystyle x_{i}^{(k+1)}={\frac {1}{a_{ii}}}\left(b_{i}-\sum _{j=1}^{i-1}a_{ij}x_{j}^{(k+1)}-\sum _{j=i+1}^{n}a_{ij}x_{j}^{(k)}\right),\quad i=1,2,\dots ,n.}
7333:
4295:
3632:
1913:
1151:
1104:
770:
10169:[On a process for solving by successive approximation the equations to which the method of least squares leads as well as linear equations generally].
6768:
3898:{\displaystyle \mathbf {L} ^{-1}={\begin{bmatrix}16&0\\7&-11\end{bmatrix}}^{-1}={\begin{bmatrix}0.0625&0.0000\\0.0398&-0.0909\\\end{bmatrix}}.}
6374:{\displaystyle \mathbf {L} ^{-1}={\begin{bmatrix}2&0\\5&7\\\end{bmatrix}}^{-1}={\begin{bmatrix}0.500&0.000\\-0.357&0.143\\\end{bmatrix}}.}
3298:{\displaystyle \mathbf {A} ={\begin{bmatrix}16&3\\7&-11\\\end{bmatrix}}\quad {\text{and}}\quad \mathbf {b} ={\begin{bmatrix}11\\13\end{bmatrix}}.}
5778:{\displaystyle \mathbf {A} ={\begin{bmatrix}2&3\\5&7\\\end{bmatrix}}\quad {\text{and}}\quad \mathbf {b} ={\begin{bmatrix}11\\13\\\end{bmatrix}}.}
10441:
10409:
7564:. Then, repeat iterations until convergence is achieved, or break if the divergence in the solutions start to diverge beyond a predefined level.
10562:
57:. Though it can be applied to any matrix with non-zero elements on the diagonals, convergence is only guaranteed if the matrix is either
10587:
2259:. The procedure is generally continued until the changes made by an iteration are below some tolerance, such as a sufficiently small
10613:
2840:
10582:
10434:
10079:
6040:{\displaystyle \mathbf {T} =-\mathbf {L} ^{-1}\mathbf {U} \quad {\text{and}}\quad \mathbf {c} =\mathbf {L} ^{-1}\mathbf {b} .}
3560:{\displaystyle \mathbf {T} =-\mathbf {L} ^{-1}\mathbf {U} \quad {\text{and}}\quad \mathbf {c} =\mathbf {L} ^{-1}\mathbf {b} .}
10314:
10133:
2632:
2594:
1051:
5646:
3165:
84:
10608:
10541:
10427:
10572:
2688:
10499:
2435:
2179:
1023:
5621:{\displaystyle \mathbf {x} =\mathbf {A} ^{-1}\mathbf {b} \approx {\begin{bmatrix}0.8122\\-0.6650\end{bmatrix}}}
2430:
However, unlike the Jacobi method, the computations for each element are generally much harder to implement in
1818:{\displaystyle \mathbf {x} ^{(k+1)}=\mathbf {L} ^{-1}\left(\mathbf {b} -\mathbf {U} \mathbf {x} ^{(k)}\right).}
624:
10376:
2260:
10531:
2322:
2281:
1857:
797:
10371:
10224:
Bagnara, Roberto (March 1995). "A Unified Proof for the
Convergence of Jacobi and Gauss-Seidel Methods".
10171:
Abhandlungen der
Mathematisch-Physikalischen Klasse der Königlich Bayerischen Akademie der Wissenschaften
7494:
7394:
6733:
4260:
2885:
2725:
2363:
720:
567:
10485:
10414:
7872:
2559:
2446:
1514:
58:
40:
7338:
10536:
10074:
54:
10166:
10450:
10366:
10238:
5876:{\displaystyle \mathbf {x} ^{(k+1)}=\mathbf {L} ^{-1}(\mathbf {b} -\mathbf {U} \mathbf {x} ^{(k)})}
3396:{\displaystyle \mathbf {x} ^{(k+1)}=\mathbf {L} ^{-1}(\mathbf {b} -\mathbf {U} \mathbf {x} ^{(k)})}
2531:
70:
66:
20:
10084:
7211:
6708:
6686:
6664:
6236:
6097:
6075:
6053:
4235:
4213:
4191:
3757:
3595:
3573:
2920:
2572:
2539:
2511:
2486:
2442:. Furthermore, the values at each iteration are dependent on the order of the original equations.
2176:
Notice that the formula uses two summations per iteration which can be expressed as one summation
1835:
1702:
1678:
1029:
1003:
977:
838:
775:
602:
545:
523:
501:
479:
7575:
670:
10399:
8896:
The following iterative procedure produces the solution vector of a linear system of equations:
7300:{\displaystyle \mathbf {x} =\mathbf {A} ^{-1}\mathbf {b} ={\begin{bmatrix}-38\\29\end{bmatrix}}}
10388:
10233:
2452:
7230:
is neither diagonally dominant nor positive definite. Then, convergence to the exact solution
10495:
967:{\displaystyle \mathbf {L} \mathbf {x} ^{(k+1)}=\mathbf {b} -\mathbf {U} \mathbf {x} ^{(k)},}
2942:
2830:{\displaystyle \mathbf {M} \mathbf {x} ^{(k+1)}=\mathbf {N} \mathbf {x} ^{(k)}+\mathbf {b} }
10490:
10285:
10104:
8799:
8771:
8743:
8715:
7925:
7540:
7467:
7440:
7367:
2398:
2271:
The element-wise formula for the Gauss–Seidel method is related to that of the (iterative)
2235:
1916:
1109:
50:
4359:
The closer the guess to the final solution, the fewer iterations the algorithm will need.
8:
10469:
10302:
860:
10326:
10251:
10092:(a "row-oriented" method, whereas Gauss-Seidel is "column-oriented." See, for example,
7318:
3617:
2431:
1898:
1136:
1089:
755:
10342:
10310:
10129:
997:
10505:
10243:
10099:
7208:
In a test for convergence we find that the algorithm diverges. In fact, the matrix
62:
36:
2483:
The convergence properties of the Gauss–Seidel method are dependent on the matrix
10546:
10292:(in German), vol. 9, Göttingen: Köninglichen Gesellschaft der Wissenschaften
10089:
2682:
2566:
The Gauss–Seidel method may converge even if these conditions are not satisfied.
5951:{\displaystyle \mathbf {x} ^{(k+1)}=\mathbf {T} \mathbf {x} ^{(k)}+\mathbf {c} }
3471:{\displaystyle \mathbf {x} ^{(k+1)}=\mathbf {T} \mathbf {x} ^{(k)}+\mathbf {c} }
10464:
10383:
8335:
is the initial approximation, then the first approximate solution is given by:
9594:
10602:
10510:
10298:
2439:
2272:
1674:
The Gauss–Seidel method now solves the left hand side of this expression for
47:
10419:
10346:
542:
is unknown, the Gauss–Seidel method can be used to iteratively approximate
4352:{\displaystyle \mathbf {x} ^{(0)}={\begin{bmatrix}1.0\\1.0\end{bmatrix}}.}
9793:
6822:{\displaystyle \mathbf {x} ^{(0)}={\begin{bmatrix}1.1\\2.3\end{bmatrix}}}
7364:. At any step in a Gauss-Seidel iteration, solve the first equation for
10255:
69:. It was only mentioned in a private letter from Gauss to his student
10526:
10331:
10247:
10093:
10404:
44:
73:
in 1823. A publication was not delivered before 1874 by Seidel.
10577:
10567:
2569:
Golub and Van Loan give a theorem for an algorithm that splits
1722:
on the right hand side. Analytically, this may be written as
10323:
5635:
is strictly diagonally dominant, but not positive definite.
10350:
9595:
Program to solve arbitrary number of equations using Matlab
16:
Iterative method used to solve a linear system of equations
2875:{\displaystyle \mathbf {x} =\mathbf {A} ^{-1}\mathbf {b} }
2360:
that have already been computed, and only the elements of
5638:
2505:. Namely, the procedure is known to converge if either:
1832:
However, by taking advantage of the triangular form of
7310:
7273:
7130:
7098:
7069:
7030:
6972:
6940:
6911:
6872:
6798:
6618:
6589:
6550:
6492:
6453:
6414:
6334:
6286:
6187:
6136:
5751:
5700:
5594:
5551:
As expected, the algorithm converges to the solution:
5508:
5476:
5444:
5402:
5341:
5309:
5277:
5235:
5174:
5142:
5110:
5068:
5007:
4975:
4943:
4901:
4840:
4808:
4776:
4734:
4673:
4641:
4609:
4567:
4506:
4474:
4445:
4403:
4325:
4145:
4116:
4077:
4016:
3977:
3938:
3858:
3807:
3708:
3654:
3271:
3217:
2674:{\displaystyle r=\rho (\mathbf {M} ^{-1}\mathbf {N} )}
2622:{\displaystyle \mathbf {A} =\mathbf {M} -\mathbf {N} }
1474:
1347:
1328:
1177:
1079:{\displaystyle \mathbf {A} =\mathbf {L} +\mathbf {U} }
405:
332:
142:
87:
10128:(2nd ed.). Pearson Education, Inc. p. 109.
9605:
8802:
8774:
8746:
8718:
8341:
7955:
7928:
7875:
7573:
7543:
7497:
7470:
7443:
7397:
7370:
7341:
7321:
7307:
is not guaranteed and, in this case, will not occur.
7236:
7214:
6837:
6771:
6736:
6711:
6689:
6667:
6389:
6261:
6239:
6122:
6100:
6078:
6056:
5967:
5889:
5793:
5686:
5671:{\displaystyle \mathbf {A} \mathbf {x} =\mathbf {b} }
5649:
5557:
4368:
4298:
4263:
4238:
4216:
4194:
3913:
3782:
3760:
3640:
3620:
3598:
3576:
3487:
3409:
3313:
3203:
3190:{\displaystyle \mathbf {A} \mathbf {x} =\mathbf {b} }
3168:
2945:
2923:
2888:
2843:
2763:
2728:
2691:
2635:
2597:
2575:
2542:
2514:
2489:
2455:
2401:
2366:
2325:
2284:
2238:
2182:
1925:
1901:
1860:
1838:
1728:
1705:
1681:
1512:
1498:
1161:
1139:
1112:
1092:
1054:
1032:
1006:
980:
897:
863:
841:
800:
778:
758:
723:
673:
627:
605:
570:
548:
526:
504:
482:
128:
8891:
3157:
2232:
that uses the most recently calculated iteration of
1503:
The system of linear equations may be rewritten as:
3145:-loop) check if convergence is reached
9860:
8815:
8787:
8759:
8731:
8697:
8321:
7941:
7914:
7859:
7556:
7529:
7483:
7456:
7429:
7383:
7356:
7327:
7299:
7222:
7198:
6821:
6757:
6719:
6697:
6675:
6651:
6373:
6247:
6223:
6108:
6086:
6064:
6039:
5950:
5875:
5777:
5670:
5620:
5541:
4351:
4284:
4246:
4224:
4202:
4178:
3897:
3768:
3744:
3626:
3606:
3584:
3559:
3470:
3395:
3297:
3189:
2957:
2931:
2909:
2874:
2829:
2749:
2714:
2673:
2621:
2583:
2550:
2522:
2497:
2467:
2419:
2387:
2352:
2311:
2251:
2224:
2166:
1907:
1887:
1846:
1817:
1713:
1689:
1663:
1488:
1145:
1125:
1098:
1078:
1040:
1014:
988:
966:
875:
849:
827:
786:
764:
744:
709:
659:
613:
591:
556:
534:
512:
490:
466:
109:{\textstyle \mathbf {A} \mathbf {x} =\mathbf {b} }
108:
10600:
10341:This article incorporates text from the article
10435:
10395:Gauss Siedel Iteration from www.geocities.com
10297:
10268:
10211:
10199:
10187:
6072:into the sum of a lower triangular component
3592:into the sum of a lower triangular component
2715:{\displaystyle \mathbf {M} ^{-1}\mathbf {N} }
10449:
10442:
10428:
10309:(3rd ed.), Baltimore: Johns Hopkins,
2225:{\displaystyle \sum _{j\neq i}a_{ij}x_{j}}
1895:can be computed sequentially for each row
1086:. More specifically, the decomposition of
10391:From Holistic Numerical Methods Institute
10237:
891:The solution is obtained iteratively via
6094:and a strict upper triangular component
3614:and a strict upper triangular component
1827:
660:{\displaystyle \mathbf {x} _{i}^{(0)}=0}
10573:Basic Linear Algebra Subprograms (BLAS)
10223:
886:
10601:
10164:
5639:Another example for the matrix version
10423:
10324:
10284:
10148:
10123:
7437:; then solve the second equation for
10384:Gauss–Seidel from www.math-linux.com
9599:The following code uses the formula
8884:The exact solution of the system is
2353:{\displaystyle \mathbf {x} ^{(k+1)}}
2312:{\displaystyle \mathbf {x} ^{(k+1)}}
2278:In Gauss-Seidel, the computation of
1888:{\displaystyle \mathbf {x} ^{(k+1)}}
828:{\displaystyle \mathbf {x} ^{(k+1)}}
7311:An example for the equation version
2395:that have not been computed in the
13:
10325:Black, Noel & Moore, Shirley.
7530:{\displaystyle x_{3},\dots ,x_{n}}
7430:{\displaystyle x_{2},\dots ,x_{n}}
6758:{\displaystyle \mathbf {x} ^{(0)}}
4285:{\displaystyle \mathbf {x} ^{(0)}}
2910:{\displaystyle \mathbf {x} ^{(0)}}
2750:{\displaystyle \mathbf {x} ^{(k)}}
2434:, since they can have a very long
2388:{\displaystyle \mathbf {x} ^{(k)}}
1499:Why the matrix-based formula works
772:-th approximation or iteration of
745:{\displaystyle \mathbf {x} ^{(k)}}
592:{\displaystyle \mathbf {x} ^{(0)}}
14:
10625:
10359:
8892:An example using Python and NumPy
7915:{\displaystyle x_{1},x_{2},x_{3}}
3158:An example for the matrix version
2438:, and are thus most feasible for
33:method of successive displacement
10085:Iterative method: Linear systems
8991:"System of equations:"
7357:{\displaystyle \mathbf {x} _{0}}
7344:
7261:
7247:
7238:
7216:
7163:
7002:
6844:
6774:
6739:
6730:First of all, we have to choose
6713:
6691:
6669:
6534:
6395:
6264:
6241:
6175:
6124:
6102:
6080:
6058:
6030:
6016:
6007:
5995:
5981:
5969:
5944:
5924:
5918:
5892:
5854:
5848:
5840:
5823:
5796:
5739:
5688:
5664:
5656:
5651:
5582:
5568:
5559:
5374:
5207:
5040:
4873:
4706:
4539:
4375:
4301:
4266:
4240:
4218:
4196:
4061:
3919:
3785:
3762:
3696:
3642:
3600:
3578:
3550:
3536:
3527:
3515:
3501:
3489:
3464:
3444:
3438:
3412:
3374:
3368:
3360:
3343:
3316:
3259:
3205:
3183:
3175:
3170:
2925:
2891:
2868:
2854:
2845:
2823:
2803:
2797:
2771:
2765:
2731:
2708:
2694:
2664:
2650:
2615:
2607:
2599:
2577:
2544:
2516:
2491:
2369:
2328:
2287:
2275:, with an important difference:
1863:
1840:
1791:
1785:
1777:
1758:
1731:
1707:
1683:
1653:
1648:
1640:
1628:
1623:
1614:
1602:
1597:
1589:
1584:
1575:
1563:
1555:
1547:
1535:
1523:
1518:
1476:
1330:
1163:
1072:
1064:
1056:
1034:
1008:
982:
945:
939:
931:
905:
899:
843:
803:
780:
726:
630:
607:
573:
550:
528:
506:
484:
393:
320:
130:
102:
94:
89:
9791:
7335:equations and a starting point
6173:
6167:
6005:
5999:
5737:
5731:
5643:Another linear system shown as
3694:
3688:
3525:
3519:
3257:
3251:
2133:
1698:, using the previous value for
391:
318:
10614:Relaxation (iterative methods)
10262:
10217:
10205:
10193:
10181:
10158:
10142:
10117:
9778:
9772:
9725:
9713:
9628:
9616:
8660:
8651:
8637:
8631:
8538:
8532:
8518:
8504:
8418:
8404:
7174:
7168:
7013:
7007:
6855:
6849:
6785:
6779:
6750:
6744:
5935:
5929:
5909:
5897:
5870:
5865:
5859:
5836:
5813:
5801:
5385:
5379:
5218:
5212:
5051:
5045:
4884:
4878:
4717:
4711:
4550:
4544:
4386:
4380:
4312:
4306:
4277:
4271:
3455:
3449:
3429:
3417:
3390:
3385:
3379:
3356:
3333:
3321:
2902:
2896:
2814:
2808:
2788:
2776:
2742:
2736:
2668:
2645:
2478:
2414:
2402:
2380:
2374:
2345:
2333:
2304:
2292:
2120:
2114:
2056:
2044:
1948:
1936:
1880:
1868:
1802:
1796:
1748:
1736:
1559:
1543:
956:
950:
922:
910:
820:
808:
737:
731:
646:
640:
599:denotes the initial guess for
584:
578:
76:
1:
10278:
7491:just found and the remaining
6727:can be obtained iteratively.
4254:can be obtained iteratively.
2266:
7223:{\displaystyle \mathbf {A} }
6720:{\displaystyle \mathbf {x} }
6698:{\displaystyle \mathbf {c} }
6676:{\displaystyle \mathbf {T} }
6248:{\displaystyle \mathbf {L} }
6109:{\displaystyle \mathbf {U} }
6087:{\displaystyle \mathbf {L} }
6065:{\displaystyle \mathbf {A} }
4247:{\displaystyle \mathbf {x} }
4225:{\displaystyle \mathbf {c} }
4203:{\displaystyle \mathbf {T} }
3769:{\displaystyle \mathbf {L} }
3607:{\displaystyle \mathbf {L} }
3585:{\displaystyle \mathbf {A} }
2968:
2932:{\displaystyle \mathbf {M} }
2584:{\displaystyle \mathbf {A} }
2551:{\displaystyle \mathbf {A} }
2523:{\displaystyle \mathbf {A} }
2498:{\displaystyle \mathbf {A} }
2445:Gauss-Seidel is the same as
1847:{\displaystyle \mathbf {L} }
1714:{\displaystyle \mathbf {x} }
1690:{\displaystyle \mathbf {x} }
1041:{\displaystyle \mathbf {U} }
1015:{\displaystyle \mathbf {L} }
989:{\displaystyle \mathbf {A} }
850:{\displaystyle \mathbf {x} }
787:{\displaystyle \mathbf {x} }
614:{\displaystyle \mathbf {x} }
557:{\displaystyle \mathbf {x} }
535:{\displaystyle \mathbf {x} }
513:{\displaystyle \mathbf {b} }
491:{\displaystyle \mathbf {A} }
59:strictly diagonally dominant
7:
10372:Encyclopedia of Mathematics
10080:Gaussian belief propagation
10068:
8964:# initialize the RHS vector
3152:
2558:is strictly or irreducibly
710:{\displaystyle i=1,2,...,n}
10:
10630:
10486:System of linear equations
3011:until convergence
2447:successive over-relaxation
41:system of linear equations
10555:
10537:Cache-oblivious algorithm
10519:
10478:
10457:
10269:Golub & Van Loan 1996
10212:Golub & Van Loan 1996
10200:Golub & Van Loan 1996
10188:Golub & Van Loan 1996
10075:Conjugate gradient method
3162:A linear system shown as
2468:{\displaystyle \omega =1}
1024:strictly upper triangular
120:linear equations, where:
55:Philipp Ludwig von Seidel
10609:Numerical linear algebra
10588:General purpose software
10451:Numerical linear algebra
10110:
9869:
9538:
8898:
3002:Choose an initial guess
2882:for any starting vector
2591:into two parts. Suppose
43:. It is named after the
21:numerical linear algebra
10400:The Gauss-Seidel Method
10165:Seidel, Ludwig (1874).
10124:Sauer, Timothy (2006).
8922:# initialize the matrix
9862:
8817:
8789:
8761:
8733:
8699:
8323:
7943:
7916:
7861:
7558:
7531:
7485:
7458:
7431:
7385:
7358:
7329:
7301:
7224:
7200:
6823:
6759:
6721:
6699:
6677:
6653:
6375:
6249:
6225:
6110:
6088:
6066:
6041:
5952:
5877:
5779:
5672:
5622:
5543:
4353:
4286:
4248:
4226:
4204:
4180:
3899:
3770:
3746:
3628:
3608:
3586:
3561:
3472:
3397:
3299:
3191:
2959:
2958:{\displaystyle r<1}
2933:
2911:
2876:
2831:
2751:
2716:
2675:
2623:
2585:
2552:
2524:
2499:
2469:
2421:
2389:
2354:
2313:
2253:
2226:
2168:
2090:
2020:
1909:
1889:
1848:
1819:
1715:
1691:
1665:
1490:
1147:
1127:
1100:
1080:
1042:
1016:
990:
968:
877:
851:
829:
788:
766:
746:
711:
661:
615:
593:
558:
536:
514:
492:
468:
116:be a square system of
110:
10583:Specialized libraries
10496:Matrix multiplication
10491:Matrix decompositions
10327:"Gauss-Seidel Method"
10286:Gauss, Carl Friedrich
9863:
9536:Produces the output:
8818:
8816:{\displaystyle x_{4}}
8790:
8788:{\displaystyle x_{3}}
8762:
8760:{\displaystyle x_{2}}
8734:
8732:{\displaystyle x_{1}}
8700:
8324:
7944:
7942:{\displaystyle x_{4}}
7917:
7862:
7567:Consider an example:
7559:
7557:{\displaystyle x_{n}}
7532:
7486:
7484:{\displaystyle x_{1}}
7459:
7457:{\displaystyle x_{2}}
7432:
7386:
7384:{\displaystyle x_{1}}
7359:
7330:
7302:
7225:
7201:
6824:
6760:
6722:
6700:
6678:
6654:
6376:
6250:
6226:
6111:
6089:
6067:
6042:
5953:
5878:
5780:
5673:
5623:
5544:
4354:
4287:
4257:First of all, choose
4249:
4227:
4205:
4181:
3900:
3771:
3747:
3629:
3609:
3587:
3562:
3473:
3398:
3300:
3192:
2960:
2934:
2912:
2877:
2832:
2752:
2717:
2676:
2624:
2586:
2553:
2525:
2500:
2470:
2422:
2420:{\displaystyle (k+1)}
2390:
2355:
2319:uses the elements of
2314:
2254:
2252:{\displaystyle x_{j}}
2227:
2169:
2064:
1994:
1910:
1890:
1849:
1828:Element-based formula
1820:
1716:
1692:
1666:
1491:
1148:
1128:
1126:{\displaystyle L_{*}}
1101:
1081:
1043:
1017:
996:is decomposed into a
991:
969:
878:
852:
835:the approximation of
830:
789:
767:
747:
712:
662:
616:
594:
559:
537:
515:
493:
469:
111:
10303:Van Loan, Charles F.
10105:Richardson iteration
9603:
9541:System of equations:
8800:
8772:
8744:
8716:
8339:
7953:
7926:
7873:
7571:
7541:
7495:
7468:
7441:
7395:
7368:
7339:
7319:
7234:
7212:
6835:
6769:
6734:
6709:
6687:
6665:
6387:
6259:
6237:
6120:
6098:
6076:
6054:
5965:
5887:
5791:
5684:
5647:
5631:In fact, the matrix
5555:
4366:
4296:
4261:
4236:
4214:
4192:
3911:
3780:
3758:
3638:
3618:
3596:
3574:
3485:
3407:
3311:
3201:
3166:
2979:Gauss–Seidel method
2943:
2921:
2886:
2841:
2761:
2726:
2722:. Then the iterates
2689:
2633:
2629:is nonsingular. Let
2595:
2573:
2540:
2512:
2487:
2453:
2399:
2364:
2323:
2282:
2236:
2180:
1923:
1917:forward substitution
1899:
1858:
1836:
1726:
1703:
1679:
1510:
1159:
1137:
1110:
1090:
1052:
1030:
1004:
978:
895:
887:Matrix-based formula
861:
839:
798:
776:
756:
721:
671:
625:
603:
568:
546:
524:
502:
480:
126:
85:
51:Carl Friedrich Gauss
27:, also known as the
10470:Numerical stability
10343:Gauss-Seidel_method
10307:Matrix Computations
9782:
9729:
9632:
3105:-loop)
2939:is nonsingular and
2560:diagonally dominant
2124:
2060:
1952:
876:{\displaystyle k+1}
650:
25:Gauss–Seidel method
10349:that is under the
10126:Numerical Analysis
9858:
9856:
9762:
9748:
9703:
9689:
9606:
8813:
8785:
8757:
8729:
8695:
8693:
8319:
8317:
7939:
7912:
7857:
7855:
7554:
7537:; and continue to
7527:
7481:
7454:
7427:
7381:
7354:
7325:
7297:
7291:
7220:
7196:
7194:
7148:
7116:
7084:
7058:
6987:
6958:
6926:
6900:
6819:
6813:
6755:
6717:
6695:
6673:
6649:
6647:
6636:
6604:
6578:
6520:
6478:
6442:
6371:
6362:
6311:
6245:
6221:
6212:
6161:
6106:
6084:
6062:
6037:
5948:
5873:
5775:
5766:
5725:
5668:
5618:
5612:
5539:
5537:
5526:
5494:
5462:
5433:
5359:
5327:
5295:
5266:
5192:
5160:
5128:
5099:
5025:
4993:
4961:
4932:
4858:
4826:
4794:
4765:
4691:
4659:
4627:
4598:
4524:
4492:
4460:
4434:
4349:
4340:
4282:
4244:
4222:
4200:
4176:
4174:
4163:
4131:
4105:
4047:
4002:
3966:
3895:
3886:
3835:
3766:
3742:
3733:
3682:
3624:
3604:
3582:
3557:
3468:
3393:
3295:
3286:
3245:
3187:
2955:
2929:
2907:
2872:
2827:
2747:
2712:
2671:
2619:
2581:
2548:
2520:
2495:
2465:
2417:
2385:
2350:
2309:
2249:
2222:
2198:
2164:
2104:
2034:
1926:
1905:
1885:
1854:, the elements of
1844:
1815:
1711:
1687:
1661:
1659:
1486:
1482:
1480:
1471:
1463:
1336:
1334:
1325:
1317:
1143:
1123:
1096:
1076:
1038:
1012:
986:
964:
873:
847:
825:
784:
762:
742:
707:
657:
628:
611:
589:
554:
532:
510:
488:
464:
455:
382:
309:
106:
10596:
10595:
10316:978-0-8018-5414-9
10135:978-0-321-78367-7
9733:
9674:
9654:
8882:
8881:
7328:{\displaystyle n}
6171:
6003:
5787:Use the equation
5735:
3692:
3627:{\displaystyle U}
3523:
3307:Use the equation
3255:
2532:positive-definite
2183:
1974:
1908:{\displaystyle i}
1342:
1340:
1172:
1170:
1146:{\displaystyle U}
1099:{\displaystyle A}
974:where the matrix
765:{\displaystyle k}
67:positive definite
10621:
10506:Matrix splitting
10444:
10437:
10430:
10421:
10420:
10380:
10337:
10336:
10319:
10293:
10272:
10266:
10260:
10259:
10241:
10221:
10215:
10209:
10203:
10197:
10191:
10185:
10179:
10178:
10162:
10156:
10146:
10140:
10139:
10121:
10100:Matrix splitting
10064:
10061:
10058:
10055:
10052:
10049:
10046:
10043:
10040:
10037:
10034:
10031:
10028:
10025:
10022:
10019:
10016:
10013:
10010:
10007:
10004:
10001:
9998:
9995:
9992:
9989:
9986:
9983:
9980:
9977:
9974:
9971:
9968:
9965:
9962:
9959:
9956:
9953:
9950:
9947:
9944:
9941:
9938:
9935:
9932:
9929:
9926:
9923:
9920:
9917:
9914:
9911:
9908:
9905:
9902:
9899:
9896:
9893:
9890:
9887:
9883:
9880:
9877:
9873:
9867:
9865:
9864:
9859:
9857:
9787:
9783:
9781:
9770:
9761:
9760:
9747:
9728:
9711:
9702:
9701:
9688:
9670:
9669:
9655:
9653:
9652:
9637:
9631:
9614:
9590:
9587:
9584:
9581:
9578:
9575:
9572:
9569:
9566:
9563:
9560:
9557:
9554:
9551:
9548:
9545:
9542:
9532:
9529:
9526:
9523:
9520:
9517:
9514:
9511:
9508:
9505:
9502:
9499:
9496:
9493:
9490:
9487:
9484:
9481:
9478:
9475:
9472:
9469:
9466:
9463:
9460:
9457:
9454:
9453:"Solution:
9451:
9448:
9445:
9442:
9439:
9436:
9433:
9430:
9427:
9424:
9421:
9418:
9415:
9412:
9409:
9406:
9403:
9400:
9397:
9394:
9391:
9388:
9385:
9382:
9379:
9376:
9373:
9370:
9367:
9364:
9361:
9358:
9355:
9352:
9349:
9346:
9343:
9340:
9337:
9334:
9331:
9328:
9325:
9322:
9319:
9316:
9313:
9310:
9307:
9304:
9301:
9298:
9295:
9292:
9289:
9286:
9283:
9280:
9277:
9274:
9271:
9268:
9265:
9262:
9259:
9256:
9253:
9250:
9247:
9244:
9243:"Iteration
9241:
9238:
9235:
9232:
9229:
9226:
9223:
9220:
9217:
9214:
9211:
9208:
9205:
9202:
9199:
9196:
9193:
9190:
9187:
9184:
9181:
9178:
9175:
9172:
9169:
9166:
9163:
9160:
9157:
9154:
9151:
9148:
9145:
9142:
9139:
9136:
9133:
9130:
9127:
9124:
9121:
9118:
9115:
9112:
9109:
9106:
9103:
9100:
9097:
9094:
9091:
9088:
9085:
9082:
9079:
9076:
9073:
9070:
9067:
9064:
9061:
9058:
9055:
9052:
9049:
9046:
9043:
9040:
9037:
9034:
9031:
9028:
9025:
9022:
9019:
9016:
9013:
9010:
9007:
9004:
9001:
8998:
8995:
8992:
8989:
8986:
8983:
8980:
8977:
8974:
8971:
8968:
8965:
8962:
8959:
8956:
8953:
8950:
8947:
8944:
8941:
8938:
8935:
8932:
8929:
8926:
8923:
8920:
8917:
8914:
8911:
8908:
8905:
8902:
8887:
8822:
8820:
8819:
8814:
8812:
8811:
8794:
8792:
8791:
8786:
8784:
8783:
8766:
8764:
8763:
8758:
8756:
8755:
8738:
8736:
8735:
8730:
8728:
8727:
8710:
8709:
8704:
8702:
8701:
8696:
8694:
8681:
8667:
8644:
8617:
8616:
8576:
8559:
8545:
8525:
8514:
8493:
8492:
8467:
8453:
8439:
8425:
8414:
8396:
8395:
8370:
8355:
8354:
8334:
8328:
8326:
8325:
8320:
8318:
8311:
8297:
8292:
8291:
8276:
8271:
8270:
8248:
8247:
8228:
8214:
8209:
8208:
8193:
8188:
8187:
8172:
8167:
8166:
8147:
8146:
8127:
8113:
8108:
8107:
8089:
8084:
8083:
8068:
8063:
8062:
8046:
8045:
8026:
8012:
8007:
8006:
7991:
7986:
7985:
7969:
7968:
7948:
7946:
7945:
7940:
7938:
7937:
7921:
7919:
7918:
7913:
7911:
7910:
7898:
7897:
7885:
7884:
7866:
7864:
7863:
7858:
7856:
7844:
7843:
7826:
7825:
7811:
7810:
7797:
7779:
7778:
7764:
7763:
7746:
7745:
7731:
7730:
7703:
7702:
7685:
7684:
7670:
7669:
7652:
7651:
7625:
7623:
7622:
7605:
7604:
7590:
7589:
7563:
7561:
7560:
7555:
7553:
7552:
7536:
7534:
7533:
7528:
7526:
7525:
7507:
7506:
7490:
7488:
7487:
7482:
7480:
7479:
7463:
7461:
7460:
7455:
7453:
7452:
7436:
7434:
7433:
7428:
7426:
7425:
7407:
7406:
7390:
7388:
7387:
7382:
7380:
7379:
7363:
7361:
7360:
7355:
7353:
7352:
7347:
7334:
7332:
7331:
7326:
7306:
7304:
7303:
7298:
7296:
7295:
7264:
7259:
7258:
7250:
7241:
7229:
7227:
7226:
7221:
7219:
7205:
7203:
7202:
7197:
7195:
7178:
7177:
7166:
7153:
7152:
7121:
7120:
7089:
7088:
7063:
7062:
7017:
7016:
7005:
6992:
6991:
6963:
6962:
6931:
6930:
6905:
6904:
6859:
6858:
6847:
6831:Then calculate:
6828:
6826:
6825:
6820:
6818:
6817:
6789:
6788:
6777:
6764:
6762:
6761:
6756:
6754:
6753:
6742:
6726:
6724:
6723:
6718:
6716:
6704:
6702:
6701:
6696:
6694:
6682:
6680:
6679:
6674:
6672:
6658:
6656:
6655:
6650:
6648:
6641:
6640:
6609:
6608:
6583:
6582:
6537:
6525:
6524:
6483:
6482:
6447:
6446:
6398:
6380:
6378:
6377:
6372:
6367:
6366:
6325:
6324:
6316:
6315:
6276:
6275:
6267:
6254:
6252:
6251:
6246:
6244:
6230:
6228:
6227:
6222:
6217:
6216:
6178:
6172:
6169:
6166:
6165:
6127:
6115:
6113:
6112:
6107:
6105:
6093:
6091:
6090:
6085:
6083:
6071:
6069:
6068:
6063:
6061:
6046:
6044:
6043:
6038:
6033:
6028:
6027:
6019:
6010:
6004:
6001:
5998:
5993:
5992:
5984:
5972:
5957:
5955:
5954:
5949:
5947:
5939:
5938:
5927:
5921:
5913:
5912:
5895:
5882:
5880:
5879:
5874:
5869:
5868:
5857:
5851:
5843:
5835:
5834:
5826:
5817:
5816:
5799:
5784:
5782:
5781:
5776:
5771:
5770:
5742:
5736:
5733:
5730:
5729:
5691:
5677:
5675:
5674:
5669:
5667:
5659:
5654:
5634:
5627:
5625:
5624:
5619:
5617:
5616:
5585:
5580:
5579:
5571:
5562:
5548:
5546:
5545:
5540:
5538:
5531:
5530:
5499:
5498:
5467:
5466:
5438:
5437:
5389:
5388:
5377:
5364:
5363:
5332:
5331:
5300:
5299:
5271:
5270:
5222:
5221:
5210:
5197:
5196:
5165:
5164:
5133:
5132:
5104:
5103:
5055:
5054:
5043:
5030:
5029:
4998:
4997:
4966:
4965:
4937:
4936:
4888:
4887:
4876:
4863:
4862:
4831:
4830:
4799:
4798:
4770:
4769:
4721:
4720:
4709:
4696:
4695:
4664:
4663:
4632:
4631:
4603:
4602:
4554:
4553:
4542:
4529:
4528:
4497:
4496:
4465:
4464:
4439:
4438:
4390:
4389:
4378:
4362:Then calculate:
4358:
4356:
4355:
4350:
4345:
4344:
4316:
4315:
4304:
4291:
4289:
4288:
4283:
4281:
4280:
4269:
4253:
4251:
4250:
4245:
4243:
4231:
4229:
4228:
4223:
4221:
4209:
4207:
4206:
4201:
4199:
4185:
4183:
4182:
4177:
4175:
4168:
4167:
4136:
4135:
4110:
4109:
4064:
4052:
4051:
4007:
4006:
3971:
3970:
3922:
3904:
3902:
3901:
3896:
3891:
3890:
3849:
3848:
3840:
3839:
3797:
3796:
3788:
3775:
3773:
3772:
3767:
3765:
3751:
3749:
3748:
3743:
3738:
3737:
3699:
3693:
3690:
3687:
3686:
3645:
3633:
3631:
3630:
3625:
3613:
3611:
3610:
3605:
3603:
3591:
3589:
3588:
3583:
3581:
3566:
3564:
3563:
3558:
3553:
3548:
3547:
3539:
3530:
3524:
3521:
3518:
3513:
3512:
3504:
3492:
3477:
3475:
3474:
3469:
3467:
3459:
3458:
3447:
3441:
3433:
3432:
3415:
3402:
3400:
3399:
3394:
3389:
3388:
3377:
3371:
3363:
3355:
3354:
3346:
3337:
3336:
3319:
3304:
3302:
3301:
3296:
3291:
3290:
3262:
3256:
3253:
3250:
3249:
3208:
3196:
3194:
3193:
3188:
3186:
3178:
3173:
3137:
3094:
3036:
3007:
3000:
2964:
2962:
2961:
2956:
2938:
2936:
2935:
2930:
2928:
2916:
2914:
2913:
2908:
2906:
2905:
2894:
2881:
2879:
2878:
2873:
2871:
2866:
2865:
2857:
2848:
2836:
2834:
2833:
2828:
2826:
2818:
2817:
2806:
2800:
2792:
2791:
2774:
2768:
2756:
2754:
2753:
2748:
2746:
2745:
2734:
2721:
2719:
2718:
2713:
2711:
2706:
2705:
2697:
2680:
2678:
2677:
2672:
2667:
2662:
2661:
2653:
2628:
2626:
2625:
2620:
2618:
2610:
2602:
2590:
2588:
2587:
2582:
2580:
2557:
2555:
2554:
2549:
2547:
2529:
2527:
2526:
2521:
2519:
2504:
2502:
2501:
2496:
2494:
2474:
2472:
2471:
2466:
2426:
2424:
2423:
2418:
2394:
2392:
2391:
2386:
2384:
2383:
2372:
2359:
2357:
2356:
2351:
2349:
2348:
2331:
2318:
2316:
2315:
2310:
2308:
2307:
2290:
2258:
2256:
2255:
2250:
2248:
2247:
2231:
2229:
2228:
2223:
2221:
2220:
2211:
2210:
2197:
2173:
2171:
2170:
2165:
2129:
2125:
2123:
2112:
2103:
2102:
2089:
2084:
2059:
2042:
2033:
2032:
2019:
2008:
1990:
1989:
1975:
1973:
1972:
1957:
1951:
1934:
1914:
1912:
1911:
1906:
1894:
1892:
1891:
1886:
1884:
1883:
1866:
1853:
1851:
1850:
1845:
1843:
1824:
1822:
1821:
1816:
1811:
1807:
1806:
1805:
1794:
1788:
1780:
1770:
1769:
1761:
1752:
1751:
1734:
1720:
1718:
1717:
1712:
1710:
1696:
1694:
1693:
1688:
1686:
1670:
1668:
1667:
1662:
1660:
1656:
1651:
1643:
1631:
1626:
1617:
1605:
1600:
1592:
1587:
1578:
1566:
1558:
1550:
1538:
1526:
1521:
1495:
1493:
1492:
1487:
1481:
1479:
1472:
1467:
1416:
1415:
1384:
1383:
1364:
1363:
1335:
1333:
1326:
1321:
1314:
1313:
1294:
1293:
1279:
1278:
1230:
1229:
1218:
1217:
1189:
1188:
1166:
1152:
1150:
1149:
1144:
1132:
1130:
1129:
1124:
1122:
1121:
1105:
1103:
1102:
1097:
1085:
1083:
1082:
1077:
1075:
1067:
1059:
1047:
1045:
1044:
1039:
1037:
1021:
1019:
1018:
1013:
1011:
998:lower triangular
995:
993:
992:
987:
985:
973:
971:
970:
965:
960:
959:
948:
942:
934:
926:
925:
908:
902:
883:-th) iteration.
882:
880:
879:
874:
857:at the next (or
856:
854:
853:
848:
846:
834:
832:
831:
826:
824:
823:
806:
793:
791:
790:
785:
783:
771:
769:
768:
763:
751:
749:
748:
743:
741:
740:
729:
716:
714:
713:
708:
666:
664:
663:
658:
649:
638:
633:
620:
618:
617:
612:
610:
598:
596:
595:
590:
588:
587:
576:
563:
561:
560:
555:
553:
541:
539:
538:
533:
531:
519:
517:
516:
511:
509:
497:
495:
494:
489:
487:
473:
471:
470:
465:
460:
459:
452:
451:
431:
430:
417:
416:
396:
387:
386:
379:
378:
358:
357:
344:
343:
323:
314:
313:
306:
305:
286:
285:
271:
270:
232:
231:
212:
211:
200:
199:
186:
185:
166:
165:
154:
153:
133:
119:
115:
113:
112:
107:
105:
97:
92:
39:used to solve a
37:iterative method
10629:
10628:
10624:
10623:
10622:
10620:
10619:
10618:
10599:
10598:
10597:
10592:
10551:
10547:Multiprocessing
10515:
10511:Sparse problems
10474:
10453:
10448:
10367:"Seidel method"
10365:
10362:
10356:
10317:
10281:
10276:
10275:
10267:
10263:
10248:10.1137/1037008
10222:
10218:
10210:
10206:
10198:
10194:
10186:
10182:
10163:
10159:
10151:, p. 279;
10147:
10143:
10136:
10122:
10118:
10113:
10090:Kaczmarz method
10071:
10066:
10065:
10062:
10059:
10056:
10053:
10050:
10047:
10044:
10041:
10038:
10035:
10032:
10029:
10026:
10023:
10020:
10017:
10014:
10011:
10008:
10005:
10002:
9999:
9996:
9993:
9990:
9987:
9984:
9981:
9978:
9975:
9972:
9969:
9966:
9963:
9960:
9957:
9954:
9951:
9948:
9945:
9942:
9939:
9936:
9933:
9930:
9927:
9924:
9921:
9918:
9915:
9912:
9909:
9906:
9903:
9900:
9897:
9894:
9891:
9888:
9885:
9881:
9878:
9875:
9871:
9855:
9854:
9824:
9823:
9792:
9771:
9766:
9753:
9749:
9737:
9712:
9707:
9694:
9690:
9678:
9665:
9661:
9660:
9656:
9645:
9641:
9636:
9615:
9610:
9604:
9601:
9600:
9597:
9592:
9591:
9588:
9585:
9582:
9579:
9576:
9573:
9570:
9567:
9564:
9561:
9558:
9555:
9552:
9549:
9546:
9543:
9540:
9534:
9533:
9530:
9527:
9524:
9521:
9518:
9515:
9512:
9509:
9506:
9503:
9500:
9497:
9494:
9491:
9488:
9485:
9482:
9479:
9476:
9473:
9470:
9467:
9464:
9461:
9458:
9455:
9452:
9449:
9446:
9443:
9440:
9437:
9434:
9431:
9428:
9425:
9422:
9419:
9416:
9413:
9410:
9407:
9404:
9401:
9398:
9395:
9392:
9389:
9386:
9383:
9380:
9377:
9374:
9371:
9368:
9365:
9362:
9359:
9356:
9353:
9350:
9347:
9344:
9341:
9338:
9335:
9332:
9329:
9326:
9323:
9320:
9317:
9314:
9311:
9308:
9305:
9302:
9299:
9296:
9293:
9290:
9287:
9284:
9281:
9278:
9275:
9272:
9269:
9266:
9263:
9260:
9257:
9254:
9251:
9248:
9245:
9242:
9239:
9236:
9233:
9230:
9227:
9224:
9221:
9218:
9215:
9212:
9209:
9206:
9203:
9200:
9197:
9194:
9191:
9188:
9186:ITERATION_LIMIT
9185:
9182:
9179:
9176:
9173:
9170:
9167:
9164:
9161:
9158:
9155:
9152:
9149:
9146:
9143:
9140:
9137:
9134:
9131:
9128:
9125:
9122:
9119:
9116:
9113:
9110:
9107:
9105:" + "
9104:
9101:
9098:
9095:
9093:" = "
9092:
9089:
9086:
9083:
9080:
9077:
9074:
9071:
9068:
9065:
9062:
9059:
9056:
9053:
9050:
9047:
9044:
9041:
9038:
9035:
9032:
9029:
9026:
9023:
9020:
9017:
9014:
9011:
9008:
9005:
9002:
8999:
8996:
8993:
8990:
8987:
8984:
8981:
8978:
8975:
8972:
8969:
8966:
8963:
8960:
8957:
8954:
8951:
8948:
8945:
8942:
8939:
8936:
8933:
8930:
8927:
8924:
8921:
8918:
8915:
8913:ITERATION_LIMIT
8912:
8909:
8906:
8903:
8900:
8894:
8885:
8807:
8803:
8801:
8798:
8797:
8779:
8775:
8773:
8770:
8769:
8751:
8747:
8745:
8742:
8741:
8723:
8719:
8717:
8714:
8713:
8692:
8691:
8677:
8663:
8640:
8618:
8612:
8608:
8605:
8604:
8572:
8555:
8541:
8521:
8510:
8494:
8488:
8484:
8481:
8480:
8463:
8449:
8435:
8421:
8410:
8397:
8391:
8387:
8384:
8383:
8366:
8356:
8350:
8346:
8342:
8340:
8337:
8336:
8332:
8316:
8315:
8307:
8293:
8287:
8283:
8272:
8266:
8262:
8249:
8243:
8239:
8236:
8235:
8224:
8210:
8204:
8200:
8189:
8183:
8179:
8168:
8162:
8158:
8148:
8142:
8138:
8135:
8134:
8123:
8109:
8103:
8099:
8085:
8079:
8075:
8064:
8058:
8054:
8047:
8041:
8037:
8034:
8033:
8022:
8008:
8002:
7998:
7987:
7981:
7977:
7970:
7964:
7960:
7956:
7954:
7951:
7950:
7933:
7929:
7927:
7924:
7923:
7906:
7902:
7893:
7889:
7880:
7876:
7874:
7871:
7870:
7854:
7853:
7845:
7839:
7835:
7827:
7821:
7817:
7812:
7806:
7802:
7795:
7794:
7780:
7774:
7770:
7765:
7759:
7755:
7747:
7741:
7737:
7732:
7726:
7722:
7716:
7715:
7704:
7698:
7694:
7686:
7680:
7676:
7671:
7665:
7661:
7653:
7647:
7643:
7637:
7636:
7624:
7618:
7614:
7606:
7600:
7596:
7591:
7585:
7581:
7574:
7572:
7569:
7568:
7548:
7544:
7542:
7539:
7538:
7521:
7517:
7502:
7498:
7496:
7493:
7492:
7475:
7471:
7469:
7466:
7465:
7448:
7444:
7442:
7439:
7438:
7421:
7417:
7402:
7398:
7396:
7393:
7392:
7375:
7371:
7369:
7366:
7365:
7348:
7343:
7342:
7340:
7337:
7336:
7320:
7317:
7316:
7313:
7290:
7289:
7283:
7282:
7269:
7268:
7260:
7251:
7246:
7245:
7237:
7235:
7232:
7231:
7215:
7213:
7210:
7209:
7193:
7192:
7179:
7167:
7162:
7161:
7158:
7157:
7147:
7146:
7137:
7136:
7126:
7125:
7115:
7114:
7105:
7104:
7094:
7093:
7083:
7082:
7076:
7075:
7065:
7064:
7057:
7056:
7051:
7045:
7044:
7036:
7026:
7025:
7018:
7006:
7001:
7000:
6997:
6996:
6986:
6985:
6979:
6978:
6968:
6967:
6957:
6956:
6947:
6946:
6936:
6935:
6925:
6924:
6918:
6917:
6907:
6906:
6899:
6898:
6893:
6887:
6886:
6878:
6868:
6867:
6860:
6848:
6843:
6842:
6838:
6836:
6833:
6832:
6812:
6811:
6805:
6804:
6794:
6793:
6778:
6773:
6772:
6770:
6767:
6766:
6743:
6738:
6737:
6735:
6732:
6731:
6712:
6710:
6707:
6706:
6690:
6688:
6685:
6684:
6668:
6666:
6663:
6662:
6646:
6645:
6635:
6634:
6625:
6624:
6614:
6613:
6603:
6602:
6596:
6595:
6585:
6584:
6577:
6576:
6571:
6562:
6561:
6556:
6546:
6545:
6538:
6533:
6530:
6529:
6519:
6518:
6513:
6507:
6506:
6498:
6488:
6487:
6477:
6476:
6471:
6465:
6464:
6459:
6449:
6448:
6441:
6440:
6435:
6426:
6425:
6420:
6410:
6409:
6399:
6394:
6390:
6388:
6385:
6384:
6361:
6360:
6355:
6346:
6345:
6340:
6330:
6329:
6317:
6310:
6309:
6304:
6298:
6297:
6292:
6282:
6281:
6280:
6268:
6263:
6262:
6260:
6257:
6256:
6240:
6238:
6235:
6234:
6233:The inverse of
6211:
6210:
6205:
6199:
6198:
6193:
6183:
6182:
6174:
6168:
6160:
6159:
6154:
6148:
6147:
6142:
6132:
6131:
6123:
6121:
6118:
6117:
6101:
6099:
6096:
6095:
6079:
6077:
6074:
6073:
6057:
6055:
6052:
6051:
6029:
6020:
6015:
6014:
6006:
6000:
5994:
5985:
5980:
5979:
5968:
5966:
5963:
5962:
5943:
5928:
5923:
5922:
5917:
5896:
5891:
5890:
5888:
5885:
5884:
5858:
5853:
5852:
5847:
5839:
5827:
5822:
5821:
5800:
5795:
5794:
5792:
5789:
5788:
5765:
5764:
5758:
5757:
5747:
5746:
5738:
5732:
5724:
5723:
5718:
5712:
5711:
5706:
5696:
5695:
5687:
5685:
5682:
5681:
5663:
5655:
5650:
5648:
5645:
5644:
5641:
5632:
5611:
5610:
5601:
5600:
5590:
5589:
5581:
5572:
5567:
5566:
5558:
5556:
5553:
5552:
5536:
5535:
5525:
5524:
5515:
5514:
5504:
5503:
5493:
5492:
5483:
5482:
5472:
5471:
5461:
5460:
5451:
5450:
5440:
5439:
5432:
5431:
5423:
5417:
5416:
5408:
5398:
5397:
5390:
5378:
5373:
5372:
5369:
5368:
5358:
5357:
5348:
5347:
5337:
5336:
5326:
5325:
5316:
5315:
5305:
5304:
5294:
5293:
5284:
5283:
5273:
5272:
5265:
5264:
5256:
5250:
5249:
5241:
5231:
5230:
5223:
5211:
5206:
5205:
5202:
5201:
5191:
5190:
5181:
5180:
5170:
5169:
5159:
5158:
5149:
5148:
5138:
5137:
5127:
5126:
5117:
5116:
5106:
5105:
5098:
5097:
5089:
5083:
5082:
5074:
5064:
5063:
5056:
5044:
5039:
5038:
5035:
5034:
5024:
5023:
5014:
5013:
5003:
5002:
4992:
4991:
4982:
4981:
4971:
4970:
4960:
4959:
4950:
4949:
4939:
4938:
4931:
4930:
4922:
4916:
4915:
4907:
4897:
4896:
4889:
4877:
4872:
4871:
4868:
4867:
4857:
4856:
4847:
4846:
4836:
4835:
4825:
4824:
4815:
4814:
4804:
4803:
4793:
4792:
4783:
4782:
4772:
4771:
4764:
4763:
4755:
4749:
4748:
4740:
4730:
4729:
4722:
4710:
4705:
4704:
4701:
4700:
4690:
4689:
4680:
4679:
4669:
4668:
4658:
4657:
4648:
4647:
4637:
4636:
4626:
4625:
4616:
4615:
4605:
4604:
4597:
4596:
4588:
4582:
4581:
4573:
4563:
4562:
4555:
4543:
4538:
4537:
4534:
4533:
4523:
4522:
4513:
4512:
4502:
4501:
4491:
4490:
4481:
4480:
4470:
4469:
4459:
4458:
4452:
4451:
4441:
4440:
4433:
4432:
4424:
4418:
4417:
4409:
4399:
4398:
4391:
4379:
4374:
4373:
4369:
4367:
4364:
4363:
4339:
4338:
4332:
4331:
4321:
4320:
4305:
4300:
4299:
4297:
4294:
4293:
4270:
4265:
4264:
4262:
4259:
4258:
4239:
4237:
4234:
4233:
4217:
4215:
4212:
4211:
4195:
4193:
4190:
4189:
4173:
4172:
4162:
4161:
4152:
4151:
4141:
4140:
4130:
4129:
4123:
4122:
4112:
4111:
4104:
4103:
4095:
4089:
4088:
4083:
4073:
4072:
4065:
4060:
4057:
4056:
4046:
4045:
4037:
4031:
4030:
4022:
4012:
4011:
4001:
4000:
3995:
3989:
3988:
3983:
3973:
3972:
3965:
3964:
3956:
3950:
3949:
3944:
3934:
3933:
3923:
3918:
3914:
3912:
3909:
3908:
3885:
3884:
3876:
3870:
3869:
3864:
3854:
3853:
3841:
3834:
3833:
3825:
3819:
3818:
3813:
3803:
3802:
3801:
3789:
3784:
3783:
3781:
3778:
3777:
3761:
3759:
3756:
3755:
3754:The inverse of
3732:
3731:
3726:
3720:
3719:
3714:
3704:
3703:
3695:
3689:
3681:
3680:
3672:
3666:
3665:
3660:
3650:
3649:
3641:
3639:
3636:
3635:
3619:
3616:
3615:
3599:
3597:
3594:
3593:
3577:
3575:
3572:
3571:
3549:
3540:
3535:
3534:
3526:
3520:
3514:
3505:
3500:
3499:
3488:
3486:
3483:
3482:
3463:
3448:
3443:
3442:
3437:
3416:
3411:
3410:
3408:
3405:
3404:
3378:
3373:
3372:
3367:
3359:
3347:
3342:
3341:
3320:
3315:
3314:
3312:
3309:
3308:
3285:
3284:
3278:
3277:
3267:
3266:
3258:
3252:
3244:
3243:
3235:
3229:
3228:
3223:
3213:
3212:
3204:
3202:
3199:
3198:
3182:
3174:
3169:
3167:
3164:
3163:
3160:
3155:
3150:
3144:
3136:
3123:
3114:
3106:
3104:
3093:
3085:
3069:
3065:
3061:
3052:
3042:
3031:
3027:
3017:
3006:to the solution
3001:
2993:
2992:
2988:
2971:
2944:
2941:
2940:
2924:
2922:
2919:
2918:
2895:
2890:
2889:
2887:
2884:
2883:
2867:
2858:
2853:
2852:
2844:
2842:
2839:
2838:
2822:
2807:
2802:
2801:
2796:
2775:
2770:
2769:
2764:
2762:
2759:
2758:
2735:
2730:
2729:
2727:
2724:
2723:
2707:
2698:
2693:
2692:
2690:
2687:
2686:
2683:spectral radius
2663:
2654:
2649:
2648:
2634:
2631:
2630:
2614:
2606:
2598:
2596:
2593:
2592:
2576:
2574:
2571:
2570:
2543:
2541:
2538:
2537:
2515:
2513:
2510:
2509:
2490:
2488:
2485:
2484:
2481:
2454:
2451:
2450:
2440:sparse matrices
2400:
2397:
2396:
2373:
2368:
2367:
2365:
2362:
2361:
2332:
2327:
2326:
2324:
2321:
2320:
2291:
2286:
2285:
2283:
2280:
2279:
2269:
2243:
2239:
2237:
2234:
2233:
2216:
2212:
2203:
2199:
2187:
2181:
2178:
2177:
2113:
2108:
2095:
2091:
2085:
2068:
2043:
2038:
2025:
2021:
2009:
1998:
1985:
1981:
1980:
1976:
1965:
1961:
1956:
1935:
1930:
1924:
1921:
1920:
1900:
1897:
1896:
1867:
1862:
1861:
1859:
1856:
1855:
1839:
1837:
1834:
1833:
1830:
1795:
1790:
1789:
1784:
1776:
1775:
1771:
1762:
1757:
1756:
1735:
1730:
1729:
1727:
1724:
1723:
1706:
1704:
1701:
1700:
1682:
1680:
1677:
1676:
1658:
1657:
1652:
1647:
1639:
1632:
1627:
1622:
1619:
1618:
1613:
1606:
1601:
1596:
1588:
1583:
1580:
1579:
1574:
1567:
1562:
1554:
1546:
1540:
1539:
1534:
1527:
1522:
1517:
1513:
1511:
1508:
1507:
1501:
1475:
1473:
1462:
1461:
1456:
1451:
1446:
1440:
1439:
1434:
1429:
1424:
1418:
1417:
1408:
1404:
1402:
1397:
1392:
1386:
1385:
1376:
1372:
1370:
1365:
1359:
1355:
1353:
1343:
1341:
1329:
1327:
1316:
1315:
1306:
1302:
1300:
1295:
1286:
1282:
1280:
1271:
1267:
1264:
1263:
1258:
1253:
1248:
1242:
1241:
1236:
1231:
1225:
1221:
1219:
1213:
1209:
1206:
1205:
1200:
1195:
1190:
1184:
1180:
1173:
1171:
1162:
1160:
1157:
1156:
1138:
1135:
1134:
1117:
1113:
1111:
1108:
1107:
1091:
1088:
1087:
1071:
1063:
1055:
1053:
1050:
1049:
1033:
1031:
1028:
1027:
1007:
1005:
1002:
1001:
981:
979:
976:
975:
949:
944:
943:
938:
930:
909:
904:
903:
898:
896:
893:
892:
889:
862:
859:
858:
842:
840:
837:
836:
807:
802:
801:
799:
796:
795:
779:
777:
774:
773:
757:
754:
753:
730:
725:
724:
722:
719:
718:
672:
669:
668:
639:
634:
629:
626:
623:
622:
606:
604:
601:
600:
577:
572:
571:
569:
566:
565:
549:
547:
544:
543:
527:
525:
522:
521:
520:are known, and
505:
503:
500:
499:
483:
481:
478:
477:
454:
453:
447:
443:
440:
439:
433:
432:
426:
422:
419:
418:
412:
408:
401:
400:
392:
381:
380:
374:
370:
367:
366:
360:
359:
353:
349:
346:
345:
339:
335:
328:
327:
319:
308:
307:
298:
294:
292:
287:
278:
274:
272:
263:
259:
256:
255:
250:
245:
240:
234:
233:
224:
220:
218:
213:
207:
203:
201:
195:
191:
188:
187:
178:
174:
172:
167:
161:
157:
155:
149:
145:
138:
137:
129:
127:
124:
123:
117:
101:
93:
88:
86:
83:
82:
79:
29:Liebmann method
17:
12:
11:
5:
10627:
10617:
10616:
10611:
10594:
10593:
10591:
10590:
10585:
10580:
10575:
10570:
10565:
10559:
10557:
10553:
10552:
10550:
10549:
10544:
10539:
10534:
10529:
10523:
10521:
10517:
10516:
10514:
10513:
10508:
10503:
10493:
10488:
10482:
10480:
10476:
10475:
10473:
10472:
10467:
10465:Floating point
10461:
10459:
10455:
10454:
10447:
10446:
10439:
10432:
10424:
10418:
10417:
10415:C code example
10412:
10407:
10402:
10397:
10392:
10386:
10381:
10361:
10360:External links
10358:
10339:
10338:
10321:
10315:
10299:Golub, Gene H.
10295:
10280:
10277:
10274:
10273:
10261:
10239:10.1.1.26.5207
10216:
10204:
10202:, eqn (10.1.3)
10192:
10190:, p. 511.
10180:
10157:
10141:
10134:
10115:
10114:
10112:
10109:
10108:
10107:
10102:
10097:
10087:
10082:
10077:
10070:
10067:
9884:A, b, x, iters
9870:
9853:
9850:
9847:
9844:
9841:
9838:
9835:
9832:
9829:
9826:
9825:
9822:
9819:
9816:
9813:
9810:
9807:
9804:
9801:
9798:
9795:
9794:
9790:
9786:
9780:
9777:
9774:
9769:
9765:
9759:
9756:
9752:
9746:
9743:
9740:
9736:
9732:
9727:
9724:
9721:
9718:
9715:
9710:
9706:
9700:
9697:
9693:
9687:
9684:
9681:
9677:
9673:
9668:
9664:
9659:
9651:
9648:
9644:
9640:
9635:
9630:
9627:
9624:
9621:
9618:
9613:
9609:
9596:
9593:
9583:Iteration 10:
9539:
8899:
8893:
8890:
8880:
8879:
8876:
8873:
8870:
8866:
8865:
8862:
8859:
8856:
8852:
8851:
8848:
8845:
8842:
8838:
8837:
8834:
8831:
8828:
8824:
8823:
8810:
8806:
8795:
8782:
8778:
8767:
8754:
8750:
8739:
8726:
8722:
8690:
8687:
8684:
8680:
8676:
8673:
8670:
8666:
8662:
8659:
8656:
8653:
8650:
8647:
8643:
8639:
8636:
8633:
8630:
8627:
8624:
8621:
8619:
8615:
8611:
8607:
8606:
8603:
8600:
8597:
8594:
8591:
8588:
8585:
8582:
8579:
8575:
8571:
8568:
8565:
8562:
8558:
8554:
8551:
8548:
8544:
8540:
8537:
8534:
8531:
8528:
8524:
8520:
8517:
8513:
8509:
8506:
8503:
8500:
8497:
8495:
8491:
8487:
8483:
8482:
8479:
8476:
8473:
8470:
8466:
8462:
8459:
8456:
8452:
8448:
8445:
8442:
8438:
8434:
8431:
8428:
8424:
8420:
8417:
8413:
8409:
8406:
8403:
8400:
8398:
8394:
8390:
8386:
8385:
8382:
8379:
8376:
8373:
8369:
8365:
8362:
8359:
8357:
8353:
8349:
8345:
8344:
8314:
8310:
8306:
8303:
8300:
8296:
8290:
8286:
8282:
8279:
8275:
8269:
8265:
8261:
8258:
8255:
8252:
8250:
8246:
8242:
8238:
8237:
8234:
8231:
8227:
8223:
8220:
8217:
8213:
8207:
8203:
8199:
8196:
8192:
8186:
8182:
8178:
8175:
8171:
8165:
8161:
8157:
8154:
8151:
8149:
8145:
8141:
8137:
8136:
8133:
8130:
8126:
8122:
8119:
8116:
8112:
8106:
8102:
8098:
8095:
8092:
8088:
8082:
8078:
8074:
8071:
8067:
8061:
8057:
8053:
8050:
8048:
8044:
8040:
8036:
8035:
8032:
8029:
8025:
8021:
8018:
8015:
8011:
8005:
8001:
7997:
7994:
7990:
7984:
7980:
7976:
7973:
7971:
7967:
7963:
7959:
7958:
7936:
7932:
7909:
7905:
7901:
7896:
7892:
7888:
7883:
7879:
7852:
7849:
7846:
7842:
7838:
7834:
7831:
7828:
7824:
7820:
7816:
7813:
7809:
7805:
7801:
7798:
7796:
7793:
7790:
7787:
7784:
7781:
7777:
7773:
7769:
7766:
7762:
7758:
7754:
7751:
7748:
7744:
7740:
7736:
7733:
7729:
7725:
7721:
7718:
7717:
7714:
7711:
7708:
7705:
7701:
7697:
7693:
7690:
7687:
7683:
7679:
7675:
7672:
7668:
7664:
7660:
7657:
7654:
7650:
7646:
7642:
7639:
7638:
7635:
7632:
7629:
7626:
7621:
7617:
7613:
7610:
7607:
7603:
7599:
7595:
7592:
7588:
7584:
7580:
7577:
7576:
7551:
7547:
7524:
7520:
7516:
7513:
7510:
7505:
7501:
7478:
7474:
7451:
7447:
7424:
7420:
7416:
7413:
7410:
7405:
7401:
7378:
7374:
7351:
7346:
7324:
7315:Suppose given
7312:
7309:
7294:
7288:
7285:
7284:
7281:
7278:
7275:
7274:
7272:
7267:
7263:
7257:
7254:
7249:
7244:
7240:
7218:
7191:
7188:
7185:
7182:
7180:
7176:
7173:
7170:
7165:
7160:
7159:
7156:
7151:
7145:
7142:
7139:
7138:
7135:
7132:
7131:
7129:
7124:
7119:
7113:
7110:
7107:
7106:
7103:
7100:
7099:
7097:
7092:
7087:
7081:
7078:
7077:
7074:
7071:
7070:
7068:
7061:
7055:
7052:
7050:
7047:
7046:
7043:
7040:
7037:
7035:
7032:
7031:
7029:
7024:
7021:
7019:
7015:
7012:
7009:
7004:
6999:
6998:
6995:
6990:
6984:
6981:
6980:
6977:
6974:
6973:
6971:
6966:
6961:
6955:
6952:
6949:
6948:
6945:
6942:
6941:
6939:
6934:
6929:
6923:
6920:
6919:
6916:
6913:
6912:
6910:
6903:
6897:
6894:
6892:
6889:
6888:
6885:
6882:
6879:
6877:
6874:
6873:
6871:
6866:
6863:
6861:
6857:
6854:
6851:
6846:
6841:
6840:
6816:
6810:
6807:
6806:
6803:
6800:
6799:
6797:
6792:
6787:
6784:
6781:
6776:
6765:, for example
6752:
6749:
6746:
6741:
6715:
6693:
6671:
6644:
6639:
6633:
6630:
6627:
6626:
6623:
6620:
6619:
6617:
6612:
6607:
6601:
6598:
6597:
6594:
6591:
6590:
6588:
6581:
6575:
6572:
6570:
6567:
6564:
6563:
6560:
6557:
6555:
6552:
6551:
6549:
6544:
6541:
6539:
6536:
6532:
6531:
6528:
6523:
6517:
6514:
6512:
6509:
6508:
6505:
6502:
6499:
6497:
6494:
6493:
6491:
6486:
6481:
6475:
6472:
6470:
6467:
6466:
6463:
6460:
6458:
6455:
6454:
6452:
6445:
6439:
6436:
6434:
6431:
6428:
6427:
6424:
6421:
6419:
6416:
6415:
6413:
6408:
6405:
6402:
6400:
6397:
6393:
6392:
6370:
6365:
6359:
6356:
6354:
6351:
6348:
6347:
6344:
6341:
6339:
6336:
6335:
6333:
6328:
6323:
6320:
6314:
6308:
6305:
6303:
6300:
6299:
6296:
6293:
6291:
6288:
6287:
6285:
6279:
6274:
6271:
6266:
6243:
6220:
6215:
6209:
6206:
6204:
6201:
6200:
6197:
6194:
6192:
6189:
6188:
6186:
6181:
6177:
6164:
6158:
6155:
6153:
6150:
6149:
6146:
6143:
6141:
6138:
6137:
6135:
6130:
6126:
6104:
6082:
6060:
6048:
6047:
6036:
6032:
6026:
6023:
6018:
6013:
6009:
5997:
5991:
5988:
5983:
5978:
5975:
5971:
5946:
5942:
5937:
5934:
5931:
5926:
5920:
5916:
5911:
5908:
5905:
5902:
5899:
5894:
5872:
5867:
5864:
5861:
5856:
5850:
5846:
5842:
5838:
5833:
5830:
5825:
5820:
5815:
5812:
5809:
5806:
5803:
5798:
5774:
5769:
5763:
5760:
5759:
5756:
5753:
5752:
5750:
5745:
5741:
5728:
5722:
5719:
5717:
5714:
5713:
5710:
5707:
5705:
5702:
5701:
5699:
5694:
5690:
5666:
5662:
5658:
5653:
5640:
5637:
5615:
5609:
5606:
5603:
5602:
5599:
5596:
5595:
5593:
5588:
5584:
5578:
5575:
5570:
5565:
5561:
5534:
5529:
5523:
5520:
5517:
5516:
5513:
5510:
5509:
5507:
5502:
5497:
5491:
5488:
5485:
5484:
5481:
5478:
5477:
5475:
5470:
5465:
5459:
5456:
5453:
5452:
5449:
5446:
5445:
5443:
5436:
5430:
5427:
5424:
5422:
5419:
5418:
5415:
5412:
5409:
5407:
5404:
5403:
5401:
5396:
5393:
5391:
5387:
5384:
5381:
5376:
5371:
5370:
5367:
5362:
5356:
5353:
5350:
5349:
5346:
5343:
5342:
5340:
5335:
5330:
5324:
5321:
5318:
5317:
5314:
5311:
5310:
5308:
5303:
5298:
5292:
5289:
5286:
5285:
5282:
5279:
5278:
5276:
5269:
5263:
5260:
5257:
5255:
5252:
5251:
5248:
5245:
5242:
5240:
5237:
5236:
5234:
5229:
5226:
5224:
5220:
5217:
5214:
5209:
5204:
5203:
5200:
5195:
5189:
5186:
5183:
5182:
5179:
5176:
5175:
5173:
5168:
5163:
5157:
5154:
5151:
5150:
5147:
5144:
5143:
5141:
5136:
5131:
5125:
5122:
5119:
5118:
5115:
5112:
5111:
5109:
5102:
5096:
5093:
5090:
5088:
5085:
5084:
5081:
5078:
5075:
5073:
5070:
5069:
5067:
5062:
5059:
5057:
5053:
5050:
5047:
5042:
5037:
5036:
5033:
5028:
5022:
5019:
5016:
5015:
5012:
5009:
5008:
5006:
5001:
4996:
4990:
4987:
4984:
4983:
4980:
4977:
4976:
4974:
4969:
4964:
4958:
4955:
4952:
4951:
4948:
4945:
4944:
4942:
4935:
4929:
4926:
4923:
4921:
4918:
4917:
4914:
4911:
4908:
4906:
4903:
4902:
4900:
4895:
4892:
4890:
4886:
4883:
4880:
4875:
4870:
4869:
4866:
4861:
4855:
4852:
4849:
4848:
4845:
4842:
4841:
4839:
4834:
4829:
4823:
4820:
4817:
4816:
4813:
4810:
4809:
4807:
4802:
4797:
4791:
4788:
4785:
4784:
4781:
4778:
4777:
4775:
4768:
4762:
4759:
4756:
4754:
4751:
4750:
4747:
4744:
4741:
4739:
4736:
4735:
4733:
4728:
4725:
4723:
4719:
4716:
4713:
4708:
4703:
4702:
4699:
4694:
4688:
4685:
4682:
4681:
4678:
4675:
4674:
4672:
4667:
4662:
4656:
4653:
4650:
4649:
4646:
4643:
4642:
4640:
4635:
4630:
4624:
4621:
4618:
4617:
4614:
4611:
4610:
4608:
4601:
4595:
4592:
4589:
4587:
4584:
4583:
4580:
4577:
4574:
4572:
4569:
4568:
4566:
4561:
4558:
4556:
4552:
4549:
4546:
4541:
4536:
4535:
4532:
4527:
4521:
4518:
4515:
4514:
4511:
4508:
4507:
4505:
4500:
4495:
4489:
4486:
4483:
4482:
4479:
4476:
4475:
4473:
4468:
4463:
4457:
4454:
4453:
4450:
4447:
4446:
4444:
4437:
4431:
4428:
4425:
4423:
4420:
4419:
4416:
4413:
4410:
4408:
4405:
4404:
4402:
4397:
4394:
4392:
4388:
4385:
4382:
4377:
4372:
4371:
4348:
4343:
4337:
4334:
4333:
4330:
4327:
4326:
4324:
4319:
4314:
4311:
4308:
4303:
4292:, for example
4279:
4276:
4273:
4268:
4242:
4220:
4198:
4171:
4166:
4160:
4157:
4154:
4153:
4150:
4147:
4146:
4144:
4139:
4134:
4128:
4125:
4124:
4121:
4118:
4117:
4115:
4108:
4102:
4099:
4096:
4094:
4091:
4090:
4087:
4084:
4082:
4079:
4078:
4076:
4071:
4068:
4066:
4063:
4059:
4058:
4055:
4050:
4044:
4041:
4038:
4036:
4033:
4032:
4029:
4026:
4023:
4021:
4018:
4017:
4015:
4010:
4005:
3999:
3996:
3994:
3991:
3990:
3987:
3984:
3982:
3979:
3978:
3976:
3969:
3963:
3960:
3957:
3955:
3952:
3951:
3948:
3945:
3943:
3940:
3939:
3937:
3932:
3929:
3926:
3924:
3921:
3917:
3916:
3894:
3889:
3883:
3880:
3877:
3875:
3872:
3871:
3868:
3865:
3863:
3860:
3859:
3857:
3852:
3847:
3844:
3838:
3832:
3829:
3826:
3824:
3821:
3820:
3817:
3814:
3812:
3809:
3808:
3806:
3800:
3795:
3792:
3787:
3764:
3741:
3736:
3730:
3727:
3725:
3722:
3721:
3718:
3715:
3713:
3710:
3709:
3707:
3702:
3698:
3685:
3679:
3676:
3673:
3671:
3668:
3667:
3664:
3661:
3659:
3656:
3655:
3653:
3648:
3644:
3623:
3602:
3580:
3568:
3567:
3556:
3552:
3546:
3543:
3538:
3533:
3529:
3517:
3511:
3508:
3503:
3498:
3495:
3491:
3466:
3462:
3457:
3454:
3451:
3446:
3440:
3436:
3431:
3428:
3425:
3422:
3419:
3414:
3392:
3387:
3384:
3381:
3376:
3370:
3366:
3362:
3358:
3353:
3350:
3345:
3340:
3335:
3332:
3329:
3326:
3323:
3318:
3294:
3289:
3283:
3280:
3279:
3276:
3273:
3272:
3270:
3265:
3261:
3248:
3242:
3239:
3236:
3234:
3231:
3230:
3227:
3224:
3222:
3219:
3218:
3216:
3211:
3207:
3185:
3181:
3177:
3172:
3159:
3156:
3154:
3151:
3142:
3132:
3119:
3110:
3102:
3089:
3081:
3063:
3059:
3050:
3040:
3025:
3015:
2990:
2986:
2975:
2970:
2967:
2954:
2951:
2948:
2927:
2904:
2901:
2898:
2893:
2870:
2864:
2861:
2856:
2851:
2847:
2825:
2821:
2816:
2813:
2810:
2805:
2799:
2795:
2790:
2787:
2784:
2781:
2778:
2773:
2767:
2744:
2741:
2738:
2733:
2710:
2704:
2701:
2696:
2670:
2666:
2660:
2657:
2652:
2647:
2644:
2641:
2638:
2617:
2613:
2609:
2605:
2601:
2579:
2564:
2563:
2546:
2535:
2518:
2493:
2480:
2477:
2464:
2461:
2458:
2416:
2413:
2410:
2407:
2404:
2382:
2379:
2376:
2371:
2347:
2344:
2341:
2338:
2335:
2330:
2306:
2303:
2300:
2297:
2294:
2289:
2268:
2265:
2246:
2242:
2219:
2215:
2209:
2206:
2202:
2196:
2193:
2190:
2186:
2163:
2160:
2157:
2154:
2151:
2148:
2145:
2142:
2139:
2136:
2132:
2128:
2122:
2119:
2116:
2111:
2107:
2101:
2098:
2094:
2088:
2083:
2080:
2077:
2074:
2071:
2067:
2063:
2058:
2055:
2052:
2049:
2046:
2041:
2037:
2031:
2028:
2024:
2018:
2015:
2012:
2007:
2004:
2001:
1997:
1993:
1988:
1984:
1979:
1971:
1968:
1964:
1960:
1955:
1950:
1947:
1944:
1941:
1938:
1933:
1929:
1904:
1882:
1879:
1876:
1873:
1870:
1865:
1842:
1829:
1826:
1814:
1810:
1804:
1801:
1798:
1793:
1787:
1783:
1779:
1774:
1768:
1765:
1760:
1755:
1750:
1747:
1744:
1741:
1738:
1733:
1709:
1685:
1672:
1671:
1655:
1650:
1646:
1642:
1638:
1635:
1633:
1630:
1625:
1621:
1620:
1616:
1612:
1609:
1607:
1604:
1599:
1595:
1591:
1586:
1582:
1581:
1577:
1573:
1570:
1568:
1565:
1561:
1557:
1553:
1549:
1545:
1542:
1541:
1537:
1533:
1530:
1528:
1525:
1520:
1516:
1515:
1500:
1497:
1485:
1478:
1470:
1466:
1460:
1457:
1455:
1452:
1450:
1447:
1445:
1442:
1441:
1438:
1435:
1433:
1430:
1428:
1425:
1423:
1420:
1419:
1414:
1411:
1407:
1403:
1401:
1398:
1396:
1393:
1391:
1388:
1387:
1382:
1379:
1375:
1371:
1369:
1366:
1362:
1358:
1354:
1352:
1349:
1348:
1346:
1339:
1332:
1324:
1320:
1312:
1309:
1305:
1301:
1299:
1296:
1292:
1289:
1285:
1281:
1277:
1274:
1270:
1266:
1265:
1262:
1259:
1257:
1254:
1252:
1249:
1247:
1244:
1243:
1240:
1237:
1235:
1232:
1228:
1224:
1220:
1216:
1212:
1208:
1207:
1204:
1201:
1199:
1196:
1194:
1191:
1187:
1183:
1179:
1178:
1176:
1169:
1165:
1142:
1120:
1116:
1095:
1074:
1070:
1066:
1062:
1058:
1036:
1010:
984:
963:
958:
955:
952:
947:
941:
937:
933:
929:
924:
921:
918:
915:
912:
907:
901:
888:
885:
872:
869:
866:
845:
822:
819:
816:
813:
810:
805:
782:
761:
739:
736:
733:
728:
706:
703:
700:
697:
694:
691:
688:
685:
682:
679:
676:
656:
653:
648:
645:
642:
637:
632:
609:
586:
583:
580:
575:
552:
530:
508:
486:
463:
458:
450:
446:
442:
441:
438:
435:
434:
429:
425:
421:
420:
415:
411:
407:
406:
404:
399:
395:
390:
385:
377:
373:
369:
368:
365:
362:
361:
356:
352:
348:
347:
342:
338:
334:
333:
331:
326:
322:
317:
312:
304:
301:
297:
293:
291:
288:
284:
281:
277:
273:
269:
266:
262:
258:
257:
254:
251:
249:
246:
244:
241:
239:
236:
235:
230:
227:
223:
219:
217:
214:
210:
206:
202:
198:
194:
190:
189:
184:
181:
177:
173:
171:
168:
164:
160:
156:
152:
148:
144:
143:
141:
136:
132:
104:
100:
96:
91:
78:
75:
48:mathematicians
15:
9:
6:
4:
3:
2:
10626:
10615:
10612:
10610:
10607:
10606:
10604:
10589:
10586:
10584:
10581:
10579:
10576:
10574:
10571:
10569:
10566:
10564:
10561:
10560:
10558:
10554:
10548:
10545:
10543:
10540:
10538:
10535:
10533:
10530:
10528:
10525:
10524:
10522:
10518:
10512:
10509:
10507:
10504:
10501:
10497:
10494:
10492:
10489:
10487:
10484:
10483:
10481:
10477:
10471:
10468:
10466:
10463:
10462:
10460:
10456:
10452:
10445:
10440:
10438:
10433:
10431:
10426:
10425:
10422:
10416:
10413:
10411:
10408:
10406:
10403:
10401:
10398:
10396:
10393:
10390:
10387:
10385:
10382:
10378:
10374:
10373:
10368:
10364:
10363:
10357:
10354:
10352:
10348:
10344:
10334:
10333:
10328:
10322:
10318:
10312:
10308:
10304:
10300:
10296:
10291:
10287:
10283:
10282:
10270:
10265:
10257:
10253:
10249:
10245:
10240:
10235:
10231:
10227:
10220:
10214:, Thm 10.1.2.
10213:
10208:
10201:
10196:
10189:
10184:
10176:
10173:(in German).
10172:
10168:
10161:
10154:
10150:
10145:
10137:
10131:
10127:
10120:
10116:
10106:
10103:
10101:
10098:
10095:
10091:
10088:
10086:
10083:
10081:
10078:
10076:
10073:
10072:
9868:
9851:
9848:
9845:
9842:
9839:
9836:
9833:
9830:
9827:
9820:
9817:
9814:
9811:
9808:
9805:
9802:
9799:
9796:
9788:
9784:
9775:
9767:
9763:
9757:
9754:
9750:
9744:
9741:
9738:
9734:
9730:
9722:
9719:
9716:
9708:
9704:
9698:
9695:
9691:
9685:
9682:
9679:
9675:
9671:
9666:
9662:
9657:
9649:
9646:
9642:
9638:
9633:
9625:
9622:
9619:
9611:
9607:
9580:Iteration 9:
9577:Iteration 8:
9574:Iteration 7:
9571:Iteration 6:
9568:Iteration 5:
9565:Iteration 4:
9562:Iteration 3:
9559:Iteration 2:
9556:Iteration 1:
9537:
9516:"Error:
8897:
8889:
8886:(1, 2, −1, 1)
8877:
8874:
8871:
8868:
8867:
8863:
8860:
8857:
8854:
8853:
8849:
8846:
8843:
8840:
8839:
8835:
8832:
8829:
8826:
8825:
8808:
8804:
8796:
8780:
8776:
8768:
8752:
8748:
8740:
8724:
8720:
8712:
8711:
8708:
8705:
8688:
8685:
8682:
8678:
8674:
8671:
8668:
8664:
8657:
8654:
8648:
8645:
8641:
8634:
8628:
8625:
8622:
8620:
8613:
8609:
8601:
8598:
8595:
8592:
8589:
8586:
8583:
8580:
8577:
8573:
8569:
8566:
8563:
8560:
8556:
8552:
8549:
8546:
8542:
8535:
8529:
8526:
8522:
8515:
8511:
8507:
8501:
8498:
8496:
8489:
8485:
8477:
8474:
8471:
8468:
8464:
8460:
8457:
8454:
8450:
8446:
8443:
8440:
8436:
8432:
8429:
8426:
8422:
8415:
8411:
8407:
8401:
8399:
8392:
8388:
8380:
8377:
8374:
8371:
8367:
8363:
8360:
8358:
8351:
8347:
8329:
8312:
8308:
8304:
8301:
8298:
8294:
8288:
8284:
8280:
8277:
8273:
8267:
8263:
8259:
8256:
8253:
8251:
8244:
8240:
8232:
8229:
8225:
8221:
8218:
8215:
8211:
8205:
8201:
8197:
8194:
8190:
8184:
8180:
8176:
8173:
8169:
8163:
8159:
8155:
8152:
8150:
8143:
8139:
8131:
8128:
8124:
8120:
8117:
8114:
8110:
8104:
8100:
8096:
8093:
8090:
8086:
8080:
8076:
8072:
8069:
8065:
8059:
8055:
8051:
8049:
8042:
8038:
8030:
8027:
8023:
8019:
8016:
8013:
8009:
8003:
7999:
7995:
7992:
7988:
7982:
7978:
7974:
7972:
7965:
7961:
7934:
7930:
7907:
7903:
7899:
7894:
7890:
7886:
7881:
7877:
7867:
7850:
7847:
7840:
7836:
7832:
7829:
7822:
7818:
7814:
7807:
7803:
7799:
7791:
7788:
7785:
7782:
7775:
7771:
7767:
7760:
7756:
7752:
7749:
7742:
7738:
7734:
7727:
7723:
7719:
7712:
7709:
7706:
7699:
7695:
7691:
7688:
7681:
7677:
7673:
7666:
7662:
7658:
7655:
7648:
7644:
7640:
7633:
7630:
7627:
7619:
7615:
7611:
7608:
7601:
7597:
7593:
7586:
7582:
7578:
7565:
7549:
7545:
7522:
7518:
7514:
7511:
7508:
7503:
7499:
7476:
7472:
7449:
7445:
7422:
7418:
7414:
7411:
7408:
7403:
7399:
7376:
7372:
7349:
7322:
7308:
7292:
7286:
7279:
7276:
7270:
7265:
7255:
7252:
7242:
7206:
7189:
7186:
7183:
7181:
7171:
7154:
7149:
7143:
7140:
7133:
7127:
7122:
7117:
7111:
7108:
7101:
7095:
7090:
7085:
7079:
7072:
7066:
7059:
7053:
7048:
7041:
7038:
7033:
7027:
7022:
7020:
7010:
6993:
6988:
6982:
6975:
6969:
6964:
6959:
6953:
6950:
6943:
6937:
6932:
6927:
6921:
6914:
6908:
6901:
6895:
6890:
6883:
6880:
6875:
6869:
6864:
6862:
6852:
6829:
6814:
6808:
6801:
6795:
6790:
6782:
6747:
6728:
6659:
6642:
6637:
6631:
6628:
6621:
6615:
6610:
6605:
6599:
6592:
6586:
6579:
6573:
6568:
6565:
6558:
6553:
6547:
6542:
6540:
6526:
6521:
6515:
6510:
6503:
6500:
6495:
6489:
6484:
6479:
6473:
6468:
6461:
6456:
6450:
6443:
6437:
6432:
6429:
6422:
6417:
6411:
6406:
6403:
6401:
6381:
6368:
6363:
6357:
6352:
6349:
6342:
6337:
6331:
6326:
6321:
6318:
6312:
6306:
6301:
6294:
6289:
6283:
6277:
6272:
6269:
6231:
6218:
6213:
6207:
6202:
6195:
6190:
6184:
6179:
6162:
6156:
6151:
6144:
6139:
6133:
6128:
6034:
6024:
6021:
6011:
5989:
5986:
5976:
5973:
5961:
5960:
5959:
5940:
5932:
5914:
5906:
5903:
5900:
5862:
5844:
5831:
5828:
5818:
5810:
5807:
5804:
5785:
5772:
5767:
5761:
5754:
5748:
5743:
5726:
5720:
5715:
5708:
5703:
5697:
5692:
5679:
5678:is given by:
5660:
5636:
5629:
5613:
5607:
5604:
5597:
5591:
5586:
5576:
5573:
5563:
5549:
5532:
5527:
5521:
5518:
5511:
5505:
5500:
5495:
5489:
5486:
5479:
5473:
5468:
5463:
5457:
5454:
5447:
5441:
5434:
5428:
5425:
5420:
5413:
5410:
5405:
5399:
5394:
5392:
5382:
5365:
5360:
5354:
5351:
5344:
5338:
5333:
5328:
5322:
5319:
5312:
5306:
5301:
5296:
5290:
5287:
5280:
5274:
5267:
5261:
5258:
5253:
5246:
5243:
5238:
5232:
5227:
5225:
5215:
5198:
5193:
5187:
5184:
5177:
5171:
5166:
5161:
5155:
5152:
5145:
5139:
5134:
5129:
5123:
5120:
5113:
5107:
5100:
5094:
5091:
5086:
5079:
5076:
5071:
5065:
5060:
5058:
5048:
5031:
5026:
5020:
5017:
5010:
5004:
4999:
4994:
4988:
4985:
4978:
4972:
4967:
4962:
4956:
4953:
4946:
4940:
4933:
4927:
4924:
4919:
4912:
4909:
4904:
4898:
4893:
4891:
4881:
4864:
4859:
4853:
4850:
4843:
4837:
4832:
4827:
4821:
4818:
4811:
4805:
4800:
4795:
4789:
4786:
4779:
4773:
4766:
4760:
4757:
4752:
4745:
4742:
4737:
4731:
4726:
4724:
4714:
4697:
4692:
4686:
4683:
4676:
4670:
4665:
4660:
4654:
4651:
4644:
4638:
4633:
4628:
4622:
4619:
4612:
4606:
4599:
4593:
4590:
4585:
4578:
4575:
4570:
4564:
4559:
4557:
4547:
4530:
4525:
4519:
4516:
4509:
4503:
4498:
4493:
4487:
4484:
4477:
4471:
4466:
4461:
4455:
4448:
4442:
4435:
4429:
4426:
4421:
4414:
4411:
4406:
4400:
4395:
4393:
4383:
4360:
4346:
4341:
4335:
4328:
4322:
4317:
4309:
4274:
4255:
4186:
4169:
4164:
4158:
4155:
4148:
4142:
4137:
4132:
4126:
4119:
4113:
4106:
4100:
4097:
4092:
4085:
4080:
4074:
4069:
4067:
4053:
4048:
4042:
4039:
4034:
4027:
4024:
4019:
4013:
4008:
4003:
3997:
3992:
3985:
3980:
3974:
3967:
3961:
3958:
3953:
3946:
3941:
3935:
3930:
3927:
3925:
3905:
3892:
3887:
3881:
3878:
3873:
3866:
3861:
3855:
3850:
3845:
3842:
3836:
3830:
3827:
3822:
3815:
3810:
3804:
3798:
3793:
3790:
3752:
3739:
3734:
3728:
3723:
3716:
3711:
3705:
3700:
3683:
3677:
3674:
3669:
3662:
3657:
3651:
3646:
3621:
3554:
3544:
3541:
3531:
3509:
3506:
3496:
3493:
3481:
3480:
3479:
3460:
3452:
3434:
3426:
3423:
3420:
3382:
3364:
3351:
3348:
3338:
3330:
3327:
3324:
3305:
3292:
3287:
3281:
3274:
3268:
3263:
3246:
3240:
3237:
3232:
3225:
3220:
3214:
3209:
3197:is given by:
3179:
3148:
3140:
3135:
3131:
3127:
3122:
3118:
3113:
3109:
3100:
3097:
3092:
3088:
3084:
3080:
3076:
3072:
3068:
3058:
3055:
3049:
3045:
3039:
3034:
3030:
3024:
3020:
3014:
3010:
3005:
2999:
2996:
2985:
2982:
2978:
2974:
2966:
2952:
2949:
2946:
2899:
2862:
2859:
2849:
2819:
2811:
2793:
2785:
2782:
2779:
2739:
2702:
2699:
2684:
2658:
2655:
2642:
2639:
2636:
2611:
2603:
2567:
2561:
2536:
2533:
2530:is symmetric
2508:
2507:
2506:
2476:
2462:
2459:
2456:
2448:
2443:
2441:
2437:
2436:critical path
2433:
2428:
2411:
2408:
2405:
2377:
2342:
2339:
2336:
2301:
2298:
2295:
2276:
2274:
2273:Jacobi method
2264:
2262:
2244:
2240:
2217:
2213:
2207:
2204:
2200:
2194:
2191:
2188:
2184:
2174:
2161:
2158:
2155:
2152:
2149:
2146:
2143:
2140:
2137:
2134:
2130:
2126:
2117:
2109:
2105:
2099:
2096:
2092:
2086:
2081:
2078:
2075:
2072:
2069:
2065:
2061:
2053:
2050:
2047:
2039:
2035:
2029:
2026:
2022:
2016:
2013:
2010:
2005:
2002:
1999:
1995:
1991:
1986:
1982:
1977:
1969:
1966:
1962:
1958:
1953:
1945:
1942:
1939:
1931:
1927:
1918:
1902:
1877:
1874:
1871:
1825:
1812:
1808:
1799:
1781:
1772:
1766:
1763:
1753:
1745:
1742:
1739:
1721:
1697:
1644:
1636:
1634:
1610:
1608:
1593:
1571:
1569:
1551:
1531:
1529:
1506:
1505:
1504:
1496:
1483:
1468:
1464:
1458:
1453:
1448:
1443:
1436:
1431:
1426:
1421:
1412:
1409:
1405:
1399:
1394:
1389:
1380:
1377:
1373:
1367:
1360:
1356:
1350:
1344:
1337:
1322:
1318:
1310:
1307:
1303:
1297:
1290:
1287:
1283:
1275:
1272:
1268:
1260:
1255:
1250:
1245:
1238:
1233:
1226:
1222:
1214:
1210:
1202:
1197:
1192:
1185:
1181:
1174:
1167:
1154:
1153:is given by:
1140:
1118:
1114:
1093:
1068:
1060:
1025:
999:
961:
953:
935:
927:
919:
916:
913:
884:
870:
867:
864:
817:
814:
811:
759:
734:
704:
701:
698:
695:
692:
689:
686:
683:
680:
677:
674:
654:
651:
643:
635:
581:
564:. The vector
474:
461:
456:
448:
444:
436:
427:
423:
413:
409:
402:
397:
388:
383:
375:
371:
363:
354:
350:
340:
336:
329:
324:
315:
310:
302:
299:
295:
289:
282:
279:
275:
267:
264:
260:
252:
247:
242:
237:
228:
225:
221:
215:
208:
204:
196:
192:
182:
179:
175:
169:
162:
158:
150:
146:
139:
134:
121:
98:
74:
72:
68:
64:
60:
56:
52:
49:
46:
42:
38:
34:
30:
26:
22:
10458:Key concepts
10389:Gauss–Seidel
10370:
10355:
10340:
10330:
10306:
10289:
10271:, Thm 10.1.2
10264:
10232:(1): 93–97.
10229:
10225:
10219:
10207:
10195:
10183:
10177:(3): 81–108.
10174:
10170:
10160:
10144:
10125:
10119:
9879:gauss_seidel
9598:
9535:
8895:
8883:
8706:
8333:(0, 0, 0, 0)
8330:
7869:Solving for
7868:
7566:
7464:in terms of
7391:in terms of
7314:
7207:
6830:
6729:
6705:the vectors
6660:
6382:
6232:
6049:
5883:in the form
5786:
5680:
5642:
5630:
5550:
4361:
4256:
4232:the vectors
4187:
3906:
3753:
3569:
3403:in the form
3306:
3161:
3146:
3138:
3133:
3129:
3125:
3120:
3116:
3111:
3107:
3098:
3095:
3090:
3086:
3082:
3078:
3074:
3070:
3066:
3056:
3053:
3047:
3043:
3037:
3032:
3028:
3022:
3018:
3012:
3008:
3003:
2997:
2994:
2983:
2980:
2976:
2972:
2837:converge to
2568:
2565:
2482:
2444:
2429:
2277:
2270:
2175:
1831:
1699:
1675:
1673:
1502:
1155:
890:
717:. Denote by
475:
122:
80:
32:
28:
24:
18:
10410:Matlab code
10226:SIAM Review
10153:direct link
2757:defined by
2479:Convergence
77:Description
10603:Categories
10500:algorithms
10279:References
10149:Gauss 1903
10094:this paper
9586:Solution:
9204:zeros_like
9141:zeros_like
6383:Now find:
6050:Decompose
3907:Now find:
3570:Decompose
2267:Discussion
1048:such that
1026:component
1000:component
10527:CPU cache
10377:EMS Press
10353:license.
10332:MathWorld
10234:CiteSeerX
9852:…
9815:…
9735:∑
9731:−
9676:∑
9672:−
8864:0.998351
8850:0.984341
8836:0.878864
8833:−0.987273
8655:−
8626:−
8596:−
8587:−
8567:−
8550:−
8502:−
8257:−
8219:−
8156:−
8094:−
7996:−
7815:−
7786:−
7768:−
7735:−
7674:−
7641:−
7594:−
7512:…
7412:…
7277:−
7253:−
7187:⋯
7141:−
7109:−
7039:−
6951:−
6881:−
6629:−
6566:−
6501:−
6430:−
6407:−
6350:−
6319:−
6270:−
6022:−
5987:−
5977:−
5845:−
5829:−
5605:−
5587:≈
5574:−
5519:−
5487:−
5455:−
5426:−
5411:−
5352:−
5320:−
5288:−
5259:−
5244:−
5185:−
5153:−
5121:−
5092:−
5077:−
5018:−
4986:−
4954:−
4925:−
4910:−
4851:−
4819:−
4787:−
4758:−
4743:−
4684:−
4652:−
4620:−
4591:−
4576:−
4517:−
4485:−
4427:−
4412:−
4156:−
4098:−
4040:−
4025:−
3959:−
3931:−
3879:−
3843:−
3828:−
3791:−
3675:−
3542:−
3507:−
3497:−
3365:−
3349:−
3238:−
3149:(repeat)
2977:algorithm
2969:Algorithm
2860:−
2700:−
2656:−
2643:ρ
2612:−
2457:ω
2192:≠
2185:∑
2153:…
2066:∑
2062:−
2014:−
1996:∑
1992:−
1782:−
1764:−
1645:−
1469:⏟
1454:⋯
1437:⋮
1432:⋱
1427:⋮
1422:⋮
1400:⋯
1368:⋯
1323:⏟
1298:⋯
1261:⋮
1256:⋱
1251:⋮
1246:⋮
1234:⋯
1198:⋯
1119:∗
936:−
794:, and by
437:⋮
364:⋮
290:⋯
253:⋮
248:⋱
243:⋮
238:⋮
216:⋯
170:⋯
63:symmetric
35:, is an
10556:Software
10520:Hardware
10479:Problems
10347:CFD-Wiki
10305:(1996),
10288:(1903),
10069:See also
9872:function
9402:allclose
9249:it_count
9168:it_count
8878:0.99985
8875:−1.00031
8861:−1.00253
8847:−1.01446
8331:Suppose
3153:Examples
2432:parallel
2261:residual
1022:, and a
621:, often
10405:Bickson
10379:, 2001
10256:2132758
9991:'.*
9589:Error:
8869:1.00086
8858:2.00356
8855:1.00659
8844:2.03694
8841:1.03018
8830:2.32727
8689:0.8789.
8584:0.23272
7949:gives:
5958:where:
3478:where:
2995:output:
2984:inputs:
2681:be the
71:Gerling
31:or the
10578:LAPACK
10568:MATLAB
10313:
10254:
10236:
10132:
9528:"
9465:"
9267:"
9228:float_
9159:float_
9099:format
9057:"
8901:import
8872:2.0003
8658:0.9873
8635:2.3272
8599:0.9873
8536:2.3272
8475:2.3272
5608:0.6650
5598:0.8122
5522:0.6650
5512:0.8122
5490:0.7443
5480:0.6875
5458:0.6650
5448:0.8122
5429:0.1193
5414:0.1875
5355:0.6650
5345:0.8122
5323:0.7443
5313:0.6875
5291:0.6650
5281:0.8121
5262:0.1193
5247:0.1875
5188:0.6650
5178:0.8121
5156:0.7443
5146:0.6875
5124:0.6646
5114:0.8127
5095:0.1193
5080:0.1875
5021:0.6646
5011:0.8127
4989:0.7443
4979:0.6875
4957:0.6678
4947:0.8077
4928:0.1193
4913:0.1875
4854:0.6678
4844:0.8077
4822:0.7443
4812:0.6875
4790:0.6413
4780:0.8494
4761:0.1193
4746:0.1875
4687:0.6413
4677:0.8494
4655:0.7443
4645:0.6875
4623:0.8636
4613:0.5000
4594:0.1193
4579:0.1875
4520:0.8636
4510:0.5000
4488:0.7443
4478:0.6875
4430:0.1193
4415:0.1875
4159:0.7439
4149:0.6875
4101:0.0909
4093:0.0398
4086:0.0000
4081:0.0625
4043:0.1194
4028:0.1875
3962:0.0909
3954:0.0398
3947:0.0000
3942:0.0625
3882:0.0909
3874:0.0398
3867:0.0000
3862:0.0625
3096:end if
3009:repeat
1915:using
45:German
23:, the
10563:ATLAS
10290:Werke
10252:JSTOR
10111:Notes
9904:iters
9522:error
9507:print
9471:error
9444:print
9441:x_new
9432:break
9414:x_new
9360:x_new
9324:x_new
9294:shape
9282:range
9234:print
9216:dtype
9192:x_new
9174:range
9087:print
9081:shape
9069:range
9018:shape
9006:range
8985:print
8979:array
8937:array
8904:numpy
7144:1.651
7134:4.911
7112:2.071
7102:5.500
7080:0.393
7073:2.050
7054:1.071
7042:1.500
6983:0.393
6976:2.050
6954:2.071
6944:5.500
6896:1.071
6884:1.500
6661:With
6632:2.071
6622:5.500
6574:0.143
6569:0.357
6559:0.000
6554:0.500
6516:1.071
6511:0.000
6504:1.500
6496:0.000
6438:0.143
6433:0.357
6423:0.000
6418:0.500
6358:0.143
6353:0.357
6343:0.000
6338:0.500
5421:0.000
5406:0.000
5254:0.000
5239:0.000
5087:0.000
5072:0.000
4920:0.000
4905:0.000
4753:0.000
4738:0.000
4586:0.000
4571:0.000
4422:0.000
4407:0.000
4188:With
4035:0.000
4020:0.000
3048:until
3023:until
2449:with
1106:into
476:When
61:, or
10542:SIMD
10351:GFDL
10311:ISBN
10130:ISBN
9922:size
9742:>
9683:<
9426:1e-8
9420:rtol
9111:join
8919:1000
7922:and
6683:and
6255:is:
4210:and
3776:is:
3128:) /
3067:then
3044:from
3019:from
2950:<
2534:, or
1133:and
752:the
667:for
498:and
81:Let
65:and
53:and
10532:TLB
10345:on
10244:doi
10063:end
10060:end
10057:end
9988:,:)
9973:sum
9907:for
9889:for
9483:dot
9342:dot
9312:dot
9273:for
9165:for
9117:row
9060:for
9024:row
8997:for
8827:0.6
8590:1.1
8378:0.6
7851:15.
6922:2.3
6915:1.1
6809:2.3
6802:1.1
6170:and
6002:and
5734:and
4456:1.0
4449:1.0
4336:1.0
4329:1.0
3691:and
3522:and
3254:and
3147:end
3139:end
3115:← (
3099:end
3038:for
3035:← 0
3013:for
2917:if
2685:of
19:In
10605::
10375:,
10369:,
10329:.
10301:;
10250:.
10242:.
10230:37
10228:.
10175:11
10096:.)
10054:);
10033:))
9553:=
9550:=
9547:=
9544:=
9477:np
9429:):
9396:np
9393:if
9381:s2
9375:s1
9336:np
9330:s2
9306:np
9300:s1
9297:):
9279:in
9255::
9222:np
9198:np
9189:):
9171:in
9153:np
9135:np
9126:))
9120:),
9084:)]
9066:in
9039:*x
9033:3g
9021:):
9003:in
8982:()
8973:np
8931:np
8910:np
8907:as
8888:.
8675:15
8578:25
8561:10
8553:11
8547:10
8469:11
8461:25
8455:55
8441:11
8433:25
8427:11
8313:8.
8305:15
8230:10
8222:11
8216:10
8195:10
8129:11
8121:25
8115:11
8091:11
8070:11
7993:10
7789:11
7753:10
7710:25
7659:11
7579:10
7287:29
7280:38
6600:13
6593:11
6116::
5762:13
5755:11
5628:.
4127:13
4120:11
3831:11
3811:16
3678:11
3658:16
3634::
3282:13
3275:11
3241:11
3221:16
3134:ii
3124:−
3083:ij
3077:+
3073:←
3062:≠
3057:if
3054:do
3046:1
3029:do
3021:1
2989:,
2981:is
2965:.
2475:.
2263:.
1919::
1361:12
1227:22
1215:21
1186:11
209:22
197:21
163:12
151:11
10502:)
10498:(
10443:e
10436:t
10429:v
10335:.
10320:.
10294:.
10258:.
10246::
10155:.
10138:.
10051:j
10048:,
10045:j
10042:(
10039:A
10036:/
10030:j
10027:(
10024:x
10021:*
10018:)
10015:j
10012:,
10009:j
10006:(
10003:A
10000:+
9997:)
9994:x
9985:j
9982:(
9979:A
9976:(
9970:-
9967:)
9964:j
9961:(
9958:b
9955:(
9952:=
9949:)
9946:j
9943:(
9940:x
9937:)
9934:1
9931:,
9928:A
9925:(
9919::
9916:1
9913:=
9910:j
9901::
9898:1
9895:=
9892:i
9886:)
9882:(
9876:=
9874:x
9849:,
9846:2
9843:,
9840:1
9837:,
9834:0
9831:=
9828:k
9821:n
9818:,
9812:,
9809:2
9806:,
9803:1
9800:=
9797:i
9789:,
9785:)
9779:)
9776:k
9773:(
9768:j
9764:x
9758:j
9755:i
9751:a
9745:i
9739:j
9726:)
9723:1
9720:+
9717:k
9714:(
9709:j
9705:x
9699:j
9696:i
9692:a
9686:i
9680:j
9667:i
9663:b
9658:(
9650:i
9647:i
9643:a
9639:1
9634:=
9629:)
9626:1
9623:+
9620:k
9617:(
9612:i
9608:x
9531:)
9525:}
9519:{
9513:f
9510:(
9504:b
9501:-
9498:)
9495:x
9492:,
9489:A
9486:(
9480:.
9474:=
9468:)
9462:}
9459:x
9456:{
9450:f
9447:(
9438:=
9435:x
9423:=
9417:,
9411:,
9408:x
9405:(
9399:.
9390:A
9387:/
9384:)
9378:-
9372:-
9369:b
9366:(
9363:=
9357:)
9354:x
9351:,
9348:A
9345:(
9339:.
9333:=
9327:)
9321:,
9318:A
9315:(
9309:.
9303:=
9291:.
9288:A
9285:(
9276:i
9270:)
9264:}
9261:x
9258:{
9252:}
9246:{
9240:f
9237:(
9231:)
9225:.
9219:=
9213:,
9210:x
9207:(
9201:.
9195:=
9183:,
9180:1
9177:(
9162:)
9156:.
9150:,
9147:b
9144:(
9138:.
9132:=
9129:x
9123:b
9114:(
9108:.
9102:(
9096:.
9090:(
9078:.
9075:A
9072:(
9063:j
9054:}
9051:1
9048:+
9045:j
9042:{
9036:}
9030::
9027:=
9015:.
9012:A
9009:(
9000:i
8994:)
8988:(
8976:.
8970:=
8967:b
8961:)
8958:]
8955:,
8952:,
8949:,
8946:,
8943:[
8940:(
8934:.
8928:=
8925:A
8916:=
8809:4
8805:x
8781:3
8777:x
8753:2
8749:x
8725:1
8721:x
8686:=
8683:8
8679:/
8672:+
8669:8
8665:/
8661:)
8652:(
8649:+
8646:8
8642:/
8638:)
8632:(
8629:3
8623:=
8614:4
8610:x
8602:,
8593:=
8581:+
8574:/
8570:3
8564:=
8557:/
8543:/
8539:)
8533:(
8530:+
8527:5
8523:/
8519:)
8516:5
8512:/
8508:3
8505:(
8499:=
8490:3
8486:x
8478:,
8472:=
8465:/
8458:+
8451:/
8447:3
8444:=
8437:/
8430:+
8423:/
8419:)
8416:5
8412:/
8408:3
8405:(
8402:=
8393:2
8389:x
8381:,
8375:=
8372:5
8368:/
8364:3
8361:=
8352:1
8348:x
8309:/
8302:+
8299:8
8295:/
8289:3
8285:x
8281:+
8278:8
8274:/
8268:2
8264:x
8260:3
8254:=
8245:4
8241:x
8233:,
8226:/
8212:/
8206:4
8202:x
8198:+
8191:/
8185:2
8181:x
8177:+
8174:5
8170:/
8164:1
8160:x
8153:=
8144:3
8140:x
8132:,
8125:/
8118:+
8111:/
8105:4
8101:x
8097:3
8087:/
8081:3
8077:x
8073:+
8066:/
8060:1
8056:x
8052:=
8043:2
8039:x
8031:,
8028:5
8024:/
8020:3
8017:+
8014:5
8010:/
8004:3
8000:x
7989:/
7983:2
7979:x
7975:=
7966:1
7962:x
7935:4
7931:x
7908:3
7904:x
7900:,
7895:2
7891:x
7887:,
7882:1
7878:x
7848:=
7841:4
7837:x
7833:8
7830:+
7823:3
7819:x
7808:2
7804:x
7800:3
7792:,
7783:=
7776:4
7772:x
7761:3
7757:x
7750:+
7743:2
7739:x
7728:1
7724:x
7720:2
7713:,
7707:=
7700:4
7696:x
7692:3
7689:+
7682:3
7678:x
7667:2
7663:x
7656:+
7649:1
7645:x
7634:,
7631:6
7628:=
7620:3
7616:x
7612:2
7609:+
7602:2
7598:x
7587:1
7583:x
7550:n
7546:x
7523:n
7519:x
7515:,
7509:,
7504:3
7500:x
7477:1
7473:x
7450:2
7446:x
7423:n
7419:x
7415:,
7409:,
7404:2
7400:x
7377:1
7373:x
7350:0
7345:x
7323:n
7293:]
7271:[
7266:=
7262:b
7256:1
7248:A
7243:=
7239:x
7217:A
7190:.
7184:=
7175:)
7172:3
7169:(
7164:x
7155:.
7150:]
7128:[
7123:=
7118:]
7096:[
7091:+
7086:]
7067:[
7060:]
7049:0
7034:0
7028:[
7023:=
7014:)
7011:2
7008:(
7003:x
6994:.
6989:]
6970:[
6965:=
6960:]
6938:[
6933:+
6928:]
6909:[
6902:]
6891:0
6876:0
6870:[
6865:=
6856:)
6853:1
6850:(
6845:x
6815:]
6796:[
6791:=
6786:)
6783:0
6780:(
6775:x
6751:)
6748:0
6745:(
6740:x
6714:x
6692:c
6670:T
6643:.
6638:]
6616:[
6611:=
6606:]
6587:[
6580:]
6548:[
6543:=
6535:c
6527:,
6522:]
6490:[
6485:=
6480:]
6474:0
6469:0
6462:3
6457:0
6451:[
6444:]
6412:[
6404:=
6396:T
6369:.
6364:]
6332:[
6327:=
6322:1
6313:]
6307:7
6302:5
6295:0
6290:2
6284:[
6278:=
6273:1
6265:L
6242:L
6219:.
6214:]
6208:0
6203:0
6196:3
6191:0
6185:[
6180:=
6176:U
6163:]
6157:7
6152:5
6145:0
6140:2
6134:[
6129:=
6125:L
6103:U
6081:L
6059:A
6035:.
6031:b
6025:1
6017:L
6012:=
6008:c
5996:U
5990:1
5982:L
5974:=
5970:T
5945:c
5941:+
5936:)
5933:k
5930:(
5925:x
5919:T
5915:=
5910:)
5907:1
5904:+
5901:k
5898:(
5893:x
5871:)
5866:)
5863:k
5860:(
5855:x
5849:U
5841:b
5837:(
5832:1
5824:L
5819:=
5814:)
5811:1
5808:+
5805:k
5802:(
5797:x
5773:.
5768:]
5749:[
5744:=
5740:b
5727:]
5721:7
5716:5
5709:3
5704:2
5698:[
5693:=
5689:A
5665:b
5661:=
5657:x
5652:A
5633:A
5614:]
5592:[
5583:b
5577:1
5569:A
5564:=
5560:x
5533:.
5528:]
5506:[
5501:=
5496:]
5474:[
5469:+
5464:]
5442:[
5435:]
5400:[
5395:=
5386:)
5383:7
5380:(
5375:x
5366:.
5361:]
5339:[
5334:=
5329:]
5307:[
5302:+
5297:]
5275:[
5268:]
5233:[
5228:=
5219:)
5216:6
5213:(
5208:x
5199:.
5194:]
5172:[
5167:=
5162:]
5140:[
5135:+
5130:]
5108:[
5101:]
5066:[
5061:=
5052:)
5049:5
5046:(
5041:x
5032:.
5027:]
5005:[
5000:=
4995:]
4973:[
4968:+
4963:]
4941:[
4934:]
4899:[
4894:=
4885:)
4882:4
4879:(
4874:x
4865:.
4860:]
4838:[
4833:=
4828:]
4806:[
4801:+
4796:]
4774:[
4767:]
4732:[
4727:=
4718:)
4715:3
4712:(
4707:x
4698:.
4693:]
4671:[
4666:=
4661:]
4639:[
4634:+
4629:]
4607:[
4600:]
4565:[
4560:=
4551:)
4548:2
4545:(
4540:x
4531:.
4526:]
4504:[
4499:=
4494:]
4472:[
4467:+
4462:]
4443:[
4436:]
4401:[
4396:=
4387:)
4384:1
4381:(
4376:x
4347:.
4342:]
4323:[
4318:=
4313:)
4310:0
4307:(
4302:x
4278:)
4275:0
4272:(
4267:x
4241:x
4219:c
4197:T
4170:.
4165:]
4143:[
4138:=
4133:]
4114:[
4107:]
4075:[
4070:=
4062:c
4054:,
4049:]
4014:[
4009:=
4004:]
3998:0
3993:0
3986:3
3981:0
3975:[
3968:]
3936:[
3928:=
3920:T
3893:.
3888:]
3856:[
3851:=
3846:1
3837:]
3823:7
3816:0
3805:[
3799:=
3794:1
3786:L
3763:L
3740:.
3735:]
3729:0
3724:0
3717:3
3712:0
3706:[
3701:=
3697:U
3684:]
3670:7
3663:0
3652:[
3647:=
3643:L
3622:U
3601:L
3579:A
3555:.
3551:b
3545:1
3537:L
3532:=
3528:c
3516:U
3510:1
3502:L
3494:=
3490:T
3465:c
3461:+
3456:)
3453:k
3450:(
3445:x
3439:T
3435:=
3430:)
3427:1
3424:+
3421:k
3418:(
3413:x
3391:)
3386:)
3383:k
3380:(
3375:x
3369:U
3361:b
3357:(
3352:1
3344:L
3339:=
3334:)
3331:1
3328:+
3325:k
3322:(
3317:x
3293:.
3288:]
3269:[
3264:=
3260:b
3247:]
3233:7
3226:3
3215:[
3210:=
3206:A
3184:b
3180:=
3176:x
3171:A
3143:i
3141:(
3130:a
3126:σ
3121:i
3117:b
3112:i
3108:φ
3103:j
3101:(
3091:j
3087:φ
3079:a
3075:σ
3071:σ
3064:i
3060:j
3051:n
3041:j
3033:σ
3026:n
3016:i
3004:φ
2998:φ
2991:b
2987:A
2953:1
2947:r
2926:M
2903:)
2900:0
2897:(
2892:x
2869:b
2863:1
2855:A
2850:=
2846:x
2824:b
2820:+
2815:)
2812:k
2809:(
2804:x
2798:N
2794:=
2789:)
2786:1
2783:+
2780:k
2777:(
2772:x
2766:M
2743:)
2740:k
2737:(
2732:x
2709:N
2703:1
2695:M
2669:)
2665:N
2659:1
2651:M
2646:(
2640:=
2637:r
2616:N
2608:M
2604:=
2600:A
2578:A
2562:.
2545:A
2517:A
2492:A
2463:1
2460:=
2415:)
2412:1
2409:+
2406:k
2403:(
2381:)
2378:k
2375:(
2370:x
2346:)
2343:1
2340:+
2337:k
2334:(
2329:x
2305:)
2302:1
2299:+
2296:k
2293:(
2288:x
2245:j
2241:x
2218:j
2214:x
2208:j
2205:i
2201:a
2195:i
2189:j
2162:.
2159:n
2156:,
2150:,
2147:2
2144:,
2141:1
2138:=
2135:i
2131:,
2127:)
2121:)
2118:k
2115:(
2110:j
2106:x
2100:j
2097:i
2093:a
2087:n
2082:1
2079:+
2076:i
2073:=
2070:j
2057:)
2054:1
2051:+
2048:k
2045:(
2040:j
2036:x
2030:j
2027:i
2023:a
2017:1
2011:i
2006:1
2003:=
2000:j
1987:i
1983:b
1978:(
1970:i
1967:i
1963:a
1959:1
1954:=
1949:)
1946:1
1943:+
1940:k
1937:(
1932:i
1928:x
1903:i
1881:)
1878:1
1875:+
1872:k
1869:(
1864:x
1841:L
1813:.
1809:)
1803:)
1800:k
1797:(
1792:x
1786:U
1778:b
1773:(
1767:1
1759:L
1754:=
1749:)
1746:1
1743:+
1740:k
1737:(
1732:x
1708:x
1684:x
1654:x
1649:U
1641:b
1637:=
1629:x
1624:L
1615:b
1611:=
1603:x
1598:U
1594:+
1590:x
1585:L
1576:b
1572:=
1564:x
1560:)
1556:U
1552:+
1548:L
1544:(
1536:b
1532:=
1524:x
1519:A
1484:.
1477:U
1465:]
1459:0
1449:0
1444:0
1413:n
1410:2
1406:a
1395:0
1390:0
1381:n
1378:1
1374:a
1357:a
1351:0
1345:[
1338:+
1331:L
1319:]
1311:n
1308:n
1304:a
1291:2
1288:n
1284:a
1276:1
1273:n
1269:a
1239:0
1223:a
1211:a
1203:0
1193:0
1182:a
1175:[
1168:=
1164:A
1141:U
1115:L
1094:A
1073:U
1069:+
1065:L
1061:=
1057:A
1035:U
1009:L
983:A
962:,
957:)
954:k
951:(
946:x
940:U
932:b
928:=
923:)
920:1
917:+
914:k
911:(
906:x
900:L
871:1
868:+
865:k
844:x
821:)
818:1
815:+
812:k
809:(
804:x
781:x
760:k
738:)
735:k
732:(
727:x
705:n
702:,
699:.
696:.
693:.
690:,
687:2
684:,
681:1
678:=
675:i
655:0
652:=
647:)
644:0
641:(
636:i
631:x
608:x
585:)
582:0
579:(
574:x
551:x
529:x
507:b
485:A
462:.
457:]
449:n
445:b
428:2
424:b
414:1
410:b
403:[
398:=
394:b
389:,
384:]
376:n
372:x
355:2
351:x
341:1
337:x
330:[
325:=
321:x
316:,
311:]
303:n
300:n
296:a
283:2
280:n
276:a
268:1
265:n
261:a
229:n
226:2
222:a
205:a
193:a
183:n
180:1
176:a
159:a
147:a
140:[
135:=
131:A
118:n
103:b
99:=
95:x
90:A
Text is available under the Creative Commons Attribution-ShareAlike License. Additional terms may apply.