Gustavo Alberto Ochoa Diaz

Análisis Numérico
1) Calcule el error absoluto y relativo en las aproximaciones de p mediante 𝑝∗., para:

a) p=e10 ; p∗¿ 22000

b) p=8! ; p∗¿ 39900

Error absoluto = |𝑝 − 𝑝∗|

Error relativo = ¿ p− p∗¿ /| p|

A) Si p=e10 y p∗¿ 22000, entonces:

Error absoluto = ¿ e 10−22000∨¿

Error absoluto = 22025.4657948 − 22000

Error absoluto = 25.4657948

Error relativo = ¿ e 10−22000∨¿∨e10∨¿

Error relativo = 25.4657 / 22026.465

Error relativo = 0.001156

2) Realice los siguientes cálculos (i) en forma exacta a CINCO cifras, (ii) mediante una

aproximación a CINCO cifras y (iii) con una aritmética de truncamiento a CINCO cifras, (iv)

calcule los errores relativos en los incisos (ii) y (iii), comparando con la solución exacta de I.

5 2
( 4
a. 3 − 13 + 9 )
MATLAB

syms x

f= ( 53 − 132 )+ 49 ;
(i)

exacta = vpa(f, 5);

exacta =1.9573
(ii)

>> aprox = round(exacta, 5)

aprox = 1.95726

(iii)

fix ( exacta∗105 )
trunc= 5
10
trunc = 1.9573

(iv)

>> error_aproximado = |(exacta−aprox)|/exacta;

>> error_truncamiento = |(exacta−trunc)|/exacta ;

>> error_aproximado

error_aproximado = 0.00000253275

>> error_truncamiento

error_truncamiento = 0.0000179039

1 2 1
(
b. 5 – 7 – 4 . )
(i)

>> syms x

f = (1/5 - 2/7) - 1/4;

exacta = vpa(f, 5)

exacta = -0.33571

(ii)

>> aprox = round(exacta, 5)

aprox = -0.33571
(iii)

>> trunc = fix(exacta∗105 )/ 105

trunc = -0.3357

(iv)

>> error_aproximado = |(exacta−aprox)|/exacta

error_aproximado = -0.00001276

>> error_truncamiento = |(exacta−trunc)|/exacta

error_truncamiento = -0.00004255

3) Encuentre en cada caso la aproximación que se pide, utilizando el teorema de Taylor y estime
el error.
a. 𝑪𝒐𝒔 (𝟔𝟏º) 𝒑𝒐𝒍𝒊𝒏𝒐𝒎𝒊𝒐 𝒅𝒆 𝒈𝒓𝒂𝒅𝒐 9
>> syms x ;

f =cos (x);
P 9=taylor (f , x ,' Order ' , 10);
>> x=61∗pi/180 ;

cos 61=double(subs (P 9 , x ));

>> cos 61

cos 61=0.484 8
Error absoluto de cos (61°) está dado por:

¿ R 10(x )∨¿ 1∗¿ 61∗pi/180−0 ∨(10) ¿(10)!≤2.2105e-1 4

b. √ 82 𝒑𝒐𝒍𝒊𝒏𝒐𝒎𝒊𝒐 𝒅𝒆 𝒈𝒓𝒂𝒅𝒐 5
3

1
3
f (x)=x
−2
f ' (x)=()
1
3
x 3
 −5
f ' ' (x)= ( )
−2
9
x 3

−8
f ' ' ' ( x)= ( )
10
27
x 3

−11
f ' ' ' ' (x)= ( )
−80
81
x 3

−14
f ' ' ' ' ' (x )=( )
440
243
x 3

P 5(x)=f ( 82)+ f '(82)( x−82)+ ( 12 ) f ( 82) ( x−82) +( 16 ) f ( 82 )( x−82) +( 241 ) f ( 82) ( x−82) +( 1201 ) f ( 82)
'' 2 '' ' 3 '' ' ' 4 ' '' ' '

P 5(x)=4 +(
246 )
( x−82 )−(
29,568 )
( x−82 ) + (
1,120,896 )
( x−82 ) −(
477,757,696 )
( x−82 ) +(
1 1 5 2 40 3 4 440
2,213,837,312

Rn(x)= ( n+11 ) f ( c )( x−a )

( n +1) n+1

R 5(x )=( ) f ( c )( x−82 )

1 (6 ) 6
6
−17
(6)
f (x )=
−20
81
x ( ) 3

−17
f (6) (x )= ( )
−20
81
x 3

P 5(∛ 82)≈ 3.999998

el error de la aproximación es menor que 1.69 x 10−6

4) Ejercicios: Calcula los polinomios de Taylor de orden 2, 3 y 4 alrededor del punto a para las
siguientes funciones f.

a. 𝒇(𝒙) = 𝒍𝒐𝒈(𝒙 + 𝟏); 𝒂 = 0

Polinomio de Taylor de orden 2:

'
P 2 ( x )=f ( a )+ f ( a ) ( x −a ) + ( 12 ) f ( a )( x−a)
'' 2
Polinomio de Taylor de orden 3:

P 3 ( x ) =f ( a )+ f ' ( a ) ( x−a )+ ( 12 ) f ( a )( x−a) +( 16 ) f

'' 2 ' ''
( a ) ( x−a )3

Polinomio de Taylor de orden 4:

P 4 ( x )=f ( a ) + f ' ( a ) ( x −a ) + ( 12 ) f ( a )( x−a) +( 16 ) f

'' 2 '' '
( a ) ( x−a )3 + ( 241 ) f '' ''
( a ) ( x−a )4

f ( x )=log ( x+1 )

' 1
f ( x )=
x+1
'' −1
f ( x )=
( x +1 )2
'' ' 2
f ( x )=
( x +1 )3
'' ' ' −6
f ( x )=
( x+1 )4
Entonces, evaluando en el punto a = 0, se tiene:

f ( 0 )=log ( 1 )=0
' 1
f ( 0 )= =1
1
'' −1
f ( 0 )= 2
=−1
1
'' ' 2
f ( 0 )= 3
=2
1
'' ' ' −6
f (0)= 4
=−6
1
Por lo tanto, los polinomios de Taylor son:

Polinomio de Taylor de orden 2:

()
2
1 ( )( 2 x
P 2 ( x )=0+1 ( x −0 ) + −1 x−0 ) =x −
2 2
Polinomio de Taylor de orden 3:
 () ()
2 3
1 ( )( 2 1 ( 3 x x
P 3 ( x ) =0+1 ( x−0 ) + −1 x−0 ) + 2 x −0 ) =x− +
2 6 2 3
Polinomio de Taylor de orden 4:

() () ( )
2 3 4
1 ( )( 2 1 ( 3 1 ( )( 4 x x x
P 4 ( x )=0+1 ( x −0 ) + −1 x−0 ) + 2 x−0 ) + −6 x−0 ) =x− + −
2 6 24 2 3 4

b. 𝒇(𝒙) = √𝒙 + 𝟏; 𝒂 = 𝟎.
Polinomio de Taylor de orden 2:

'
P 2 ( x )=f ( a )+ f ( a ) ( x −a ) + ( 12 ) f ( a )( x−a)
'' 2

Polinomio de Taylor de orden 3:

P 3 ( x ) =f ( a )+ f ' ( a ) ( x−a )+ ( 12 ) f ( a )( x−a) +( 16 ) f

'' 2 ' ''
( a ) ( x−a )3

Polinomio de Taylor de orden 4:

P 4 ( x )=f ( a ) + f ' ( a ) ( x −a ) + ( 12 ) f ( a )( x−a) +( 16 ) f

'' 2 '' '
( a ) ( x−a )3 + ( 241 ) f '' ''
( a ) ( x−a )4

f ( x )= √ x +1
' 1
f ( x )=
2∗√ x +1
'' 1
f ( x )= 3
4∗( x +1 ) 2
'' ' 3
f ( x )= 5
8∗( x+ 1 ) 2
'' ' ' −15
f ( x )= 7
16∗ ( x +1 ) 2

Si evaluamos en el punto a = 0:

' 1 1
f ( 0 )=√ x+ 1=1 f ( 0 )= =
2∗√ x+1 2
'' 1 −1
f ( 0 )= =
3
4
4∗ ( x +1 ) 2
 '' ' 3 3
f ( 0 )= =
5
8
8∗( x +1 ) 2

'' ' ' −15 −15

f (0)= =
7
16
16∗( x +1 ) 2

Por lo tanto, los polinomios de Taylor son:

Polinomio de Taylor de orden 2:

( )( )
2
1 1 −1 ( 2 x x
P 2 ( x )=1+ ( x−0 ) + x−0 ) =1+ −
2 2 4 2 8
Polinomio de Taylor de orden 3:

( )( ) ()
2 3
1 1 −1 ( 2 3 ( 3 x x x
P 3 ( x ) =1+ ( x −0 ) + x−0 ) + 2 x−0 ) =1+ − −
2 2 4 8 2 8 16

Polinomio de Taylor de orden 4:

( )( ) () ( )
2 3
1 1 −1 ( 2 3 ( 3 −15 ( ) ( 4 x x x
P 4 ( x )=1+ ( x−0 ) + x−0 ) + 2 x −0 ) + −6 x −0 ) =1+ − −
2 2 4 8 16 2 8 16

5) Encuentre en cada caso un polinomio que satisface:

a. 𝑷(𝟎) = 𝟕, 𝑷′(𝟎) = 𝟑, 𝑷′′(𝟎) = 𝟖, 𝑷′′′(𝟎) = 𝟓4

'' '' '
f ( c ) ( x−c ) 2 f ( c ) ( x−c ) 3
P ( x )=f ( c ) + f ' ( c ) ( x−c ) + + +⋯
2! 3!
8 x 2 54 x 3
P ( x )=7+3 x + +
2! 3!
P ( x )=7+3 x +4 x 2+9 x 3

𝑷(𝟏) = 𝟏, 𝑷′(𝟏) = 𝟓, 𝑷 ′′(𝟏) = 𝟑𝟐, 𝑷′′′(𝟏) = 𝟒𝟐

3 2
P(x )=( x−1 ) +5 ( x−1 ) +16(x−1)+42
 3 2
P(1)=( 1−1 ) +5 ( 1−1 ) +16 (1−1)+42=1
2
P '(x )=3 ( x−1 ) +10(x−1)+16
2
P '(1)=3 ( 1−1 ) +10(1−1)+16=5
P ' '(x )=6 (x−1)+10
P ' '(1)=6 (1−1)+10=32
P ' ' '(x )=6
P ' ' '(1)=6
3 2
P(x )=( x−1 ) +5 ( x−1 ) +16(x−1)+42

6. Desarrolla los siguientes ítems.

a. Expresa el polinomio 𝒑(𝒙) = 𝒙 𝟒 − 𝟑𝒙 𝟐 + 𝟓𝐱 − 𝟏 en potencia de (𝒙 − 𝟏).

4 3 2
P ( x )=a ( x−1 ) + b ( x−1 ) +c ( x−1 ) + d ( x−1 ) + c
' 3 3 2
P ( x )=4 x −6 x +5=4 x ( x−1 ) +3 b ( x−1 ) + 2 c ( x−1 ) +d
'' 2 2
P ( x )=12 x −6=12a ( x−1 ) +6 b ( x −1 )+ 2 c
' ''
P ( x )=24 a ( x−1 ) +6 b
' '' '
P ( x )=24=24 a
24
a= =1 ; b=4 ; c=3 ; d=3 ; c=2
24
4 3 2
P ( x )=( x−1 ) + 4 ( x−1 ) +3 ( x−1 ) +3 ( x−1 ) +2

b. Expresa el polinomio 𝒑(𝒙) = 𝒙 𝟓 + 𝟓𝒙 𝟒 − 𝟑𝒙 𝟑 + 𝟐𝒙 𝟐 − 𝟒𝒙 + 𝟓 en potencias de (𝒙 + 𝟏)

p(x )=x 5+5 x 4−3 x 3+ 2 x 2−4 x+5

p(x +1)=(x +1)5+ 5(x+ 1)4−3 ( x+1)3+2(x +1)2−4( x+1)+5
p(x +1)=x 5+5 x 4 +10 x 3+10 x 2+5 x +5

7) Encuentre el o los intervalos donde existe al menos una raíz en a cada una de las siguientes
ecuaciones. (Describa muy bien el intervalo que se escoja para mostrar la raíz)
a) f (x)=cos ( x + √ 2 )−x ( 2x +√ 2)=0
Para encontrar los intervalos donde existe al menos una raíz, podemos utilizar el método de la
bisección. Empecemos por considerar el intervalo [-2, -1]:

f (−2)=cos (−2+ √ 2 ) −(−2 ) ( −22 + √2)=−0.1987

f (−1)=cos (−1+ √ 2 ) −(−1 ) ( + √ 2 )=0.705 2
−1
2
xi+ xi−1
xi+1=
2

Obteniendo:

−2+ (−1 )
x 2= =−1.5
2
−1+ (−1.5 )
x 3= =−1.2 5
2
−1.5 + (−1.25 )
x 4= =−1.37 5
2
−1.25+ (−1.375 )
x 5= =−1.312 5
2
la raíz de la ecuación se encuentra en el intervalo (-1.3125, -1.25), y aproximadamente tiene un
valor de x ≈ -1.281.

b) f (x)=3 ln ( 2+ x ) +2 x 3−x 2−2 x – 20

Podemos empezar por considerar el intervalo [-4, -3]:

3 2
f (−4 )=3 ln(2−4 )+ 2 (−4 ) −(−4 ) −2(−4)−20=−107.560 6
3 2
f (−3)=3 ln (2−3)+2 (−3 ) −(−3 ) −2(−3)−20=2.428 2
xi+ xi−1
xi+1=
2
Obteniendo:

−4+ (−3 )
x 2= =−3. 5
2
−3+ (−3.5 )
x 3= =−3.2 5
2
−3.5+ (−3.25 )
x 4= =−3.37 5
2
−3.25+ (−3.375 )
x 5= =−3.3125
2
una de las raíces de la ecuación se encuentra en el intervalo (-3.3125, -3.25), y aproximadamente
tiene un valor de x ≈ -3.28. Raíz cerca de x ≈ -0.8: Intervalo [-1, 0]

Raíz cerca de x ≈ 1.5: Intervalo [1, 2]

8) Aplique los métodos siguientes para obtener una solución con exactitud de 𝟏𝟎−𝟒 para los
problemas: Grafique, escoja el intervalo donde sea necesario hacerlo.

a. f ( x )=e x −x 2+ 3 x−2

>> x = linspace(-2, 3, 1000);

>> y = exp(x) - x^2 + 3*x - 2;

Plot(x,y)

La gráfica muestra que hay una raíz entre x=0 y x=1.5:

−4
a=0 ; b=1.5 ; tol=10
 n a xi b f ( a) f ( xi ) f (b )
1 0 1.5 4.7317 -1
2 0.75 1.125 1.125 1.8045 3.1896 4.7317
3 0.9375 0.9375 1.125 2.4872 3.1896 2.8348
4 0.9375 1.0312 1.031 2.4872 2.6598 2.8339
5 0.9843 0.9843 1.031 2.6600 2.7467 2.8339
6 0.9843 1.007 1.007 2.6600 2.7021 2.7443
7 0.9960 1.0019 1.007 2.7034 2.7239 2.7443

La solución obtenida es x ≈ 0.9996

b. f ( x )=1−0 , 6 x /x