Bisektio-esimerkki-istunto
Testiesimerkkinä katsomme viimeistä kommenttiesimerkkiä. Rajoitamme iteraatioiden
lukumäärän kymmeneen tilan säästämiseksi.
format long % Tulostustarkkuus maksimiarvoonsa
f=inline('exp(x)-1.5-atan(x)','x')
ab=bisekt(f,[-2,5],0.00001,10)
format short % Takaisin oletusarvoon
virheet=diff(ab') % Muodostetaan välien pituudet.
[(1:10)' ab] % Laitetaan iteraatiokierrosnumero 1. sarakkeeksi.
ab =
-2.00000000000000 5.00000000000000
-2.00000000000000 1.50000000000000
-0.25000000000000 1.50000000000000
0.62500000000000 1.50000000000000
0.62500000000000 1.06250000000000
0.62500000000000 0.84375000000000
0.73437500000000 0.84375000000000
0.73437500000000 0.78906250000000
0.76171875000000 0.78906250000000
0.76171875000000 0.77539062500000
virheet =
Columns 1 through 7
7.0000 3.5000 1.7500 0.8750 0.4375 0.2188 0.1094
Columns 8 through 10
0.0547 0.0273 0.0137
ans =
1.0000 -2.0000 5.0000
2.0000 -2.0000 1.5000
3.0000 -0.2500 1.5000
4.0000 0.6250 1.5000
5.0000 0.6250 1.0625
6.0000 0.6250 0.8438
7.0000 0.7344 0.8438
8.0000 0.7344 0.7891
9.0000 0.7617 0.7891
10.0000 0.7617 0.7754