Page 216 - Chapra y Canale. Metodos Numericos para Ingenieros 5edición_Neat

192 RAÍCES DE POLINOMIOS

Solución. Bajo las mismas condiciones iniciales del ejemplo 5.6, se usa MATLAB
para determinar la raíz positiva.
>> x0=[0 1.3];
>> x=fzero(inline(‘x^10–1’),x0)

x =
1

De manera semejante, se emplean los valores iniciales –1.3 y 0 para determinar la
raíz negativa

>> x0=[–1.3 0];
>> x=fzero(inline(‘x^10–1’),x0)
x =
–1
Se puede usar un valor único; resulta un caso interesante cuando se usa el valor
inicial 0
>> x0=0;
>> x=fzero(inline(‘x^10–1’),x0)

x =
–1

Se tiene que para ese valor el algoritmo llevará a la raíz a su valor negativo.
El uso de optimset se ilustra al mostrar en pantalla la forma en que las iteraciones
conducen a la solución
>> x0=0;
>> option=optimset(‘DISP’,’ITER’);
>> x=fzero(inline(‘x^10–1’),x0,option)

Func–count x f(x) Procedure
1 0 –1 initial
2 –0.0282843 –1 search
3 0.0282843 –1 search
4 –0.04 –1 search
•
•
•
21 0.64 –0.988471 search
22 –0.905097 –0.631065 search
23 0.905097 –0.631065 search
24 –1.28 10.8059 search

Looking for a zero in the interval [–1.28], 0.9051]

25 0.784528 –0.911674 interpolation
26 –0.247736 –0.999999 bisection
27 –0.763868 –0.932363 bisection

6/12/06 13:51:28
Chapra-07.indd 192 6/12/06 13:51:28
Chapra-07.indd 192