3485:
22:
4433:
4301:
4169:
3949:
3845:
1500:
4538:
1983:
1000:
894:
812:
728:
3679:
1052:
948:
676:
594:
435:
383:
329:
247:
4032:
4307:
4175:
4043:
1368:
1655:
124:
3736:
1731:
1563:
1176:
1116:
1082:
624:
3549:
1889:
1238:
1146:
842:
760:
529:
499:
465:
277:
1817:
1379:
1280:
3479:
2017:
3603:
3576:
3449:
2076:
1844:
1778:
1751:
1682:
1590:
2041:
3851:
3747:
3497:
At the cost of adding an additional implicit relation, these methods can be adapted to have time reversal symmetry. In these methods, the averaged position
32:
4615:
90:
62:
2082:
to converge arbitrarily close to the true solution. Below is a Matlab function which implements the Gauss-Legendre method of order four.
4439:
69:
4634:
1894:
954:
848:
766:
682:
3611:
1006:
902:
630:
548:
389:
337:
283:
201:
76:
3481:
in as few as 2 Newton steps. The only extra work compared to explicit Runge-Kutta methods is the computation of the
Jacobian.
2051:
The Gauss-Legendre methods are implicit, so in general they cannot be applied exactly. Instead one makes an educated guess of
3955:
1194:
47:
136:
4428:{\displaystyle k_{3}=f\left({\frac {x_{f}+x_{i}}{2}}+{\frac {\sqrt {15}}{30}}hk_{1}+{\frac {\sqrt {15}}{15}}hk_{2}\right)}
4296:{\displaystyle k_{2}=f\left({\frac {x_{f}+x_{i}}{2}}+{\frac {\sqrt {15}}{24}}hk_{1}-{\frac {\sqrt {15}}{24}}hk_{3}\right)}
4164:{\displaystyle k_{1}=f\left({\frac {x_{f}+x_{i}}{2}}-{\frac {\sqrt {15}}{15}}hk_{2}-{\frac {\sqrt {15}}{30}}hk_{3}\right)}
58:
1285:
4601:
1185:
The computational cost of higher-order Gauss–Legendre methods is usually excessive, and thus, they are rarely used.
4607:
1985:
is invoked to accomplish this. This is an implicit constraint that must be solved by a root finding algorithm like
1595:
128:
83:
3484:
3685:
1687:
1507:
1152:
1092:
1058:
600:
3500:
1849:
1199:
1122:
818:
736:
505:
475:
441:
253:
116:
1783:
1371:
132:
112:
2079:
1986:
1243:
4611:
3454:
3605:
in standard Runge-Kutta methods. The method of order 2 is just an implicit midpoint method.
1992:
39:
3581:
3554:
3427:
2054:
1822:
1756:
1736:
1660:
1568:
162:
158:
3944:{\displaystyle k_{2}=f\left({\frac {x_{f}+x_{i}}{2}}+{\frac {\sqrt {3}}{6}}hk_{1}\right)}
3840:{\displaystyle k_{1}=f\left({\frac {x_{f}+x_{i}}{2}}-{\frac {\sqrt {3}}{6}}hk_{2}\right)}
2026:
4628:
4597:
2020:
1495:{\displaystyle x(h)=x(0)+{\frac {h}{2}}\sum _{i=1}^{\ell }w_{i}k_{i}+O(h^{2\ell }),}
151:
21:
1780:
is still impossible to evaluate. To maintain truncation error of order
4603:
A First Course in the
Numerical Analysis of Differential Equations
3483:
4533:{\displaystyle x_{f}=x_{i}+{\frac {h}{18}}(5k_{1}+8k_{2}+5k_{3})}
2319:% return a time derivative and a Jacobian of that time derivative
1978:{\textstyle k_{i}=f{\left(x_{0}+h\sum _{j}a_{ij}k_{j}\right)}}
15:
995:{\displaystyle {\tfrac {5}{36}}+{\tfrac {1}{30}}{\sqrt {15}}}
889:{\displaystyle {\tfrac {5}{36}}-{\tfrac {1}{24}}{\sqrt {15}}}
807:{\displaystyle {\tfrac {5}{36}}+{\tfrac {1}{24}}{\sqrt {15}}}
723:{\displaystyle {\tfrac {5}{36}}-{\tfrac {1}{30}}{\sqrt {15}}}
191:
The Gauss–Legendre method of order four has
Butcher tableau:
3674:{\displaystyle k_{1}=f\left({\frac {x_{f}+x_{i}}{2}}\right)}
1047:{\displaystyle {\tfrac {2}{9}}+{\tfrac {1}{15}}{\sqrt {15}}}
943:{\displaystyle {\tfrac {1}{2}}+{\tfrac {1}{10}}{\sqrt {15}}}
671:{\displaystyle {\tfrac {2}{9}}-{\tfrac {1}{15}}{\sqrt {15}}}
589:{\displaystyle {\tfrac {1}{2}}-{\tfrac {1}{10}}{\sqrt {15}}}
538:
The Gauss–Legendre method of order six has
Butcher tableau:
1282:. The distinguishing feature of GLRK is the estimation of
430:{\displaystyle {\tfrac {1}{4}}+{\tfrac {1}{6}}{\sqrt {3}}}
378:{\displaystyle {\tfrac {1}{2}}+{\tfrac {1}{6}}{\sqrt {3}}}
324:{\displaystyle {\tfrac {1}{4}}-{\tfrac {1}{6}}{\sqrt {3}}}
242:{\displaystyle {\tfrac {1}{2}}-{\tfrac {1}{6}}{\sqrt {3}}}
4027:{\displaystyle x_{f}=x_{i}+{\frac {h}{2}}(k_{1}+k_{2})}
43:
1897:
1598:
1288:
1157:
1127:
1097:
1063:
1026:
1011:
974:
959:
922:
907:
868:
853:
823:
786:
771:
741:
702:
687:
650:
635:
605:
568:
553:
510:
480:
446:
409:
394:
357:
342:
303:
288:
258:
221:
206:
4442:
4310:
4178:
4046:
3958:
3854:
3750:
3688:
3614:
3584:
3557:
3503:
3457:
3430:
2057:
2029:
1995:
1852:
1825:
1786:
1759:
1739:
1690:
1663:
1571:
1510:
1382:
1246:
1202:
1155:
1125:
1095:
1061:
1009:
957:
905:
851:
821:
769:
739:
685:
633:
603:
551:
508:
478:
444:
392:
340:
286:
256:
204:
125:
numerical methods for ordinary differential equations
2523:
x, dynamics, dt, threshold, damping, max_iterations
542:
195:
169:
4037:The method of order 6 with 3 stages is as follows.
3741:The method of order 4 with 2 stages is as follows.
3424:This algorithm is surprisingly cheap. The error in
3365:'Newton did not converge by %d iterations.'
1193:Gauss-Legendre Runge-Kutta (GLRK) methods solve an
4532:
4427:
4295:
4163:
4026:
3943:
3839:
3730:
3673:
3597:
3570:
3543:
3473:
3443:
2070:
2035:
2011:
1977:
1883:
1838:
1811:
1772:
1745:
1725:
1676:
1649:
1584:
1557:
1494:
1362:
1274:
1232:
1170:
1140:
1110:
1076:
1046:
994:
942:
888:
836:
806:
754:
722:
670:
618:
588:
523:
493:
459:
429:
377:
323:
271:
241:
1363:{\textstyle x(h)-x_{0}=\int _{0}^{h}dt\,f(x(t))}
3488:An integrated orbit near the Lorenz attractor.
157:The Gauss–Legendre method of order two is the
2957:'Newton iteration %d: error is %f.\n'
2600:% Use explicit Euler steps as initial guesses
8:
2591:'damping should be between 0 and 1.'
48:introducing citations to additional sources
1989:. The values of the Runge-Kutta parameters
1650:{\textstyle c_{i}={\frac {1}{2}}h(1+r_{i})}
4521:
4505:
4489:
4469:
4460:
4447:
4441:
4414:
4395:
4386:
4367:
4352:
4339:
4332:
4315:
4309:
4282:
4263:
4254:
4235:
4220:
4207:
4200:
4183:
4177:
4150:
4131:
4122:
4103:
4088:
4075:
4068:
4051:
4045:
4015:
4002:
3985:
3976:
3963:
3957:
3930:
3911:
3896:
3883:
3876:
3859:
3853:
3826:
3807:
3792:
3779:
3772:
3755:
3749:
3722:
3706:
3693:
3687:
3655:
3642:
3635:
3619:
3613:
3589:
3583:
3562:
3556:
3533:
3524:
3511:
3502:
3462:
3456:
3435:
3429:
2062:
2056:
2028:
2000:
1994:
1963:
1950:
1940:
1924:
1914:
1902:
1896:
1863:
1851:
1830:
1824:
1797:
1785:
1764:
1758:
1738:
1708:
1695:
1689:
1668:
1662:
1638:
1612:
1603:
1597:
1576:
1570:
1543:
1515:
1509:
1477:
1458:
1448:
1438:
1427:
1413:
1381:
1338:
1326:
1321:
1308:
1287:
1266:
1245:
1204:
1203:
1201:
1156:
1154:
1126:
1124:
1096:
1094:
1062:
1060:
1037:
1025:
1010:
1008:
985:
973:
958:
956:
933:
921:
906:
904:
879:
867:
852:
850:
822:
820:
797:
785:
770:
768:
740:
738:
713:
701:
686:
684:
661:
649:
634:
632:
604:
602:
579:
567:
552:
550:
509:
507:
479:
477:
445:
443:
420:
408:
393:
391:
368:
356:
341:
339:
314:
302:
287:
285:
257:
255:
232:
220:
205:
203:
1753:. A further approximation is needed, as
38:Relevant discussion may be found on the
4580:
4568:
4556:
4549:
1891:. The Runge-Kutta implicit definition
127:. Gauss–Legendre methods are implicit
3578:instead of just the initial position
1733:of the Legendre polynomial of degree
139:. The Gauss–Legendre method based on
7:
3731:{\displaystyle x_{f}=x_{i}+hk_{1}}
1726:{\displaystyle P_{\ell }(r_{i})=0}
1558:{\displaystyle k_{i}=f(x(hc_{i}))}
14:
1171:{\displaystyle {\tfrac {5}{18}}}
1111:{\displaystyle {\tfrac {5}{18}}}
1077:{\displaystyle {\tfrac {5}{36}}}
619:{\displaystyle {\tfrac {5}{36}}}
31:relies largely or entirely on a
20:
3544:{\displaystyle (x_{f}+x_{i})/2}
1884:{\displaystyle O(h^{2\ell -1})}
1233:{\displaystyle {\dot {x}}=f(x)}
1141:{\displaystyle {\tfrac {4}{9}}}
837:{\displaystyle {\tfrac {2}{9}}}
755:{\displaystyle {\tfrac {1}{2}}}
524:{\displaystyle {\tfrac {1}{2}}}
494:{\displaystyle {\tfrac {1}{2}}}
460:{\displaystyle {\tfrac {1}{4}}}
272:{\displaystyle {\tfrac {1}{4}}}
150:All Gauss–Legendre methods are
4527:
4479:
4021:
3995:
3530:
3504:
1878:
1856:
1806:
1790:
1714:
1701:
1644:
1625:
1552:
1549:
1533:
1527:
1486:
1470:
1407:
1401:
1392:
1386:
1357:
1354:
1348:
1342:
1298:
1292:
1256:
1250:
1227:
1221:
1195:ordinary differential equation
131:. More specifically, they are
1:
1812:{\displaystyle O(h^{2\ell })}
1592:are the quadrature weights,
1565:are the sampled velocities,
4651:
4608:Cambridge University Press
2294:'Lorenz Attractor'
1275:{\displaystyle x(0)=x_{0}}
2019:can be determined from a
137:Gauss–Legendre quadrature
3474:{\displaystyle 10^{-12}}
2084:
3493:Time-symmetric variants
1657:are the abscissas, and
135:based on the points of
59:"Gauss–Legendre method"
4534:
4429:
4297:
4165:
4028:
3945:
3841:
3732:
3675:
3599:
3572:
3545:
3489:
3475:
3445:
2072:
2037:
2013:
2012:{\displaystyle a_{ij}}
1979:
1885:
1840:
1813:
1774:
1747:
1727:
1678:
1651:
1586:
1559:
1496:
1443:
1364:
1276:
1234:
1172:
1142:
1112:
1078:
1048:
996:
944:
890:
838:
808:
756:
724:
672:
620:
590:
525:
495:
461:
431:
379:
325:
273:
243:
159:implicit midpoint rule
121:Gauss–Legendre methods
4535:
4430:
4298:
4166:
4029:
3946:
3842:
3733:
3676:
3600:
3598:{\displaystyle x_{i}}
3573:
3571:{\displaystyle k_{i}}
3551:is used in computing
3546:
3487:
3476:
3446:
3444:{\displaystyle k_{i}}
2073:
2071:{\displaystyle k_{i}}
2038:
2014:
1980:
1886:
1841:
1839:{\displaystyle k_{i}}
1814:
1775:
1773:{\displaystyle k_{i}}
1748:
1746:{\displaystyle \ell }
1728:
1679:
1677:{\displaystyle r_{i}}
1652:
1587:
1585:{\displaystyle w_{i}}
1560:
1497:
1423:
1365:
1277:
1235:
1173:
1143:
1113:
1079:
1049:
997:
945:
891:
839:
809:
757:
725:
673:
621:
591:
526:
496:
462:
432:
380:
326:
274:
244:
4440:
4308:
4176:
4044:
3956:
3852:
3748:
3686:
3612:
3582:
3555:
3501:
3455:
3428:
2055:
2027:
1993:
1895:
1850:
1823:
1784:
1757:
1737:
1688:
1661:
1596:
1569:
1508:
1380:
1286:
1244:
1200:
1153:
1123:
1093:
1059:
1007:
955:
903:
849:
819:
767:
737:
683:
631:
601:
549:
506:
476:
442:
390:
338:
284:
254:
202:
117:scientific computing
44:improve this article
4635:Runge–Kutta methods
1372:Gaussian quadrature
1331:
133:collocation methods
129:Runge–Kutta methods
4530:
4425:
4293:
4161:
4024:
3941:
3837:
3728:
3671:
3595:
3568:
3541:
3490:
3471:
3441:
2068:
2033:
2009:
1975:
1945:
1881:
1836:
1809:
1770:
1743:
1723:
1674:
1647:
1582:
1555:
1492:
1360:
1317:
1272:
1230:
1168:
1166:
1138:
1136:
1108:
1106:
1074:
1072:
1044:
1035:
1020:
992:
983:
968:
940:
931:
916:
886:
877:
862:
834:
832:
804:
795:
780:
752:
750:
720:
711:
696:
668:
659:
644:
616:
614:
586:
577:
562:
521:
519:
491:
489:
457:
455:
427:
418:
403:
375:
366:
351:
321:
312:
297:
269:
267:
239:
230:
215:
143:points has order 2
113:numerical analysis
4617:978-0-521-55655-2
4477:
4405:
4401:
4377:
4373:
4362:
4273:
4269:
4245:
4241:
4230:
4141:
4137:
4113:
4109:
4098:
3993:
3921:
3917:
3906:
3817:
3813:
3802:
3665:
2309:= lorenz_dynamics
2047:Practical example
2036:{\displaystyle h}
1936:
1620:
1421:
1212:
1181:
1180:
1165:
1135:
1105:
1071:
1042:
1034:
1019:
990:
982:
967:
938:
930:
915:
884:
876:
861:
831:
802:
794:
779:
749:
718:
710:
695:
666:
658:
643:
613:
584:
576:
561:
534:
533:
518:
488:
454:
425:
417:
402:
373:
365:
350:
319:
311:
296:
266:
237:
229:
214:
187:
186:
109:
108:
94:
4642:
4620:
4584:
4578:
4572:
4566:
4560:
4554:
4539:
4537:
4536:
4531:
4526:
4525:
4510:
4509:
4494:
4493:
4478:
4470:
4465:
4464:
4452:
4451:
4434:
4432:
4431:
4426:
4424:
4420:
4419:
4418:
4406:
4397:
4396:
4391:
4390:
4378:
4369:
4368:
4363:
4358:
4357:
4356:
4344:
4343:
4333:
4320:
4319:
4302:
4300:
4299:
4294:
4292:
4288:
4287:
4286:
4274:
4265:
4264:
4259:
4258:
4246:
4237:
4236:
4231:
4226:
4225:
4224:
4212:
4211:
4201:
4188:
4187:
4170:
4168:
4167:
4162:
4160:
4156:
4155:
4154:
4142:
4133:
4132:
4127:
4126:
4114:
4105:
4104:
4099:
4094:
4093:
4092:
4080:
4079:
4069:
4056:
4055:
4033:
4031:
4030:
4025:
4020:
4019:
4007:
4006:
3994:
3986:
3981:
3980:
3968:
3967:
3950:
3948:
3947:
3942:
3940:
3936:
3935:
3934:
3922:
3913:
3912:
3907:
3902:
3901:
3900:
3888:
3887:
3877:
3864:
3863:
3846:
3844:
3843:
3838:
3836:
3832:
3831:
3830:
3818:
3809:
3808:
3803:
3798:
3797:
3796:
3784:
3783:
3773:
3760:
3759:
3737:
3735:
3734:
3729:
3727:
3726:
3711:
3710:
3698:
3697:
3680:
3678:
3677:
3672:
3670:
3666:
3661:
3660:
3659:
3647:
3646:
3636:
3624:
3623:
3604:
3602:
3601:
3596:
3594:
3593:
3577:
3575:
3574:
3569:
3567:
3566:
3550:
3548:
3547:
3542:
3537:
3529:
3528:
3516:
3515:
3480:
3478:
3477:
3472:
3470:
3469:
3450:
3448:
3447:
3442:
3440:
3439:
3420:
3417:
3414:
3411:
3408:
3405:
3402:
3399:
3396:
3393:
3390:
3387:
3384:
3381:
3378:
3375:
3372:
3369:
3366:
3363:
3360:
3357:
3354:
3351:
3348:
3345:
3342:
3339:
3336:
3333:
3330:
3327:
3324:
3321:
3318:
3315:
3312:
3309:
3306:
3303:
3300:
3297:
3294:
3291:
3288:
3285:
3282:
3279:
3276:
3273:
3270:
3267:
3264:
3261:
3258:
3255:
3252:
3249:
3246:
3243:
3240:
3237:
3234:
3231:
3228:
3225:
3222:
3219:
3216:
3213:
3210:
3207:
3204:
3201:
3198:
3195:
3192:
3189:
3186:
3183:
3180:
3177:
3174:
3171:
3168:
3165:
3162:
3159:
3156:
3153:
3150:
3147:
3144:
3141:
3138:
3135:
3132:
3129:
3126:
3123:
3120:
3117:
3114:
3111:
3108:
3105:
3102:
3099:
3096:
3093:
3090:
3087:
3084:
3081:
3078:
3075:
3072:
3069:
3066:
3063:
3060:
3057:
3054:
3051:
3048:
3045:
3042:
3039:
3036:
3033:
3030:
3027:
3024:
3021:
3018:
3015:
3012:
3009:
3006:
3003:
3000:
2997:
2994:
2991:
2988:
2985:
2982:
2979:
2976:
2973:
2970:
2967:
2964:
2961:
2958:
2955:
2952:
2949:
2946:
2943:
2940:
2937:
2934:
2931:
2928:
2925:
2922:
2919:
2916:
2913:
2910:
2907:
2904:
2901:
2898:
2895:
2892:
2889:
2886:
2883:
2880:
2877:
2874:
2871:
2868:
2865:
2862:
2859:
2856:
2853:
2850:
2847:
2844:
2841:
2838:
2835:
2832:
2829:
2826:
2823:
2820:
2817:
2814:
2811:
2808:
2805:
2802:
2799:
2796:
2793:
2790:
2787:
2784:
2781:
2778:
2775:
2772:
2769:
2766:
2763:
2760:
2757:
2754:
2751:
2748:
2745:
2742:
2739:
2736:
2733:
2730:
2727:
2724:
2721:
2718:
2715:
2712:
2709:
2706:
2703:
2700:
2697:
2694:
2691:
2688:
2685:
2682:
2679:
2676:
2673:
2670:
2667:
2664:
2661:
2658:
2655:
2652:
2649:
2646:
2643:
2640:
2637:
2634:
2631:
2628:
2625:
2622:
2619:
2616:
2613:
2610:
2607:
2604:
2601:
2598:
2595:
2592:
2589:
2586:
2583:
2580:
2577:
2574:
2571:
2568:
2565:
2562:
2559:
2556:
2553:
2550:
2547:
2544:
2541:
2538:
2535:
2532:
2529:
2526:
2522:
2519:
2516:
2512:
2509:
2506:
2503:
2500:
2497:
2494:
2491:
2488:
2485:
2482:
2479:
2476:
2473:
2470:
2467:
2464:
2461:
2458:
2455:
2452:
2449:
2446:
2443:
2440:
2437:
2434:
2431:
2428:
2425:
2422:
2419:
2416:
2413:
2410:
2407:
2404:
2401:
2398:
2395:
2392:
2389:
2386:
2383:
2380:
2377:
2374:
2371:
2368:
2365:
2362:
2359:
2356:
2353:
2350:
2347:
2344:
2341:
2338:
2335:
2332:
2329:
2326:
2323:
2320:
2317:
2313:
2310:
2307:
2304:
2301:
2298:
2295:
2292:
2289:
2286:
2283:
2280:
2277:
2274:
2271:
2268:
2265:
2262:
2259:
2256:
2253:
2250:
2247:
2244:
2241:
2238:
2235:
2232:
2229:
2226:
2223:
2220:
2217:
2214:
2211:
2208:
2205:
2202:
2199:
2196:
2193:
2190:
2187:
2184:
2181:
2178:
2175:
2172:
2169:
2166:
2163:
2160:
2157:
2154:
2151:
2148:
2145:
2142:
2139:
2136:
2133:
2130:
2127:
2124:
2121:
2118:
2115:
2112:
2109:
2106:
2103:
2100:
2097:
2094:
2091:
2088:
2087:% starting point
2078:, and then uses
2077:
2075:
2074:
2069:
2067:
2066:
2042:
2040:
2039:
2034:
2018:
2016:
2015:
2010:
2008:
2007:
1984:
1982:
1981:
1976:
1974:
1973:
1969:
1968:
1967:
1958:
1957:
1944:
1929:
1928:
1907:
1906:
1890:
1888:
1887:
1882:
1877:
1876:
1845:
1843:
1842:
1837:
1835:
1834:
1818:
1816:
1815:
1810:
1805:
1804:
1779:
1777:
1776:
1771:
1769:
1768:
1752:
1750:
1749:
1744:
1732:
1730:
1729:
1724:
1713:
1712:
1700:
1699:
1683:
1681:
1680:
1675:
1673:
1672:
1656:
1654:
1653:
1648:
1643:
1642:
1621:
1613:
1608:
1607:
1591:
1589:
1588:
1583:
1581:
1580:
1564:
1562:
1561:
1556:
1548:
1547:
1520:
1519:
1501:
1499:
1498:
1493:
1485:
1484:
1463:
1462:
1453:
1452:
1442:
1437:
1422:
1414:
1369:
1367:
1366:
1361:
1330:
1325:
1313:
1312:
1281:
1279:
1278:
1273:
1271:
1270:
1239:
1237:
1236:
1231:
1214:
1213:
1205:
1177:
1175:
1174:
1169:
1167:
1158:
1147:
1145:
1144:
1139:
1137:
1128:
1117:
1115:
1114:
1109:
1107:
1098:
1083:
1081:
1080:
1075:
1073:
1064:
1053:
1051:
1050:
1045:
1043:
1038:
1036:
1027:
1021:
1012:
1001:
999:
998:
993:
991:
986:
984:
975:
969:
960:
949:
947:
946:
941:
939:
934:
932:
923:
917:
908:
895:
893:
892:
887:
885:
880:
878:
869:
863:
854:
843:
841:
840:
835:
833:
824:
813:
811:
810:
805:
803:
798:
796:
787:
781:
772:
761:
759:
758:
753:
751:
742:
729:
727:
726:
721:
719:
714:
712:
703:
697:
688:
677:
675:
674:
669:
667:
662:
660:
651:
645:
636:
625:
623:
622:
617:
615:
606:
595:
593:
592:
587:
585:
580:
578:
569:
563:
554:
543:
530:
528:
527:
522:
520:
511:
500:
498:
497:
492:
490:
481:
466:
464:
463:
458:
456:
447:
436:
434:
433:
428:
426:
421:
419:
410:
404:
395:
384:
382:
381:
376:
374:
369:
367:
358:
352:
343:
330:
328:
327:
322:
320:
315:
313:
304:
298:
289:
278:
276:
275:
270:
268:
259:
248:
246:
245:
240:
238:
233:
231:
222:
216:
207:
196:
170:
123:are a family of
104:
101:
95:
93:
52:
24:
16:
4650:
4649:
4645:
4644:
4643:
4641:
4640:
4639:
4625:
4624:
4618:
4596:
4593:
4588:
4587:
4579:
4575:
4567:
4563:
4555:
4551:
4546:
4517:
4501:
4485:
4456:
4443:
4438:
4437:
4410:
4382:
4348:
4335:
4334:
4331:
4327:
4311:
4306:
4305:
4278:
4250:
4216:
4203:
4202:
4199:
4195:
4179:
4174:
4173:
4146:
4118:
4084:
4071:
4070:
4067:
4063:
4047:
4042:
4041:
4011:
3998:
3972:
3959:
3954:
3953:
3926:
3892:
3879:
3878:
3875:
3871:
3855:
3850:
3849:
3822:
3788:
3775:
3774:
3771:
3767:
3751:
3746:
3745:
3718:
3702:
3689:
3684:
3683:
3651:
3638:
3637:
3631:
3615:
3610:
3609:
3585:
3580:
3579:
3558:
3553:
3552:
3520:
3507:
3499:
3498:
3495:
3458:
3453:
3452:
3451:can fall below
3431:
3426:
3425:
3422:
3421:
3418:
3415:
3412:
3409:
3406:
3403:
3400:
3397:
3394:
3391:
3388:
3385:
3382:
3379:
3376:
3373:
3370:
3367:
3364:
3361:
3358:
3355:
3352:
3349:
3346:
3343:
3340:
3337:
3334:
3331:
3328:
3325:
3322:
3319:
3316:
3313:
3310:
3307:
3304:
3301:
3298:
3295:
3292:
3289:
3286:
3283:
3280:
3277:
3274:
3271:
3268:
3265:
3262:
3259:
3256:
3253:
3250:
3247:
3244:
3241:
3238:
3235:
3232:
3229:
3226:
3223:
3220:
3217:
3214:
3211:
3208:
3205:
3202:
3199:
3196:
3193:
3190:
3187:
3184:
3181:
3178:
3175:
3172:
3169:
3166:
3163:
3160:
3157:
3154:
3151:
3148:
3145:
3142:
3139:
3136:
3133:
3130:
3127:
3124:
3121:
3118:
3115:
3112:
3109:
3106:
3103:
3100:
3097:
3094:
3091:
3088:
3085:
3082:
3079:
3076:
3073:
3070:
3067:
3064:
3061:
3058:
3055:
3052:
3049:
3046:
3043:
3040:
3037:
3034:
3031:
3028:
3025:
3022:
3019:
3016:
3013:
3010:
3007:
3004:
3001:
2998:
2995:
2992:
2989:
2986:
2983:
2980:
2977:
2974:
2971:
2968:
2965:
2962:
2959:
2956:
2953:
2950:
2947:
2944:
2941:
2938:
2935:
2932:
2929:
2926:
2923:
2920:
2917:
2914:
2911:
2908:
2905:
2902:
2899:
2896:
2893:
2890:
2887:
2884:
2881:
2878:
2875:
2872:
2869:
2866:
2863:
2860:
2857:
2854:
2851:
2848:
2845:
2842:
2839:
2836:
2833:
2830:
2827:
2824:
2821:
2818:
2815:
2812:
2809:
2806:
2803:
2800:
2797:
2794:
2791:
2788:
2785:
2782:
2779:
2776:
2773:
2770:
2767:
2764:
2761:
2758:
2755:
2752:
2749:
2746:
2743:
2740:
2737:
2734:
2731:
2728:
2725:
2722:
2719:
2716:
2713:
2710:
2707:
2704:
2701:
2698:
2695:
2692:
2689:
2686:
2683:
2680:
2677:
2674:
2671:
2668:
2665:
2662:
2659:
2656:
2653:
2650:
2647:
2644:
2641:
2638:
2635:
2632:
2629:
2626:
2623:
2620:
2617:
2614:
2611:
2608:
2605:
2602:
2599:
2596:
2593:
2590:
2587:
2584:
2581:
2578:
2575:
2572:
2569:
2566:
2563:
2560:
2557:
2554:
2551:
2548:
2545:
2542:
2539:
2536:
2533:
2530:
2527:
2524:
2520:
2517:
2514:
2510:
2507:
2504:
2501:
2498:
2495:
2492:
2489:
2486:
2483:
2480:
2477:
2474:
2471:
2468:
2465:
2462:
2459:
2456:
2453:
2450:
2447:
2444:
2441:
2438:
2435:
2432:
2429:
2426:
2423:
2420:
2417:
2414:
2411:
2408:
2405:
2402:
2399:
2396:
2393:
2390:
2387:
2384:
2381:
2378:
2375:
2372:
2369:
2366:
2363:
2360:
2357:
2354:
2351:
2348:
2345:
2342:
2339:
2336:
2333:
2330:
2327:
2324:
2321:
2318:
2315:
2311:
2308:
2305:
2302:
2299:
2296:
2293:
2290:
2287:
2284:
2282:'ztick'
2281:
2278:
2276:'ytick'
2275:
2272:
2270:'xtick'
2269:
2266:
2263:
2260:
2257:
2254:
2251:
2248:
2245:
2242:
2239:
2236:
2233:
2230:
2227:
2224:
2221:
2218:
2215:
2212:
2209:
2206:
2203:
2200:
2197:
2194:
2191:
2188:
2185:
2182:
2179:
2176:
2173:
2171:lorenz_dynamics
2170:
2167:
2164:
2161:
2158:
2155:
2152:
2149:
2146:
2143:
2140:
2137:
2134:
2131:
2128:
2125:
2122:
2119:
2116:
2113:
2110:
2107:
2104:
2101:
2098:
2095:
2092:
2089:
2086:
2080:Newton's method
2058:
2053:
2052:
2049:
2025:
2024:
1996:
1991:
1990:
1987:Newton's method
1959:
1946:
1920:
1919:
1915:
1898:
1893:
1892:
1859:
1848:
1847:
1826:
1821:
1820:
1819:, we only need
1793:
1782:
1781:
1760:
1755:
1754:
1735:
1734:
1704:
1691:
1686:
1685:
1664:
1659:
1658:
1634:
1599:
1594:
1593:
1572:
1567:
1566:
1539:
1511:
1506:
1505:
1473:
1454:
1444:
1378:
1377:
1304:
1284:
1283:
1262:
1242:
1241:
1198:
1197:
1191:
1151:
1150:
1121:
1120:
1091:
1090:
1057:
1056:
1005:
1004:
953:
952:
901:
900:
847:
846:
817:
816:
765:
764:
735:
734:
681:
680:
629:
628:
599:
598:
547:
546:
504:
503:
474:
473:
440:
439:
388:
387:
336:
335:
282:
281:
252:
251:
200:
199:
163:Butcher tableau
105:
99:
96:
53:
51:
37:
25:
12:
11:
5:
4648:
4646:
4638:
4637:
4627:
4626:
4623:
4622:
4616:
4598:Iserles, Arieh
4592:
4589:
4586:
4585:
4573:
4561:
4548:
4547:
4545:
4542:
4541:
4540:
4529:
4524:
4520:
4516:
4513:
4508:
4504:
4500:
4497:
4492:
4488:
4484:
4481:
4476:
4473:
4468:
4463:
4459:
4455:
4450:
4446:
4435:
4423:
4417:
4413:
4409:
4404:
4400:
4394:
4389:
4385:
4381:
4376:
4372:
4366:
4361:
4355:
4351:
4347:
4342:
4338:
4330:
4326:
4323:
4318:
4314:
4303:
4291:
4285:
4281:
4277:
4272:
4268:
4262:
4257:
4253:
4249:
4244:
4240:
4234:
4229:
4223:
4219:
4215:
4210:
4206:
4198:
4194:
4191:
4186:
4182:
4171:
4159:
4153:
4149:
4145:
4140:
4136:
4130:
4125:
4121:
4117:
4112:
4108:
4102:
4097:
4091:
4087:
4083:
4078:
4074:
4066:
4062:
4059:
4054:
4050:
4035:
4034:
4023:
4018:
4014:
4010:
4005:
4001:
3997:
3992:
3989:
3984:
3979:
3975:
3971:
3966:
3962:
3951:
3939:
3933:
3929:
3925:
3920:
3916:
3910:
3905:
3899:
3895:
3891:
3886:
3882:
3874:
3870:
3867:
3862:
3858:
3847:
3835:
3829:
3825:
3821:
3816:
3812:
3806:
3801:
3795:
3791:
3787:
3782:
3778:
3770:
3766:
3763:
3758:
3754:
3739:
3738:
3725:
3721:
3717:
3714:
3709:
3705:
3701:
3696:
3692:
3681:
3669:
3664:
3658:
3654:
3650:
3645:
3641:
3634:
3630:
3627:
3622:
3618:
3592:
3588:
3565:
3561:
3540:
3536:
3532:
3527:
3523:
3519:
3514:
3510:
3506:
3494:
3491:
3468:
3465:
3461:
3438:
3434:
3371:max_iterations
2945:max_iterations
2085:
2065:
2061:
2048:
2045:
2032:
2006:
2003:
1999:
1972:
1966:
1962:
1956:
1953:
1949:
1943:
1939:
1935:
1932:
1927:
1923:
1918:
1913:
1910:
1905:
1901:
1880:
1875:
1872:
1869:
1866:
1862:
1858:
1855:
1833:
1829:
1808:
1803:
1800:
1796:
1792:
1789:
1767:
1763:
1742:
1722:
1719:
1716:
1711:
1707:
1703:
1698:
1694:
1684:are the roots
1671:
1667:
1646:
1641:
1637:
1633:
1630:
1627:
1624:
1619:
1616:
1611:
1606:
1602:
1579:
1575:
1554:
1551:
1546:
1542:
1538:
1535:
1532:
1529:
1526:
1523:
1518:
1514:
1491:
1488:
1483:
1480:
1476:
1472:
1469:
1466:
1461:
1457:
1451:
1447:
1441:
1436:
1433:
1430:
1426:
1420:
1417:
1412:
1409:
1406:
1403:
1400:
1397:
1394:
1391:
1388:
1385:
1359:
1356:
1353:
1350:
1347:
1344:
1341:
1337:
1334:
1329:
1324:
1320:
1316:
1311:
1307:
1303:
1300:
1297:
1294:
1291:
1269:
1265:
1261:
1258:
1255:
1252:
1249:
1229:
1226:
1223:
1220:
1217:
1211:
1208:
1190:
1187:
1183:
1182:
1179:
1178:
1164:
1161:
1148:
1134:
1131:
1118:
1104:
1101:
1088:
1085:
1084:
1070:
1067:
1054:
1041:
1033:
1030:
1024:
1018:
1015:
1002:
989:
981:
978:
972:
966:
963:
950:
937:
929:
926:
920:
914:
911:
897:
896:
883:
875:
872:
866:
860:
857:
844:
830:
827:
814:
801:
793:
790:
784:
778:
775:
762:
748:
745:
731:
730:
717:
709:
706:
700:
694:
691:
678:
665:
657:
654:
648:
642:
639:
626:
612:
609:
596:
583:
575:
572:
566:
560:
557:
536:
535:
532:
531:
517:
514:
501:
487:
484:
471:
468:
467:
453:
450:
437:
424:
416:
413:
407:
401:
398:
385:
372:
364:
361:
355:
349:
346:
332:
331:
318:
310:
307:
301:
295:
292:
279:
265:
262:
249:
236:
228:
225:
219:
213:
210:
189:
188:
185:
184:
181:
178:
177:
174:
107:
106:
42:. Please help
28:
26:
19:
13:
10:
9:
6:
4:
3:
2:
4647:
4636:
4633:
4632:
4630:
4619:
4613:
4609:
4605:
4604:
4599:
4595:
4594:
4590:
4582:
4577:
4574:
4570:
4565:
4562:
4558:
4553:
4550:
4543:
4522:
4518:
4514:
4511:
4506:
4502:
4498:
4495:
4490:
4486:
4482:
4474:
4471:
4466:
4461:
4457:
4453:
4448:
4444:
4436:
4421:
4415:
4411:
4407:
4402:
4398:
4392:
4387:
4383:
4379:
4374:
4370:
4364:
4359:
4353:
4349:
4345:
4340:
4336:
4328:
4324:
4321:
4316:
4312:
4304:
4289:
4283:
4279:
4275:
4270:
4266:
4260:
4255:
4251:
4247:
4242:
4238:
4232:
4227:
4221:
4217:
4213:
4208:
4204:
4196:
4192:
4189:
4184:
4180:
4172:
4157:
4151:
4147:
4143:
4138:
4134:
4128:
4123:
4119:
4115:
4110:
4106:
4100:
4095:
4089:
4085:
4081:
4076:
4072:
4064:
4060:
4057:
4052:
4048:
4040:
4039:
4038:
4016:
4012:
4008:
4003:
3999:
3990:
3987:
3982:
3977:
3973:
3969:
3964:
3960:
3952:
3937:
3931:
3927:
3923:
3918:
3914:
3908:
3903:
3897:
3893:
3889:
3884:
3880:
3872:
3868:
3865:
3860:
3856:
3848:
3833:
3827:
3823:
3819:
3814:
3810:
3804:
3799:
3793:
3789:
3785:
3780:
3776:
3768:
3764:
3761:
3756:
3752:
3744:
3743:
3742:
3723:
3719:
3715:
3712:
3707:
3703:
3699:
3694:
3690:
3682:
3667:
3662:
3656:
3652:
3648:
3643:
3639:
3632:
3628:
3625:
3620:
3616:
3608:
3607:
3606:
3590:
3586:
3563:
3559:
3538:
3534:
3525:
3521:
3517:
3512:
3508:
3492:
3486:
3482:
3466:
3463:
3459:
3436:
3432:
2083:
2081:
2063:
2059:
2046:
2044:
2030:
2023:expansion in
2022:
2021:Taylor series
2004:
2001:
1997:
1988:
1970:
1964:
1960:
1954:
1951:
1947:
1941:
1937:
1933:
1930:
1925:
1921:
1916:
1911:
1908:
1903:
1899:
1873:
1870:
1867:
1864:
1860:
1853:
1831:
1827:
1801:
1798:
1794:
1787:
1765:
1761:
1740:
1720:
1717:
1709:
1705:
1696:
1692:
1669:
1665:
1639:
1635:
1631:
1628:
1622:
1617:
1614:
1609:
1604:
1600:
1577:
1573:
1544:
1540:
1536:
1530:
1524:
1521:
1516:
1512:
1502:
1489:
1481:
1478:
1474:
1467:
1464:
1459:
1455:
1449:
1445:
1439:
1434:
1431:
1428:
1424:
1418:
1415:
1410:
1404:
1398:
1395:
1389:
1383:
1375:
1373:
1351:
1345:
1339:
1335:
1332:
1327:
1322:
1318:
1314:
1309:
1305:
1301:
1295:
1289:
1267:
1263:
1259:
1253:
1247:
1224:
1218:
1215:
1209:
1206:
1196:
1188:
1186:
1162:
1159:
1149:
1132:
1129:
1119:
1102:
1099:
1089:
1087:
1086:
1068:
1065:
1055:
1039:
1031:
1028:
1022:
1016:
1013:
1003:
987:
979:
976:
970:
964:
961:
951:
935:
927:
924:
918:
912:
909:
899:
898:
881:
873:
870:
864:
858:
855:
845:
828:
825:
815:
799:
791:
788:
782:
776:
773:
763:
746:
743:
733:
732:
715:
707:
704:
698:
692:
689:
679:
663:
655:
652:
646:
640:
637:
627:
610:
607:
597:
581:
573:
570:
564:
558:
555:
545:
544:
541:
540:
539:
515:
512:
502:
485:
482:
472:
470:
469:
451:
448:
438:
422:
414:
411:
405:
399:
396:
386:
370:
362:
359:
353:
347:
344:
334:
333:
316:
308:
305:
299:
293:
290:
280:
263:
260:
250:
234:
226:
223:
217:
211:
208:
198:
197:
194:
193:
192:
182:
180:
179:
175:
172:
171:
168:
167:
166:
164:
160:
155:
153:
148:
146:
142:
138:
134:
130:
126:
122:
118:
114:
103:
92:
89:
85:
82:
78:
75:
71:
68:
64:
61: –
60:
56:
55:Find sources:
49:
45:
41:
35:
34:
33:single source
29:This article
27:
23:
18:
17:
4602:
4583:, p. 47
4581:Iserles 1996
4576:
4571:, p. 63
4569:Iserles 1996
4564:
4559:, p. 47
4557:Iserles 1996
4552:
4036:
3740:
3496:
3423:
2050:
1503:
1376:
1192:
1184:
537:
190:
156:
149:
144:
140:
120:
110:
100:January 2023
97:
87:
80:
73:
66:
54:
30:
4591:References
2936:&&
2518:gauss_step
2156:gauss_step
70:newspapers
4261:−
4129:−
4101:−
3805:−
3464:−
3356:threshold
2990:iteration
2984:iteration
2963:iteration
2939:iteration
2933:threshold
2900:iteration
1938:∑
1871:−
1868:ℓ
1846:to order
1802:ℓ
1741:ℓ
1697:ℓ
1482:ℓ
1440:ℓ
1425:∑
1319:∫
1302:−
1210:˙
1189:Intuition
865:−
699:−
647:−
565:−
300:−
218:−
40:talk page
4629:Category
4600:(1996),
3236:linsolve
3056:dynamics
3005:dynamics
2750:x2_guess
2744:dynamics
2735:x1_guess
2729:dynamics
2672:x2_guess
2618:x1_guess
2606:dynamics
2511:function
2306:function
2243:x_series
2231:x_series
2219:x_series
2201:x_series
2123:x_series
152:A-stable
3230:damping
2951:fprintf
2576:damping
2564:damping
84:scholar
4614:
3380:x_next
3284:k_next
3260:k_next
3221:k_next
2513:x_next
2300:return
1504:where
161:. Its
119:, the
86:
79:
72:
65:
57:
4544:Notes
3359:error
3317:error
2912:while
2882:error
2852:error
2585:error
2579:<=
2445:sigma
2439:sigma
2376:sigma
2364:state
2346:state
2328:state
2314:state
2288:title
2213:plot3
2117:10000
1370:with
1240:with
91:JSTOR
77:books
4612:ISBN
3353:>
3341:norm
2969:norm
2942:<
2930:>
2918:norm
2567:>
2549:sqrt
2531:size
2502:beta
2388:beta
2240:,:),
2228:,:),
2183:1e-7
2105:0.01
176:1/2
165:is:
115:and
63:news
3419:end
3377:end
3335:end
3308:));
3209:a22
3188:eye
3176:a21
3155:a12
3134:a11
3113:eye
3083:a22
3071:a21
3032:a12
3020:a11
2834:a22
2822:sq3
2804:a21
2792:sq3
2774:a12
2756:a11
2699:sq3
2645:sq3
2597:end
2543:sq3
2508:end
2457:rho
2406:rho
2285:,);
2264:gca
2258:set
2252:,:)
2210:end
2195:100
2132:for
173:1/2
111:In
46:by
4631::
4610:,
4606:,
4475:18
4403:15
4399:15
4375:30
4371:15
4271:24
4267:15
4243:24
4239:15
4139:30
4135:15
4111:15
4107:15
3467:12
3460:10
3416:);
3413:k2
3407:k1
3392:dt
3374:);
3347:er
3338:if
3332:);
3329:k2
3323:k1
3311:er
3278:k2
3275:);
3254:k1
3251:);
3248:er
3218:];
3215:j2
3203:dt
3182:j2
3170:dt
3161:j1
3149:dt
3140:j1
3128:dt
3101:);
3098:dt
3089:k2
3077:k1
3050:);
3047:dt
3038:k2
3026:k1
2981:);
2975:er
2924:er
2897:);
2894:k2
2888:k1
2876:er
2867:k2
2861:k1
2858:@(
2753:);
2738:);
2714:dt
2660:dt
2615:);
2573:||
2561:if
2558:);
2540:);
2505:];
2418:td
2412:28
2382:10
2373:);
2355:);
2337:);
2297:);
2279:,,
2273:,,
2255:);
2198:);
2177:dt
2099:dt
2043:.
1374:.
1163:18
1103:18
1069:36
1040:15
1032:15
988:15
980:30
965:36
936:15
928:10
882:15
874:24
859:36
800:15
792:24
777:36
716:15
708:30
693:36
664:15
656:15
611:36
582:15
574:10
183:1
154:.
147:.
4621:.
4528:)
4523:3
4519:k
4515:5
4512:+
4507:2
4503:k
4499:8
4496:+
4491:1
4487:k
4483:5
4480:(
4472:h
4467:+
4462:i
4458:x
4454:=
4449:f
4445:x
4422:)
4416:2
4412:k
4408:h
4393:+
4388:1
4384:k
4380:h
4365:+
4360:2
4354:i
4350:x
4346:+
4341:f
4337:x
4329:(
4325:f
4322:=
4317:3
4313:k
4290:)
4284:3
4280:k
4276:h
4256:1
4252:k
4248:h
4233:+
4228:2
4222:i
4218:x
4214:+
4209:f
4205:x
4197:(
4193:f
4190:=
4185:2
4181:k
4158:)
4152:3
4148:k
4144:h
4124:2
4120:k
4116:h
4096:2
4090:i
4086:x
4082:+
4077:f
4073:x
4065:(
4061:f
4058:=
4053:1
4049:k
4022:)
4017:2
4013:k
4009:+
4004:1
4000:k
3996:(
3991:2
3988:h
3983:+
3978:i
3974:x
3970:=
3965:f
3961:x
3938:)
3932:1
3928:k
3924:h
3919:6
3915:3
3909:+
3904:2
3898:i
3894:x
3890:+
3885:f
3881:x
3873:(
3869:f
3866:=
3861:2
3857:k
3834:)
3828:2
3824:k
3820:h
3815:6
3811:3
3800:2
3794:i
3790:x
3786:+
3781:f
3777:x
3769:(
3765:f
3762:=
3757:1
3753:k
3724:1
3720:k
3716:h
3713:+
3708:i
3704:x
3700:=
3695:f
3691:x
3668:)
3663:2
3657:i
3653:x
3649:+
3644:f
3640:x
3633:(
3629:f
3626:=
3621:1
3617:k
3591:i
3587:x
3564:i
3560:k
3539:2
3535:/
3531:)
3526:i
3522:x
3518:+
3513:f
3509:x
3505:(
3437:i
3433:k
3410:+
3404:(
3401:*
3398:2
3395:/
3389:+
3386:x
3383:=
3368:,
3362:(
3350:)
3344:(
3326:,
3320:(
3314:=
3305:d
3302::
3299:1
3296:(
3293:+
3290:d
3287:(
3281:=
3272:d
3269::
3266:1
3263:(
3257:=
3245:,
3242:j
3239:(
3233:*
3227:-
3224:=
3212:*
3206:*
3200:-
3197:)
3194:d
3191:(
3185:,
3179:*
3173:*
3167:-
3164:;
3158:*
3152:*
3146:-
3143:,
3137:*
3131:*
3125:-
3122:)
3119:d
3116:(
3110:[
3107:=
3104:j
3095:*
3092:)
3086:*
3080:+
3074:*
3068:(
3065:+
3062:x
3059:(
3053:=
3044:*
3041:)
3035:*
3029:+
3023:*
3017:(
3014:+
3011:x
3008:(
3002:=
2999:;
2996:1
2993:+
2987:=
2978:)
2972:(
2966:,
2960:,
2954:(
2948:)
2927:)
2921:(
2915:(
2909:;
2906:1
2903:=
2891:,
2885:(
2879:=
2873:;
2870:)
2864:,
2855:=
2849:;
2846:4
2843:/
2840:1
2837:=
2831:;
2828:6
2825:/
2819:+
2816:4
2813:/
2810:1
2807:=
2801:;
2798:6
2795:/
2789:-
2786:4
2783:/
2780:1
2777:=
2771:;
2768:4
2765:/
2762:1
2759:=
2747:(
2741:=
2732:(
2726:=
2723:;
2720:k
2717:*
2711:*
2708:)
2705:6
2702:/
2696:+
2693:2
2690:/
2687:1
2684:(
2681:+
2678:x
2675:=
2669:;
2666:k
2663:*
2657:*
2654:)
2651:6
2648:/
2642:-
2639:2
2636:/
2633:1
2630:(
2627:+
2624:x
2621:=
2612:x
2609:(
2603:=
2594:)
2588:(
2582:0
2570:1
2555:3
2552:(
2546:=
2537:x
2534:(
2528:=
2525:)
2521:(
2515:=
2499:-
2496:,
2493:x
2490:,
2487:y
2484:;
2481:x
2478:-
2475:,
2472:1
2469:-
2466:,
2463:z
2460:-
2454:;
2451:0
2448:,
2442:,
2436:-
2433:[
2430:=
2427:j
2424:;
2421:=
2415:;
2409:=
2403:;
2400:3
2397:/
2394:8
2391:=
2385:;
2379:=
2370:3
2367:(
2361:=
2358:z
2352:2
2349:(
2343:=
2340:y
2334:1
2331:(
2325:=
2322:x
2316:)
2312:(
2303:;
2291:(
2267:,
2261:(
2249:3
2246:(
2237:2
2234:(
2225:1
2222:(
2216:(
2207:;
2204:=
2192:,
2189:1
2186:,
2180:,
2174:,
2168:@
2165:,
2162:x
2159:(
2153:=
2150:x
2147:N
2144::
2141:1
2138:=
2135:i
2129:;
2126:=
2120:;
2114:=
2111:N
2108:;
2102:=
2096:;
2093:=
2090:x
2064:i
2060:k
2031:h
2005:j
2002:i
1998:a
1971:)
1965:j
1961:k
1955:j
1952:i
1948:a
1942:j
1934:h
1931:+
1926:0
1922:x
1917:(
1912:f
1909:=
1904:i
1900:k
1879:)
1874:1
1865:2
1861:h
1857:(
1854:O
1832:i
1828:k
1807:)
1799:2
1795:h
1791:(
1788:O
1766:i
1762:k
1721:0
1718:=
1715:)
1710:i
1706:r
1702:(
1693:P
1670:i
1666:r
1645:)
1640:i
1636:r
1632:+
1629:1
1626:(
1623:h
1618:2
1615:1
1610:=
1605:i
1601:c
1578:i
1574:w
1553:)
1550:)
1545:i
1541:c
1537:h
1534:(
1531:x
1528:(
1525:f
1522:=
1517:i
1513:k
1490:,
1487:)
1479:2
1475:h
1471:(
1468:O
1465:+
1460:i
1456:k
1450:i
1446:w
1435:1
1432:=
1429:i
1419:2
1416:h
1411:+
1408:)
1405:0
1402:(
1399:x
1396:=
1393:)
1390:h
1387:(
1384:x
1358:)
1355:)
1352:t
1349:(
1346:x
1343:(
1340:f
1336:t
1333:d
1328:h
1323:0
1315:=
1310:0
1306:x
1299:)
1296:h
1293:(
1290:x
1268:0
1264:x
1260:=
1257:)
1254:0
1251:(
1248:x
1228:)
1225:x
1222:(
1219:f
1216:=
1207:x
1160:5
1133:9
1130:4
1100:5
1066:5
1029:1
1023:+
1017:9
1014:2
977:1
971:+
962:5
925:1
919:+
913:2
910:1
871:1
856:5
829:9
826:2
789:1
783:+
774:5
747:2
744:1
705:1
690:5
653:1
641:9
638:2
608:5
571:1
559:2
556:1
516:2
513:1
486:2
483:1
452:4
449:1
423:3
415:6
412:1
406:+
400:4
397:1
371:3
363:6
360:1
354:+
348:2
345:1
317:3
309:6
306:1
294:4
291:1
264:4
261:1
235:3
227:6
224:1
212:2
209:1
145:s
141:s
102:)
98:(
88:·
81:·
74:·
67:·
50:.
36:.
Text is available under the Creative Commons Attribution-ShareAlike License. Additional terms may apply.