TAREA DE NUMÉRICO:

1) Encuentre los polinomios de Mclaurin para cos(x). Calcule

después el valor de cos (0.2), empleando n =7.
RESOLUCIÓN:
f '' ( 0 ) x2 f ' ' ' ( 0 ) x3 f IV ( 0 ) x 4 f V ( 0 ) x 5 f VI ( 0 ) x 6 f VII ( 0 ) x 7
f ( x )=f ( 0 ) + f ' ( 0 ) ( x )+ + + + + +
2! 3! 4! 5! 6! 7!

f ( x )=cos ( X )=cos ( 0 )=1

f ' (x )=−sen ( X ) =−sen (0)=0

f ' ' ( x)=−cos ( X )=−cos ( 0 )=−1

f ' ' ' (x)=sen ( X )=sen (0)=0

f IV ( x)=cos ( X )=cos(0)=1

f V ( x)=−sen ( X )=−sen(0)=0

f VI (x)=−cos ( X ) =−cos ( 0 ) =−1

f VII ( x)=sen ( X )=sen(0)=0

(−1) X 2 (1) X 4 (−1) X 6

cos(X) = 1+ + +
2! 4! 6!

X2 X4 X6
cos(0.2) = 1− + −
2! 4! 6!

0.22 0.24 0.26

cos(0.2) = 1− + −
2! 4! 6!

cos(0.2) = 0.98
 2) Usar Mclaurin para ⅇ x con n=6, ⅇ 0.2 :

f(X) = ⅇ x f(0) = ⅇ 0 = 1

f’(X) = ⅇ x f’(0) = ⅇ 0 = 1

f’’ (X) = ⅇ x f’’ (0) = ⅇ 0 = 1

f’’’(X) = ⅇ x f’’’ (0) = ⅇ 0 = 1

fIV(X) = ⅇ x fIV(0) = ⅇ 0 = 1

fV(X) = ⅇ x fV(0) = ⅇ 0 = 1

fVI(X) = ⅇ x fVI(0) = ⅇ 0 = 1

f '' ( 0 ) x2 f ' ' ' ( 0 ) x3 f IV ( 0 ) x 4 f V ( 0 ) x 5 f VI ( 0 ) x 6

f ( x )=f ( 0 ) + f ' ( 0 ) ( x )+ + + + +
2! 3! 4! 5! 6!

x x2 x3 x4 x5 x6
ⅇ =1+ X + + + + +
2 ! 3! 4 ! 5 ! 6 !

0.2 0.22 0.23 0.2 4 0.25 0.26

ⅇ =1+0.2+ + + + +
2! 3! 4 ! 5! 6!

ⅇ 0.2=1.221402

3) Encuentre P3(X) =f(X) = ⅇ 3 X; para a = 1:

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

f ( x )=f ( a ) + f ' ( a )( x−1 ) + +
2! 3!

f(X) = e5X =f(X)= e5 f(1) = 148.413

f’(X) = 5e5X f’(1) = 742.066

f’’(X) = 52e5X. f''(1) = 3710.329

f’’’(X) = 53e5X. f'''(1) = 18551.645

f IV(X) = 54e5X. fIV(1) = 92758.224

3710.329( x−1)2 18551.645(x−1)3

P3(X) = 148.413+742.066 ( X −1 )+ +
2! 3!

4) Convertir 2768 a base 10:

2768 = 2*82 + 7*81 + 6*80

2768 = 128 + 56 + 6
2768 = 190

5) Encuentre la raíz de la sgte ecuación mediante el método de

Bisección, f(X)= 3X2 -2X -10 en el intervalo de [-2.5; -1].

 [-2.5; -1]
a0 ; b0
a0 +b 0 −2.5−1
 x 0= 2
=
2
=−1.750

f(a 0)*f( x 0)>0

a 1=x 0=−1.750

b 1=b0 =−1

 [-1.750; -1]
a1 ; b1
a1 +b1 −1.750−1
 x 1= 2
=
2
=−1.375
f(a 1)*f( x 1)<0 a 2=−1.750
b 2=−1.375

 [-1.750; -1.375]
a 2; b 2
a2 +b2 −1.750−1.375
 x 2= 2
=
2
=−1.563

f(a 2)*f( x 2)>0 a 3=−1.563

b 3=−1.375

 [-1.563; -1.375]
a 3; b 3
a3 +b 3 −1.563−1.375
 x 3= 2
=
2
=−1.469

f(a 3)*f( x 3)<0 a 4=−1.563

b 4=−1.469

 [-1.563; -1.469]
a4; b4
a 4 +b 4 −1.563−1.469
 x4 = 2
=
2
=−1.516

f(a 4)*f( x 4)<0 a 5=−1.563

b 5=−1.516

n Xn |Xn+1 - Xn |
0 -1.750
1 -1.375 0.375
2 -1.563 0.188
3 -1.469 0.094
4 -1.516 0.047≤0.05
 La raíz cuyo valor estimado es menor o igual a 0.05 es x 4 =−1.516 .

6) Realizar 6 iteraciones para esta función, usando Regula Falsi,

f(X)=sen(X) – X+1 en un intervalo de [1;3]. (6 dígitos)

CONDICIÓN: f(1)*f(3)< 0 entonces sí tiene raíz

 [1;3]
a 0 ; b0

f (b¿¿ 0)a 0−f (a 0)b0

x 0= ¿
f ( b 0 )−f (a 0)

(−1.858880∗1)−( 0.841771∗3)
x 0=
(−1.858880 )−(0.841471)

x 0=1.623231

f(a 0)*f( x 0)>0

a 1=1.623231

b 1=3

 [1.623231;3]
a 1 ; b1
f (b¿¿ 1) a1−f ( a1) b1
x 1= ¿
f ( b1 ) −f ( a1)

(1.623231∗−1.858880)−(0.375395∗3)
x 1=
(−1.858880 )−(0.375395)
 x 1=1.854551

f(a 1)*f( x 1)>0

a 2=1.854551

b 2=3

 [1.854551;3]
a 2 ; b2
f (b¿¿ 2) a2 −f ( a2) b2
x 2= ¿
f ( b2 ) −f (a2 )

(−1.858880∗1.854551)−(0.105460∗3)
x 2=
(−1.858880 )−(0.105460)

x 2=1.916047

f(a 2)*f( x 2)>0

a 3=1.916047

b 3=3

 [1.916047;3]
a3; b3
f (b¿¿ 3)a 3−f (a3 )b3
x 3= ¿
f ( b 3 )−f (a3 )

(−1.858880∗1.916047)−( 0.024944∗3)
x 3=
(−1.858880 )−(0.024944)

x 3=1.850953

f(a 3)*f( x 3)>0

a 4=1.850953

b 4=3

 [1.850953;3]
a 4 ; b4
f (b¿¿ 4) a4 −f (a4 )b4
x4 = ¿
f ( b 4 ) −f ( a4 )

(−1.858880∗1.850953)−(0.110059∗3)
x4 =
(−1.858880 )−(0.110059)

x 4 =1.914666

f(a 4)*f( x 4)>0

a 5=1.914666

b 5=3

 [1.914666;3]
a5; b5
f (b¿¿ 5)a 5−f (a5 )b5
x 5= ¿
f ( b 5 )−f (a5 )

(−1.858880∗1.914666)−(1.914666∗3)
x 5=
(−1.858880 )−(1.914666)

x 5=2.465355

n Xn
0 1.623231
 1 1.854551
2 1.916047
3 1.850953
4 1.914666
5 2.465355

Con ello se realizó las 6 iteraciones.

7) Usando el método de falsa posición, para la función f(X)=cos 3(X)

√ 2 x−3 , X ∈ [2;9], iterar 6 veces.

CONDICIÓN: f(2)*f(9)< 0 entonces sí tiene raíz

 [2;9]
a 0 ; b0

f (b¿¿ 0)a 0−f (a 0)b0

x 0= ¿
f ( b 0 )−f (a 0)

(1.194680∗2)−(−3.2910∗9)
x 0=
( 1.194680 )−(−3.2910)

x 0=7.135676

f(a 0)*f( x 0)>0

a 1=7.135676

b 1=9

 [7.135676;9]
a 1 ; b1
 f (b¿¿ 1) a1−f ( a1) b1
x 1= ¿
f ( b1 ) −f ( a1)

(1.194680∗7.135676)−(−1.257240∗9)
x 1=
( 1.194680 )−(−1.257240)

x 1=8.091622

f(a 1)*f( x 1)>0

a 2=8.091622

b 2=9

 [8.091622;9]
a 2 ; b2

f (b¿¿ 2) a2 −f ( a2) b2
x 2= ¿
f ( b2 ) −f (a2 )

(1.194680∗8.091022)−(−4.696365∗9)
x 2=
( 1.19468 )−(−4.696365)

x 2=8.815785

f(a 2)*f( x 2)<0

a 3=8.091622

b 3=8.815785

 [8.091622;8.815785]
a3; b3
 f (b¿¿ 3)a 3−f (a3 )b3
x 3= ¿
f ( b 3 )−f (a3 )

(1.037141∗8.091622)−(−4.696365∗8.815785)
x 3=
( 1.037141 )−(−4.696365)

x 3=8.684790

f(a 3)*f( x 3)>0

a 4=8.684790

b 4=8.815785

 [8.68470;8.815785]
a 4 ; b4

f (b¿¿ 4) a4 −f (a4 )b4

x4 = ¿
f ( b 4 ) −f ( a4 )

(1.037141∗8.684790)−(−3.138473∗8.815785)
x4 =
(1.037141 ) −(−3.138473)

x 4 =8.783248

f(a 4)*f( x 4)>0

a 5=8.783248

b 5=¿8.815785
  [8.783248; 8.815785]
a5; b5

f (b¿¿ 5)a 5−f (a5 )b5

x 5= ¿
f ( b 5 )−f (a5 )

(1.037141∗8.783248)−(−0.722901∗8.815785)
x 5=
( 1.037141 )−(−0.722901)

x 5=8.796612

n Xn
0 7.135676
1 8.091622
2 8.815785
3 8.684790
4 8.783248
5 8.796612

Con esto se realizó las 6 iteraciones.

8) Resolver la sgte función por el método de falsa posición

f(X)=X2|sen√ X | - 5, calcular ES=0.5% con 5 dígitos. X ∈ [0;5]

CONDICIÓN: f(0)*f(5) entonces sí tiene raíz

 [0;5]
a 0 ; b0
 f (b¿¿ 0)a 0−f (a 0)b0
x 0= ¿
f ( b 0 )−f (a 0)

(14.66873∗0)−(−5∗5)
x 0=
( 14.66873 )−(−5)

x 0=1.27105
f(a 0)*f( x 0)>0
a 1=1.27105
b 1=5

 [1.27105;5]
a 1 ; b1
f (b¿¿ 1) a 1−f (a 1)b 1
x 1= ¿
f ( b1 )−f (a1)

(14.66873∗1.27105)−(−3.54065∗5)
x 1=
( 14.66873 )−(−3.54065)

x 1=1.99611
f(a 1)*f( x 1)>0
a 2=1.99611
b 2=5

 [1.99611;5]
a 2 ; b2
f (b¿¿ 2) a 2−f ( a2)b2
x 2= ¿
f ( b2 ) −f (a2)

(14.66873∗1.99611)−(−1.06515∗5)
x 2=
( 14.66873 ) −(−1.06515)

x 2=2.19370
f(a 2)*f( x 2)>0
a 3=2.19370
b 3=5

 [2.19370;5]
a3; b3
f (b¿¿ 3)a 3−f (a 3)b 3
x 3= ¿
f ( b3 ) −f (a3 )

(14.66873∗2.19370)−(−0.20702∗5)
x 3=
( 14.66873 )−(−0.20702)

x 3=2.23275
f(a 3)*f( x 3)>0
a 4=¿ 2.23275
b 4=5
 [2.23275;5]
a 4 ; b4
f (b¿¿ 4) a 4−f (a 4 )b 4
x4 = ¿
f ( b4 ) −f (a 4)

(14.66873∗2.23275)−(−0.02943∗5)
x4 =
( 14.66873 ) −(−0.02943)

x 4 =2.23829
f(a 4)*f( x 4)>0
a 5=2.23829
b 5=5

 [2.23829;5]
a5; b5
f (b¿¿ 5)a 5−f (a 5)b 5
x 5= ¿
f ( b5 ) −f (a5 )

(14.66873∗2.23829)−(−0.00403∗5)
x 5=
( 14.66873 )−(−0.00403)
 x 5=2.23905

n Xn Xn+1-Xn
0 1.27105
1 1.99611 0.72506
2 2.19370 0.19759
3 2.23275 0.03905
4 2.23829 0.00554
5 2.23905 0.00076<0.005

La raíz cuyo valor estimado es 0.005, es X5= 2.23905.

9) Resolver por método punto fijo e indicar un error aprox. % ≤
0.01, f(X)=sen(√ X )-X ,(6 dígitos)
V0=0.5
usar la sgte función de iteración: g(X)=sen(√ X )

Solución:
0.5=X0=g(X0)=sen(√ X )
X1=g(0.5)=sen(√ 0.5)=0.649637
X2=g(0.649637)=sen(√ 0.649637)=0.721524
X3=g(0.721524)=sen(√ 0.721524 )=0.750901
X4=g(0.750901)=sen(√ 0.750901)=0.762097
X5=g(0.762097)=sen(√ 0.762097)=0.766248
X6=g(0.766248)=sen(√ 0.766248)=0.767772
X7=g(0.767772)=sen(√ 0.767772)=0.768329
X8=g(0.768329)=sen(√ 0.768329)=0.768532
X9=g(0.768532)=sen(√ 0.768532)=0.768606

n Xn EA%
0 0.5
1 0.649637 23%
2 0.721534 9.95%
3 0.750901 3.91%
4 0.762097 1.46%
5 0.766248 0.54%
 6 0.767772 0.19%
7 0.768329 0.07%
8 0.768532 0.02%
9 0.768606 0.009%≤0.01%

La raíz con EA%≤0.01% , es X9=0.768606.

10) Resolver por punto fijo con exactitud estimada de 10 -2 la

sgte función f(X)= X4-3X-3en el intervalo [1;2].
X 4 −3 X 2−3
Función de iteración: g(X)= X −
4 X 3−6 X

Por teorema:
g(1)=-1.5∄[1;2]
g(2)=1.95 ∈[1; 2]
Ambas imágenes de los puntos del intervalo deben pertenecer al
intervalo cerrado [1;2], solo vemos que se cumple para un punto
con esto se concluye que no existe punto fijo en esta función
para este intervalo.

11) Resolver la sgte función con punto fijo e indicar su E s <10-2

f(X)= X4 – 3X2-3, X0=1, g(X)=(3X2+3)1/4.
Solución:
g(X)= (3X2+3)1/4
X0=1
X1= g(X0) =1.5650846
X2= g(X1) =1.7935729
X3= g(X2) =1.8859438
 X4= g(X3) =1.9228479
X5= g(X4) =1.9375076
X6= g(X5) =1.9433170
X7= g(X6) =1.9456169
X8= g(X7) =1.9465271
X9= g(X8) =1.9468872

n Xn |Xn+1 – Xn|
0 1
1 1.5650846 0.5650846
2 1.7935729 0.2284883
3 1.8859438 0.0923709
4 1.9228479 0.0369041
5 1.9375076 0.0146597
6 1.9433170 0.0058094<0.01

La raíz con Es <10-2 es X6 =1.9433170.

X
12) Resolver por punto fijo la sgte f(X)= X −sec ⁡( 2 ) dando como
resultado EA % ≤ 0.1%.
X
Función de iteración : g(X) = sec ⁡( 2 ) , V0=0.5.

Solución:
X
g(X) = sec ⁡( 2 )

X0=0.5
X1= g(X0) =1.318116072
X2= g(X1) =0.851230713
X3= g(X2) =1.131154526
 X4= g(X3) =0.969663270
X5= g(X4) =1.064625673
X6= g(X5) =1.009466020
X7= g(X6) =1.041723665
X8= g(X7) =1.022936091
X9= g(X8) =1.033903988
X10= g(X9) =1.027509915
X11= g(X10) =1.031240514

n Xn |Xn+1-Xn| EA%
0 0.5
1 1.318116072 0.818116072 62.1%
2 0.851230713 0466885359 54.8%
3 1.131154526 0.279923813 24.7%
4 0.969663270 0.161491256 16.7%
5 1.064625673 0.094962403 8.9%
6 1.009466020 0.055159653 5.5%
7 1.041723665 0.032257645 4.0%
8 1.022936091 0.018787574 1.8%
9 1.033903988 0.010967897 1.1%
10 1.027509915 0.006394073 0.6%
11 1.031240514 0.003730599 0.3%
12 1.029064927 0.002175587 0.2%
13 1.030334017 0.001269990 0.12% ≤ 0.1%

La raíz cuyo EA%≤ 0.1% es X13=1.030334017.

13) Utilizar el método Newton Raphson con un E s < 0.01% y un

X0= 5 para la sgte función: f(X)= -0.9X2+1.7X+2.5.
Solución:
f(X)= -0.9X2+1.7X+2.5
f’(X)= -1.8X+1.7X
X0= 5
f ( X0)
X1= X 0− f ’( X )
0

(−11.50000)
X1=5− (−7.30000)

X1=3.42466

(−2.23354 )
X2=3.42466−
(−4.46439 )

X2=2.92436

(−0.22528 )
X3=2.92436−
(−3.56385 )

X3=2.86115

(−0.00361 )
X4=2.86115−
(−3.45007 )

X4=2.86010

( 0.00002 )
X5=2.86010−
(−3.44818 )

X5=2.86011
 (−0.00002 )
X6=2.86011−
(−3.44820 )

X6=2.86010
i Xi │Xi+1-Xi│%
0 5
1 3.42466 157%
2 2.92436 50%
3 2.86115 6%
4 2.86010 0.105%
5 2.86011 0.001% < 0.01%
6 2.86010 0.001%
7 2.86011 0.001%

La raíz cuyo error estimado es 0.001% < 0.01% es X5=2.86011.

14) Resolver por el método Newton Raphson Modificado la sgte

función: f(X)= tan(X)-0.1X , X0=π, con una tolerancia de 0.0001
(8 dígitos).
Solución:
f(X)= tan(X)-0.1X
f’(X)= sec2(X)-0.1
f’’(X)= 2sec2(X)*tg(X)
f ( X 0 )∗f ' (X 0 )
X 1 =X 0− 2
'
[ f ( X ) ] −f ( X )∗f ' ' ( X )
0 0 0

(−0.31415927 )∗(0.9)
X 1 =π − 2
[ 0.9 ] −(−0.31415927 )∗(0.00000000)

X 1 =¿3.49065851
 ( 0.01490439 )∗(1.03247434)
X 2 =3.49065851−
[ 1.03247434 ]2 −( 0.01490439 )∗(0.82437391)

X 2 =¿3.47605458

(−0.00008751 )∗(1.02076872)
X 3 =3.47605458−
[ 1.02076872 ]2−(−0.00008751)∗(0.77897449)

X 3 =¿3.47614030

(−0.00000001)∗(1.02083551)
X 4=3.47614030−
[ 1.02083551 ] 2−(−0.00000001 )∗(0.77923628)

X 4=¿ 3.47614031

n Xn │Xn+1-Xn│
0 π
1 3.49065851 0.34906586
2 3.47605458 0.01460393
3 3.47614030 0.00008572 < 0.0001
4 3.47614031 0.00000001

La raíz cuyo valor estimado es < 0.0001, es X 3=3.47614030.

15) Resolver la sgte función por el método de newton Rapshon

modificado f(X)=2 cos( X )−tg ¿), X0=0.5, Es < 0.0040 (9 dígitos).

Solución:
 ( 1.499823203 )∗(−1.491450826)
X 1 =0.5−
[ −1.491450826 ] 2−( 1.499823203 )∗(−2.027155210)

X 1 =¿0.924880824

( 0.705367616 )∗(−2.221345160)
X 2 =0.924880824−
[ −2.221345160 ] 2−( 0.705367616 )∗(−1.826223785)

X 2 =¿1.176685809

( 0.100815771)∗(−2.569227073)
X 3 =1.176685809−
[ −2.569227073 ]2 −( 0.100815771 )∗(−1.732084742)

X 3 =¿1.214914239

( 0.001694110 )∗(−2.616289211)
X 4=1.214914239−
[ −2.616289211 ]2−( 0.001694110 )∗(−1.727881201)

X 4=¿ 1.215561486

n Xn │Xn+1-Xn│
0 0.5
1 0.924880824 0.424880824
2 1.176685809 0.251804985
3 1.214914239 0.038228430
4 1.215561486 0.000647247<0.0040

La corriente deberá de ser X4= 1.215561486 amperes, el cual posee

un Es<0.0040.
 16) Resolver por el método de newton Raphson modificado la sgte
función: f(X)= X .3 −5 X 2 +7 X −3, X0= 0, Vv=1(4 dígitos).

Solución:
f(X)= X .3 −5 X 2 +7 X −3
f’(X)=3 X 2−10 X+ 7
f’’(X)=6 X −10

(−3 )∗(7.0000)
X 1 =0−
[ 7 ] 2−(−3 )∗(−10)

X 1 =¿1.1053

(−0.0210 )∗(−0.3879)
X 2 =1.1053−
[ −0.3879 ] 2−(−0.0210 )∗(−3.3682)

X 2 =¿1.0031

( 0.0000 )∗(−0.0124)
X 3 =1.0031−
[−0.0124 ] 2−( 0.0000 )∗(−3.9814)

X 3 =¿1.0031

n Xn Er%
0 0 100%
 1 1.1053 10.53%
2 1.0031 0.31%
3 1.0031 0.31%

El error relativo % es significativo, demuestra lo más cercano que se

encuentra a la verdadera raíz, entonces la raíz con menor E r% será
X2=1.0031más cercano al verdadero valor.

17) Resolver por el método de secante modificado la sgte función:

f(X)=e− X −X , ∂ = 0.01, X0= 1.0, Vv= 0.56714329(5 dígitos).

Solución:
f(X)=e− X −X
∂ = 0.01
X0= 1.0

( 0.01000 )∗(1.00000 )∗(−0.63212)

X 1 =1.0−
(−0.64578)−(−0.63212 )

X 1 =¿0.53725

( 0.01000 )∗( 0.53725 )∗(0.04710)

X 2 =0.53725−
(0.03860)− ( 0.04710 )

X 2 =¿0.56702
 ( 0.01000 )∗( 0.56702 )∗( 0.00019)
X 3 =0.56702−
(−0.00868)−( 0.00019 )

X 3 =¿0.56714

( 0.01000 )∗( 0.56714 )∗(0.00001)

X 4=0.56714−
(−0.00887)−( 0.00001 )

X 4=¿ 0.56715

( 0.01000 )∗( 0.56715 )∗(−0.00001)

X 5 =0.56715−
(−0.00889)−(−0.00001 )

X 5 =¿0.56714

n Xn ER%
0 1.0 76.32228%
1 0.53725 5.27085%
2 0.56702 0.02174%
3 0.56714 0.00058%
4 0.56715 0.00118%

Como se puede observar en el cuadro, hallamos el E R% ya que nos

dan el valor verdadero.

18) Indicar la temperatura en la que se enfría el metal por el

método de secante modificado la sgte función: f(X)=
X 3 −12 X+ 18 X−5 , ∂ = 0.01, X0= 1.5.
Solución:
f(X)= X 3 −12 X+ 18 X−5
∂ = 0.01
X0= 1.5

( 0.01000 )∗(1.5 )∗(−1.6250000)

X 1 =1.5−
(−1.7954341)−(−1.6250000 )

X 1 =¿1.3569828

( 0.01000 )∗( 1.3569828 )∗(−0.1723862)

X 2 =1.3569828−
(−0.2965606)−(−0.1723862)

X 2 =¿1.3381444

( 0.01000 )∗( 1.3381444 )∗(−0.0048439)

X 3 =1.3381444−
(−0.1232727)−(−0.0048439 )

X 3 =¿1.3375971

( 0.01000 )∗( 1.3375971 )∗(−0.0000610)

X 4=1.3375971−
(−0.1183247)−(−0.0000610 )

X 4=¿ 1.3375902

( 0.01000 )∗( 1.3375902 )∗(−0.0000007)

X 5 =1.3375902−
(−0.1182621)− (−0.0000007 )
 X 5 =¿1.3375901

( 0.01000 )∗( 1.3375901 )∗(0.0000002)

X 6 =1.3375901−
(−0.1182612)−( 0.0000002 )

X 6 =¿1.3375901

i Xi │Xi+1-Xi│
0 1.5
1 1.3569828 0.1430172
2 1.3381444 0.0188384
3 1.3375971 0.0005473
4 1.3375902 0.0000069
5 1.3375901 0.0000001
6 1.3375901 0.0000000
Temperatura a la que se enfría el metal es X 6=1.3375901ºC, cuyo error
estimado resulta cero.

19) Resolver por el método de la secante y realizar “n” iteraciones

hasta encontrar la raíz verdadera Xn= -2 en la sgte función: f(X)=
X 3 −3 X +2, X0= -2.6, X1=-2.4 (8 dígitos).

(−4.62400000 )∗(−2.4+ 2.6 )

X 2 =−2.4−
(−4.62400000 )−(−7.77600000 )

X 2 =−2.10659898

(−1.02878220 )∗ (−2.10659898+ 2.4 )

X 3 =−2.10659898−
(−1.02878220 )− (−4.62400000 )
X 3 =−2.02264141

(−0.20686010 )∗(−2.02264141+2.10659898 )
X 4=−2.02264141−
(−0.20686010 ) + ( 1.02878220 )

X 4=−2.00151110

(−0.01361360 )∗ (−2.00151110+2.02264141 )
X 5 =−2.00151110−
(−0.01361360 ) −(−0.20686010 )

X 5 =−2.00002254

(−0.00020286 )∗(−2.00002254 +2.00151110 )

X 6 =−2.00002254−
(−0.00020286 )−(−0.01361360 )

X 6 =−2.00000002

(−0.00000018 )∗(−2.00000002+2.00002254 )
X 7 =−2.00000002−
(−0.00000018 )− (−0.00020286 )

X 7 =−2.00000000

En la séptima iteración encontramos el valor verdadero de la raíz.ç

20) Resolver mediante el método de la secante la sgte función:

f(X)= X 3 +2 X 2 +10 X −20=0, X0=0, X1=1.0, ES= 10-3. (6 DÍGITOS)
 (−7.000000 )∗( 1−0 )
X 2 =1.0−
(−7.000000 )−(−20.000000 )

X 2 =1.538462

( 3.759683 )∗( 1.538462−1.0 )

X 3 =1.538462−
( 3.759683 )−(−7.000000 )

X 3 =1.350311

(−0.388135 )∗( 1.350311−1.538462 )

X 4=1.350311−
(−0.388135 )−( 3.759683 )

X 4=1.367917

(−0.018794 )∗( 1.367917−1.350311 )

X 5 =1.367917−
(−0.018794 ) −(−0.388135 )

X 5 =1.368813

( 0.000103 )∗ (1.368813−1.367917 )
X 6 =1.368813−
( 0.000103 )− (−0.018794 )

X 6 =1.368808

i Xi │Xi+1-Xi│
0 0
1 1 1
2 1.538462 0.538462
3 1.350311 0.188151
4 1.367917 0.017606
5 1.368813 0.000896<0.001
6 1.368808 0.000005

La raíz con error estimado< 0.001 es X5=1.368813.

 21) Resolver la sgte función utilizando el método secante
modificada: f(X)= sen(X)+cos(1+X2)-1, ∂=0.02, X0=3; localize la
primera raíz positiva. (6 dígitos)

∂∗X 0∗f ( X 0 )
X 1 =X 0−
f ( X +∂∗X 0 ) −f ( X 0 )

( 0.02 )∗(3)∗(−1.697952)
X 1 =3−
(−1.509237 )−(−1.697952)

X 1 =3.539846

i Xi │Xi+1-Xi│
0 3
1 3.539846 0.539846

La primera raíz positiva resulto en X1=3.539846.

22) La resistencia de la tensión de un plástico aumenta como una

función del tiempo cuando se calienta. Se dispone de los sgtes
datos:

Puntos 0 1 2 3 4
Tiempo 10 15 20 25 30
Resistencia 4 20 18 50 33

¿Cuál será la resistencia del plástico a cierta temperatura?

Elabore la tabla de diferencias divididas:

Solución:
 i Xi f(Xi) Primera Segunda Tercera Cuarta
0 10 4 3.20000000 -0.36000000
1 15 20
-0.40000000 0.06933333
2 20 18 0.68000000
6.40000000 -0.00900000
- 25 50 -0.98000000
3 -0.11066667
-3.40000000
4 30 33

Primera diferencia:
f ( X 1 )−f (X 0) 20−4
f (X0, X1) = = =3.20000000
X 1−X 0 15−10
f ( X 2 )−f (X 1) 18−20
f (X1, X2) = X 2−X 1
=
20−15
=−0.40000000

f ( X 3 )−f ( X 2) 50−18
f (X2, X3) = X 3−X 2
=
25−20
=6.40000000

f ( X 4 )−f ( X 3) 33−50
f (X3, X4) = X 4−X 3
=
30−25
=−3.40000000

Segunda diferencia:
f ( X 1 , X 2 ) −f ( X 0 , X 1) −0.4−3.2
f (X0, X1, X2) = = =−0.36000000
X 2− X 0 20−10

f ( X 2 , X 3 )−f (X 1 , X 2) 6.4+0.4
f (X1, X2, X3) = = =0.68000000
X 3− X 1 25−15

f ( X 5 , X 4 )−f ( X 2 , X 3) −3.4−6.4
f (X2, X3, X4) = = =−0.98000000
X 4−X 2 30−20
Tercera diferencia:
f ( X 1 , X 2 , X 3 )−f ( X 0 , X 1 , X 2) 0.68+0.36
f (X0, X1, X2, X3) = = =0.06933333
X 3−X 0 25−10

f ( X 2 , X 3 , X 4 )−f ( X 1 , X 2 , X 3) −0.98−0.68
f(X1,X2,X3,X4)= = =−0.11066667
X 4−X 1 30−15

Cuarta diferencia:
f(X0,X1,X2,X3,X4)=
f ( X 1 , X 2 , X 3 , X 4 )−f (X 0 , X 1, X 2 , X 3) −0.11066667−0.06933333
= =−0.00900000
X 4−X 0 30−10

23) La producción de arroz durante el periodo 2012 al 2020 se

refleja en la sgte tabla:

AÑOS 2012 2014 2016 2018 2020

PRODUCCIÓN 1229 1302 1100 1200 1220

¿Cuál fue la producción de arroz para el año 2017, con el problema

que se tuvo del niño costero?
Solución:
i Xi f(Xi) Primera Segunda Tercera Cuarta
0 201 1229
2 36.5 -34.375000
1 201 1302 12.02083
4 3
2 201 1100 -101.0 -2.497396
6 37.750000
3 201 1200 50.0
8
 4 202 1220 -10.000000 -7.958333
0 10.0

Primera diferencia:
f ( X 1 )−f (X 0) 1302−1229
f (X0, X1) = X 1−X 0
=
2014−2012
=36.5

f ( X 2 )−f (X 1) 1100−1302
f (X1, X2) = X 2−X 1
=
2016−2014
=−101.0

f ( X 3 )−f ( X 2) 1200−1100
f (X2, X3) = X 3−X 2
=
2018−2016
=50

f ( X 4 )−f ( X 3) 1220−1200
f (X3, X4) = X 4−X 3
=
2020−2018
=10.0

Segunda diferencia:
f ( X 1 , X 2 ) −f ( X 0 , X 1) −101−36.5
f (X0, X1, X2) = X 2− X 0
=
2016−2012
=−34.375000

f ( X 2 , X 3 )−f (X 1 , X 2) 50+101
f (X1, X2, X3) = X 3− X 1
=
2018−2014
=37.750000

f ( X 3 , X 4 )−f ( X 2 , X 3) 10−50
f (X2, X3, X4) = X 4−X 2
=
2020−2016
=−10

Tercera diferencia:
f (X0, X1, X2, X3) =
f ( X 1 , X 2 , X 3 )−f ( X 0 , X 1 , X 2) 37.750000+34.375000
= =12.020833
X 3−X 0 2018−2012
 f ( X 2 , X 3 , X 4 )−f ( X 1 , X 2 , X 3) −10−37.750000
f(X1,X2,X3,X4)= = =−7.958333
X 4−X 1 2020−2014

Cuarta diferencia:

f(X0,X1,X2,X3,X4)=
f ( X 1 , X 2 , X 3 , X 4 )−f (X 0 , X 1, X 2 , X 3) −7.958333−12.020833
= =−2.497396
X 4−X 0 2020−2012

24) La tabla presenta la temperatura de un químico en diferentes

presiones.

PUNTO 0 1 2 3 4 5
S
XI 1 2 2.5 3 4 5
F(XI) 1 5 7 8 2 1

Elaborar una aproximación polinomial de Newton para la

información tabular dada.
Interpole la temperatura para una presión de 3.4, mediante
polinomios de interpolación de orden 1, 2 y 3.
Para esto realizamos la tabla de diferencias divididas para tener
los coeficientes del polinomio de Newton.
Solución:
i Xi f(Xi) Primera Segunda Tercera Cuarta Quinta
0 1 1
1 2 5 4.0 0.000000
-1.000000
-
2 2.5 7 -2.000000 0.222222
4.0

2.0 -1.666667
 3 3 8 0.455556

-6.0 -5.333333
4 4 2 1.600000
3.133333
2.500000
-1.0
5 5 1

Primera diferencia:
f ( X 1 )−f (X 0)
f (X0, X1) = X 1−X 0
=4.0

f ( X 2 )−f (X 1)
f (X1, X2) = X 2−X 1
=4.0

f ( X 3 )−f ( X 2)
f (X2, X3) = X 3−X 2
=2.0

f ( X 4 )−f ( X 3)
f (X3, X4) = X 4−X 3
=−6.0

f ( X 5 )−f ( X 4)
f (X4, X5) = X 5−X 4
=−1.0

Segunda diferencia:
f ( X 1 , X 2 ) −f ( X 0 , X 1) 4−4
f (X0, X1, X2) = X 2− X 0
=
2.5−1
=0.000000

f ( X 2 , X 3 )−f (X 1 , X 2) 2−4
f (X1, X2, X3) = X 3− X 1
=
3−2
=−2.000000

f ( X 3 , X 4 )−f ( X 2 , X 3) −6−2
f (X2, X3, X4) = X 4−X 2
=
4−2.5
=−5.333333
 f ( X 4 , X 5 )−f ( X 3 , X 4) −1+6
f (X3, X4, X5) = X 5− X 3
=
5−3
=2.500000

Tercera diferencia:
f ( X 1 , X 2 , X 3 )−f ( X 0 , X 1 , X 2) −2−0
f (X0, X1, X2, X3) = X 3−X 0
=
3−1
=−1.000000

f ( X 2 , X 3 , X 4 )−f ( X 1 , X 2 , X 3) −5.333333+2
f(X1,X2,X3,X4)= X 4−X 1
=
4−2
=−1.666667

f ( X 3 , X 4 , X 5 ) −f ( X 2, X 3 , X 4) 2.5+5.333333
f(X2,X3,X4,X5)= X 5− X 2
=
5−2.5
=3.133333

Cuarta diferencia:
f(X0,X1,X2,X3,X4)=
f ( X 1 , X 2 , X 3 , X 4 )−f (X 0 , X 1, X 2 , X 3) −1.666667+1
= =−0.222222
X 4−X 0 4−1

para n=1
P1(X)=a0+(X-X0)*a1
P1(X)=f(X0)+ (X-X0)*f(X0, X1)
P1(X)=1+4(X-1)
P1(X)=4X-3
P1(3.4)= 4*3.4-3
P1(3.4)= 10.600000

para n=2
P2(X)=a0+(X-X0)*a1+(X-X0)*(X-X1)*a2
P2(X)=f(X0)+ (X-X0)*f(X0, X1) +(X-X0)*(X-X1)*f(X0,X1,X2)
P2(X)=1+4(X-1) + (X-X0)*(X-X1)*0
P2(X)=4X-3
P2(3.4) = 4*3.4-3
P2(3.4) = 10.600000

Para n=3
P3(X)=a0+(X-X0)*a1+(X-X0)*(X-X1)*a2+(X-X0)*(X-X1)(X-X2)*a3
P3(X)=f(X0)+ (X-X0)*f(X0, X1) +(X-X0)*(X-X1)*f(X0,X1,X2)+(X-X0)(X-X1)(X-X2)f(X0,X1,X2,X3)
P3(X)=1+4(X-1) + (X-1)(X-2)(X-2.5)*-1
P3(X)=1-4X-4-(X-1)(X-2)(X-2.5)
P3(3.4) = 1+4(3.4)-4-(3.4-1)*(3.4-2)*(3.4-2.5)
P3(3.4) = 7.576000

25) Los siguientes datos representa la producción de oro refinado

en mina Pierina:
PUNTOS i 0 1 2 3 4
Año Xi 2000 2002 2004 2006 2008
Au f(Xi) 105 107 110 108 122
Elaborar una aproximación polinomial de Newton para la
información tabular dada. Interpole la producción de oro en el
año 2005, mediante polinomios de interpolación de orden 1,2 y
3.

Solución:

i Xi f(Xi) Primera Segunda Tercera Cuarta

0 200 105 1.000 0.125 -0.125 0.070
0
 1 200 107 1.500
2
2 200 110 -0.625
4 -1.000 0.438
3 200 108 2.0
6 7.000
4 200 122
8

Primera diferencia:
f ( X 1 )−f (X 0)
f (X0, X1) = =1 .00 0
X 1−X 0
f ( X 2 )−f (X 1)
f (X1, X2) = X 2−X 1 =1.500
f ( X 3 )−f ( X 2)
f (X2, X3) = X 3−X 2 =−1.000
f ( X 4 )−f ( X 3)
f (X3, X4) = X 4−X 3 =7.000

Segunda diferencia:

f ( X 1 , X 2 ) −f ( X 0 , X 1) 1.500−1.000
f (X0, X1, X2) = = =0. 125
X 2− X 0 4

f ( X 2 , X 3 )−f (X 1 , X 2) −1.000−1.500
f (X1, X2, X3) = = =−0 . 625
X 3− X 1 4

f ( X 3 , X 4 )−f ( X 2 , X 3) 7+1
f (X2, X3, X4) = = =2.000
X 4−X 2 4

Tercera diferencia:
f (X0, X1, X2, X3) =
f ( X 1 , X 2 , X 3 )−f ( X 0 , X 1 , X 2) −0.625−0.125
= =−0 . 125
X 3−X 0 6

f ( X 2 , X 3 , X 4 )−f ( X 1 , X 2 , X 3 ) 2.000+0.625
f(X1,X2,X3,X4)= = =0.4375
X 4− X 1 6

Cuarta diferencia:
 f(X0,X1,X2,X3,X4)=
f ( X 1 , X 2 , X 3 , X 4 )−f (X 0 , X 1, X 2 , X 3) 0.4375+0. 1 25
= =0.070
X 4−X 0 8

para n=1
P1(X)=f(X0)+ (X-X0)*f(X0, X1)
P1(X)=105+(X-2000)
P1(2005)=105+(2005-2000)
P1(2005)= 110

para n=2
P2(X)=f(X0)+ (X-X0)*f(X0, X1) +(X-X0)*(X-X1)*f(X0,X1,X2)
P2(X)=105+(X-X0)*1 + (X-X0)*(X-X1)*0.125
P2(2005)=105+(2005-2000)+(2005-2000)*(2005-2002)*0.125
P2(2005) = 111.875

Para n=3
P3(X)=f(X0)+(X-X0)*f(X0,X1)+(X-X0)*(X-X1)*f(X0,X1,X2)+(X-X0)(X-X1)(X-
X2)f(X0,X1,X2,X3)
P3(2005)=105+(2005-2000)*1+ (2005-2000)(2005-2002)*0.125+(2005-2000)(2005-
2002)(2005-2004)*-0.125
P3(2005)=110

26) Se sabe que la resistencia a la tensión de un plástico aumenta

como una función del tiempo cuando se calienta. Se dispone de
los siguientes datos:
Puntos 0 1 2 3 4
 Tiempo 10 15 20 25 30
Resistencia 4 14 18 50 33
a la tensión

Aproxime la función f(X) tabulada por el polinomio de Newton en

diferencias divididas hacia adelante e interpole cuando el tiempo es
20.
Solución:
i Xi f(X) Δ f(X) Δ2 f(X) Δ3 f(X) Δ4 f(X)
0 10 4 10
1 15 14 -6
4 34
2 20 18 28
32 -111
3 25 50 -77
-17 -49
4 30 33

Δ f(X):
Δ f(X0)= f(X1)-f(X0)=14-4=10
Δ f(X1)= f(X2)-f(X1)=18-14=4
Δ f(X2)= f(X3)-f(X2)=50-18=32
Δ f(X3)= f(X4)-f(X3)=33-50=-17
Δ2 f(X):
Δ2 f(X0)= Δ f(X1)- Δ f(X0)=4-10=-6
Δ2 f(X1)= Δ f(X2)- Δ f(X1)=32-4=28
Δ2 f(X2)= Δ f(X3)- Δ f(X2)=-17-32=-49

Δ3 f(X):
Δ3 f(X0)= Δ2 f(X1)- Δ2 f(X0)=28+6=34
Δ3 f(X1)= Δ2 f(X2)- Δ2 f(X1)=-49-28=-77

n=1 Primer grado

P1(X)=f(X1)+ Δf(X1)*s
P1(20)=14+4*1
P1(20)=14+4=18

n=2 Segundo grado

∆2 f ( X 1 )∗s∗( s−1)
P2(X)=f(X1)+Δf(X1)*s+
2!
28∗1∗0
P2(20)=14+4*1+ 2!

P2(20)=18

27) La aleación es una mezcla solida homogénea de dos o más

metales, o de uno o más metales con algunos elementos no
metálicos.
A continuación se muestra una tabla con la presión de una
aleación a distintas temperaturas. Aproxime la función tabulada
de Newton en diferencias divididas hacia adelante.

Puntos 0 1 2 3 4
Temperatura 10 20 30 40 50
Presión 25.24 37.81 41.25 53.14 61.07

Solución:
i Xi f(X) Δf(X) Δ2 f(X) Δ3 f(X) Δ4 f(X)

0 10 25.2
4 12.57 -9.13
17.58
1 20 37.8
1 3.44
-29.99
2 30 41.2 8.45
5 11.89
-12.41
3 40 53.1
4 7.93 -3.96
4 50 61.0
7

Δ f(X):
Δ f(X0)= f(X1)-f(X0)=37.81-25.24=12.57
Δ f(X1)= f(X2)-f(X1)=41.25-37.81=3.44
Δ f(X2)= f(X3)-f(X2)=53.14-41.25=11.89
Δ f(X3)= f(X4)-f(X3)=61.07-53.14=7.93
Δ2 f(X):
Δ2 f(X0)= Δ f(X1)- Δ f(X0)=3.44-12.57=-9.13
Δ2 f(X1)= Δ f(X2)- Δ f(X1)=11.89-3.44=8.45
Δ2 f(X2)= Δ f(X3)- Δ f(X2)=7.93-11.89=-3.96

Δ3 f(X):
Δ3 f(X0)= Δ2 f(X1)- Δ2 f(X0)=8.45-9.13=17.58
Δ3 f(X1)= Δ2 f(X2)- Δ2 f(X1)=-3.96-8.45=-12.41

Δ4 f(X):
Δ4 f(X0)= Δ3 f(X1)- Δ3 f(X0)=-12.41-17.58=-29.99

De aquí se puede obtener la aproximación para cualquier valor

de “X” , ya sea de grado 1, 2, 3 o 4.
 28)
a) Construya los polinomios de grado 1, 2 y 3 con los siguientes
datos: (8 dígitos )

f(-0.75)=-0.07181250, f(-0.5)=-0.02475000, f(-0.25)=0.33493750,

f(0)=1.10100000, empleando el método de diferencia regresiva
de Newton.
b) Interpole para un valor de X=-1/3.

Solución:
Puntos 0 1 2 3
Xi -0.75 -0.50 -0.25 0
f(Xi) -0.07181250 -0.02475000 0.33493750 1.10100000

i Xi f(Xi) ∇f(Xi) ∇ 2f(Xi) ∇ 3f(Xi)

0 - -0.07181250
0.7
0.04706250 0.31262500
5
1 -0.5 -0.02475000
0.3596875
0 0.09375000
2 - 0.33493750 0.40637500
0.2
0.76606250
5
3 0 1.10100000

∇ f ( X i ):

∇ f ( X 1 ) =f ( X 1 )−f ( X 0 )=−0.02475000−(−0.07181250)=0.04706250

∇ f ( X 2 ) =f ( X 2 )−f ( X 1 )=0.33493750−(−0.02475000)=0.35968750
∇ f ( X 3 ) =f ( X 3 ) −f ( X 2 )=1.10100000−0.33493750=0.76606250

∇2 f ( X i ) :

∇ 2 f ( X 2 ) =∇ ( ∇ f ( X 2 ) ) =∇ f ( X 2 )−∇ f ( X 1 )=0.35968750−0.04706250=0.31262500

∇ 2 f ( X 3 ) =∇ ( ∇ f ( X 3 ) ) =∇ f ( X 3 ) −∇ f ( X 2 )=0.76606250−0.3596875 0=0. 40637500

∇3 f ( X i ) :

∇ 3 f ( X 3 ) =∇ ( ∇ 2 f ( X 3 ) )=∇ 2 f ( X 3 ) −∇ 2 f ( X 2 )=0. 40637500−0. 31262500=0. 09375 00 0

Donde:
−1
−0
3
s= =−1.33333333
0.25

*Para n=1
P1(-1/3) = f ( X 3 ) +∇ f ( X 3 )∗s
P1(-1/3) = 1.10100000+0.76606250*(-1.33333333)
P1(-1/3) = 0.07958334

*Para n=2
s∗( s+1 )∗∇ 2 f (X 3 )
P2(-1/3) = f ( X 3 ) +∇ f ( X 3 )∗s+
2!
P2(-1/3) =
(−1.33333333 )∗(−0.33333333 )∗(0.40637500)
1.10100000+0.76606250∗(−1.33333333)+
2!
P2(-1/3) = 0.16988889

*Para n=3
 s∗( s+1 )∗∇ 2 f ( X 3 ) s∗( s +1 )∗( s+2 )∗∇ 3 f ( X 3 )
P3(-1/3) = f ( X 3 ) +∇ f ( X 3 )∗s+ +
2! 3!
P3(-1/3) =
(−1.33333333 )∗(−0.33333333 )∗( 0.40637500 ) 0.09375∗−1.333
1.10100000+0.76606250∗(−1.33333333 ) + +
2!
P3(-1/3) = 0.17451852

29) Interpole el valor de la presión a una temperatura de 98ºF,

empleando la tabla anterior y el polinomio de Newton en
diferencias divididas hacia atrás. (6 dígitos)

Solución:
i Xi f(Xi) ∇ f(Xi) ∇ 2f(Xi) ∇ 3f(Xi) ∇ 4 f(Xi) ∇ 5f(Xi)
0 50 24.94
60 30.11 5.170000 0.770000
1
2.010000
70 36.05 5.940000 0.350000
-0.420000
2

6.290000
-5.040000
80 42.34 1.940000 1.590000
3
8.230000
50.57 -3.030000
4 90

8.730000 0.500000 -1.440000

59.30
5 100

∇ f ( X i ):
∇ f ( X 1 ) =f ( X 1 )−f ( X 0 )=30.11−24.94=5.170000

∇ f ( X 2 ) =f ( X 2 )−f ( X 1 )=36.05−30.11=5 . 9 40000

∇ f ( X 3 ) =f ( X 3 ) −f ( X 2 )=42.34−36.05=6.290000

∇ f ( X 4 )=f ( X 4 )−f ( X 3 )=50.57−42.34=8.230000

∇ f ( X 5 ) =f ( X 5 ) −f ( X 4 )=59.30−50.57=8.73000 0

∇2 f ( X i ) :

∇ 2 f ( X 2 ) =∇ ( ∇ f ( X 2 ) ) =∇ f ( X 2 )−∇ f ( X 1 )=5.940000−5.170000=0.770000

∇ 2 f ( X 3 ) =∇ ( ∇ f ( X 3 ) ) =∇ f ( X 3 ) −∇ f ( X 2 )=6.290000−5.940000=0.350000

∇ 2 f ( X 4 )=∇ ( ∇ f ( X 4 ) ) =∇ f ( X 4 )−∇ f ( X 3 ) =8.230000−6.290000=1.940000

∇ 2 f ( X 5 ) =∇ ( ∇ f ( X 5 ) ) =∇ f ( X 5 ) −∇ f ( X 4 ) =8.730000−8.230000=0.500000

∇3 f ( X i ) :

∇ 3 f ( X 3 ) =∇ ( ∇ 2 f ( X 3 ) )=∇ 2 f ( X 3 ) −∇ 2 f ( X 2 )=0. 350000−0.770000=−0.420000

∇ 3 f ( X 4 )=∇ ( ∇2 f ( X 3 ) )=∇2 f ( X 4 ) −∇ 2 f ( X 3 )=1.940000−0.350000=1.590000

∇ 3 f ( X 5 ) =∇ ( ∇ 2 f ( X 3 ) )=∇ 2 f ( X 5 ) −∇ 2 f ( X 4 ) =0. 500000−1.940000=−1.440000

∇4 f ( X i ) :

∇ 4 f ( X 4 ) =∇ ( ∇ 3 f ( X 4 ) )=∇3 f ( X 4 ) −∇ 3 f ( X 3 )=1.590000−(−0.420000)=2.010000

∇ 4 f ( X 5 )=∇ ( ∇3 f ( X 5) )=∇3 f ( X 5 )−∇3 f ( X 4 )=−1.440000−1.590000=−3.0 30000

∇5 f ( X i ) :

∇ 5 f ( X 5 ) =∇ ( ∇ 4 f ( X 5) )=∇2 f ( X 5 )−∇2 f ( X 4 )=−3.030000−2.010000=−5.040000

Donde:
 98−10 0
s= =−0.200000
10

*Para n=1
P1(98) = f ( X 5 ) +∇ f ( X 5 )∗s
P1(98) = 59.300000+(8.73*-0.200000)
P1(98) = 57.554000

*Para n=2
s∗( s+1 )∗∇ 2 f (X 5 )
P2(98) = f ( X 5 ) +∇ f ( X 5 )∗s+
2!
( 0.500000 )∗(−0.200000 )∗( 0.800000)
P2(98) = 59.300000+ ( 8.730000∗−0.200000 ) +
2!
P2(98) = 57.514000

*Para n=3
s∗( s+1 )∗∇ 2 f ( X 5 ) s∗( s +1 )∗(s +2)∗∇ 3 f ( X 5)
P3(98) = f ( X 5 ) +∇ f ( X 5 )∗s+ +
2! 3!
P3(98)=
( 0.500000 )∗(−0.200000 )∗( 0. 800000 ) −1.440000∗−0.200000∗0.8000
59.300000+ ( 8.730000∗−0.200000 ) + +
2! 3!
P3(98) =57.583120

30) Usar el método trapezoidal para integrar numéricamente:

f ( X )=3+5 X + 9 X 2 +7 X 3 , donde a=0 y b=3. (3 dígitos)

Analíticamente:
 3
I =∫ ( 3+5 X + 9 X 2 +7 X 3 ) dx
0

5 X2 9 X3 7 X4
I =3 X +
2
+
3
+
4
0→ 3 ]
I =254.250−0.000=254.250( valor verdadero)

Aplicando el método trapezoidal:

f(0)=3.000
f(3)=288.000

I =( b−a )∗ ( f ( a) +f2 ( b ) )
288+3
I =( 3−0 )∗( )
2
I =436.500

 El cual tiene un error de:

EA=│valor verdadero – valor aproximado│
EA=│254.250 – 463.500│=182.250
182.250
ER%=│ 254.250 │∗100 %

ER%=71.68124%
 Hallamos el ET:
f ( X )=3+5 X + 9 X 2 +7 X 3

f ' (X )=5+18 X + 21 X 2

f ' ' ( X )=18+ 42 X

El valor promedio de la segunda derivada, podemos calcularla por:

 b

~
∫ f ' ' ( X)dx
f ' ' ( x )=MEDIA= a
b−a
3

~
∫ (18+ 42 X) dx
f ' ' ( x )= 0
3−0

42 X 2
18 X +
~ 2
f ' ' ( x )=
3
~
f ' ' ( x )=6 ( 3 ) +6 ¿
~
f ' ' ( x )=72.00000 0

Donde:
−1 3
ET= 12 ∗( 72.000000 )∗( 3−0 ) =│−162│

ET =162.000000

31) Usar el método trapezoidal para integrar numéricamente:

f ( X )=1+2 X +3 X 2, donde a= -2 y b=4. (6 dígitos)

Analíticamente:
4
I =∫ ( 1+2 X +3 X 2 ) dx
−2

2 X2 3 X3
I =X +
2
+
3 ]
−2 → 4

I =90.000000(valor verdadero)

Aplicando el método trapezoidal:

f(-2)=9.000000
f(4)=57.000000

I =( b−a )∗ ( f ( a) +f2 ( b ) )
 9+57
I =( 4−(−2) )∗( )
2

I =198.00000 0

El error absoluto:
EA=│valor verdadero – valor aproximado│
EA=│90.000000-198.000000│=108.000000
108.000000
ER%=│ │∗100 %
90

ER%=120%
Hallamos el ET:
f ( X )=( 1+2 X +3 X 2 )

f ' ( X ) =2+6 X

f ' ' ( X )=6

El valor promedio de la segunda derivada, podemos calcularla por:

b

~
∫ f ' ' ( X)dx
f ' ' ( x )=MEDIA= a
b−a
4

~
∫ (6) dx
''
f ( x )= −2
4−(−2)
~ 6∗(4−2)
f ' ' ( x )=
6
~
f ' ' ( x )=2.000000

Donde:
 −1 3
ET= 12 ∗( 2.000000 )∗( 6 ) =│−3 6 .000000│

ET =3 6.00000 0