Ariketa: Newton 2×2 sistema batean, iterazio-taularekin

exey+xcosy=0e^x e^y+x\cos y=0, x+y=1x+y=1 sistemaren ebazpen osoa Newton-ekin [2,1][2,-1]-etik: iteratuen taula, hondarraren eta inkrementuaren normak, eta ACOC 2an egonkortzen.

Planteamendua

exey+xcosy=0e^xe^y+x\cos y=0, x+y=1x+y=1 sistema ebazten da. Bere eraikuntza (FF eta jacobiarra) dedukzioan egin zen. Hasierako estimazioa x(0)=[2,1]Tx^{(0)}=[2,-1]^T da, tolerantzia 102010^{-20} (doitasun hedatuko aritmetika) eta gehienez 40 iterazio izanik.

Iterazioak

iterx(k)x^{(k)}F(x(k+1))\|F(x^{(k+1)})\|x(k+1)x(k)\|x^{(k+1)}-x^{(k)}\|ACOC
1[5.3247, −4.3247]0.70514.7018n/a
2[5.1697, −4.1697]0.04860.2191n/a
3[5.1573, −4.1573]0.00030.01760.8229
4[5.1572, −4.1572]1.7091·10⁻⁸0.00011.9620
5[5.1572, −4.1572]4.3406·10⁻¹⁷6.2690·10⁻⁹1.9989
6[5.1572, −4.1572]2.7997·10⁻³⁴1.5921·10⁻¹⁷2.0000
Newton-en iterazioak x(0)=[2,1]x^{(0)}=[2,-1]-etik; soluzioa α[5.15723,4.15723]\alpha\approx[5.15723,\,-4.15723] da.