“AÑO DE LA LUCHA CONTRA LA CORRUPCIÓN E IMPUNIDAD"

UNIVERSIDAD NACIONAL DE
TRUJILLO

FACULTAD DE INGENIERIA

ESCUELA DE INGENIERIA DE SISTEMAS

CURSO :

ANÁLISIS NUMÉRICO

DOCENTE :

Mg .POÉMAPE ROJAS,
GLORIA

TITULO :

PRACTICA N°1

CICLO :

VI

INTEGRANTES :

AYALA HERNÁNDEZ AZUCENA

BARDALES TANTA ESTEFANY
SORET LEIVA, JULISSA
TERÁN MEJIA LESLY

GUADALUPE -2019
2
INDICE

Ejercicio 1 ………………………………………………………………………………………………..4

El feminicidio

Ejercicio 2 ………………………………………………………………………………………………..4

Ejercicio 3 ………………………………………………………………………………………………..5

Ejercicio 4 ………………………………………………………………………………………………..5

Ejercicio 5 ………………………………………………………………………………………………..6

Ejercicio 6 ………………………………………………………………………………………………..6

Ejercicio 7 ………………………………………………………………………………………………..9

Ejercicio 8 ………………………………………………………………………………………………..9

Ejercicio 9 ………………………………………………………………………………………………..14

Ejercicio 10 ………………………………………………………………………………………………15

Ejercicio 11……………………………………………………………………………………………….17

Ejercicio 12 ………………………………………………………………………………………………18

Ejercicio 13 ………………………………………………………………………………………………23

Ejercicio 14 ………………………………………………………………………………………………24

Ejercicio 15 ………………………………………………………………………………………………24

Ejercicio 16 ………………………………………………………………………………………………27

Ejercicio 17 ………………………………………………………………………………………………27

Ejercicio 18 ………………………………………………………………………………………………35

3
EJERCICIO 1: Sea f ( x )=x∗exp ⁡( x ); se sabe que x*=1.15 y que el error de x es 0.01.
Estime cuál sería el error al evaluar la función.

 f ( x )=x∗exp ⁡( x )

df ( x )=exp ( x ) +exp ( x )∗x

df ( x )=exp (1.15) +exp( 1.15)∗(1.15)=6.7901

 E f ¿=|0.01|×|6.7901|=0.0679

 ( 6.7901 ±0.0679 )=(6.7222 ; 6.858)

EJERCICIO 2: Se tiene un rectángulo cuyo largo mide 35 ± 0.2 m y su ancho 24 ± 0.3 m

a = Largo

b = Ancho

a) El Perímetro máximo y mínimo

 P=2 a+ 2b

Pmáx =2 ( 35.2 )+2 ( 24.3 ) =70.4+48.6

Pmáx =119m

 Pmín =2 ( 34.8 ) +2 ( 23.7 )=69.6+ 47.4

Pmín =117 m

b) El perímetro promedio

Pmáx + P mín 119+117

P prom= = =118 m
2 2

c) El intervalo del perímetro

Perímetro ∈(117 ; 119)

d) El error absoluto del perímetro

EA perímetro =¿ Pmín −P prom ∨¿

EA perímetro =¿ 117 m−118 m∨¿

EA per =1 m

4
EJERCICIO 3: En el ejercicio anterior, calcular:

a) El área máxima y mínima

 A=a∗b

Amáx =( 35.2 ) ( 24.3 )

Amáx =855.36 m 2

Amín =( 34.8 )( 23.7 )

Amín =824.76 m 2

b) El área promedio

A máx + Amín 855.36+824.76 2

A prom = = =840.06 m
2 2

c) El intervalo del área

Área∈(824.76 ; 855.36)

d) El error relativo del área

ER Área =¿ A mín − A prom ∨ ¿ ¿

A prom

824.76 m2−840.06 m2
ER Área =
840.06 m2

ER Área =0.0182 m2

EJERCICIO 4: En un experimento se introducen dos líquidos en un matraz y se quiere

hallar la masa total del líquido.
Se conocen:
M1 = Masa del matraz 1 + contenido = 530 ± 10 g
M1 = Masa del matraz 1 = 75 ± 2 g
M2 = Masa del matraz 2 + contenido = 860 ± 20 g
M2 = Masa del matraz 2 = 89 ± 1 g

Solución:

 La Masa de líquido será:

M = M1 − m1 + M2 − m2 = 1226g

5
  Su error:
∆ M = ∆M + ∆m + ∆M + ∆m = 33g

 El resultado se expresará:
M = 1226 ± 33 g

EJERCICIO 5: Para medir la altura de un árbol, L, se mide la longitud de su sombra, L1,

la altura de un objeto de referencia, L2, y la longitud de su sombra, L3. Realizadas las
medidas resultan:

L3=200 ± 2 cm, L2=100. 0± 0.4 cm, L3=10.3 ± 0.2cm .

Hallar la altura del árbol con su error.

SOLUCIÖN:

L2
Por semejanza: L=L1
L3

100
Por lo tanto L=200 =2000 cm
10

EL EL1 EL2 EL3

Su error será ≈ + +
|L| |L1| |L2| |L3|

EL 2 0.4 0.2
= + + =0.034
|L| 200 100 10.3

EL=0.034|L|=0.034 (2000)

EL=68

Finalmente L=2000 ± 68 cm.

EJERCICIO 6: Dos pequeños objetos que se encuentran en un suelo horizontal

separados por una distancia de 40m son lanzados en el mismo instante
verticalmente hacia arriba con velocidades de 120 y 180m/s. Calcular la distancia
entre los dos objetos cuando el más rápido se encuentra a la máxima altura que
puede alcanzar. Si todos los datos se conocen con una precisión del 1%.Calcular el
error absoluto obtenido. Para facilitar los cálculos, tome g=10m/s^2. Se supone
ausencia de rozamiento.

6
Resolución:

Objeto A: v 0 A Tiempo invertido por B en alcanzar la altura máxima:

cuando

Objeto B: v 0 B (¿ v 0 A ) llega a la altura máxima su velocidad es nula

v0 B
v 0 B−g t B =0 tB=
g

Y xB

ŕ BA (t)

ŕ A (t)

ŕ B (t)

A B X

Tomamos la posicion inicial del objeto A como origen de coordenadas.

Posicion de ambos objetos en cualquier instante

1
r⃗ A ( t )=(0 , y A ) y A ( t )=v 0 A t− g t 2
2

1
r⃗ B ( t )=( x B , y B ) y B ( t )=v 0 B t− g t 2
2

r⃗ AB ( t )=⃗r B ( t )−⃗r A ( t )=( x B , y B− y A )

Distancia AB en cualquier instante:

|r⃗ AB ( t )|=√ x B2 +( y B− y A )2 =+ √ x B2+(v 0 B−v 0 A )2 t 2

Cuanto t=t B

1
y A ( t B ) =v 0 A t B − g t 2B
2

1
y B ( t B )=v 0 B t B− g t 2B
2

7
|r⃗ BA ( t B )|=+ √ x B2+( v 0 B−v 0 A )2 t 2B
Porcentaje error: 1 7.0 %

v0A ( ms )=¿ 120 1.2

|r⃗ BA ( t B )|=(180 ± 80)m

m |r⃗ BA ( t B )|=+ √ 402 +(180−120)2 182

v0 B ( )
s
=¿ 180 1.8
|r⃗ BA ( t B )|=1081m
x B ( m )=¿ 40 0.4

|r⃗ BA ( t B )|=(180 ± 80)m

m
g
( )
s2
=¿ 10 0.1

t B ( s )=¿ 18 0.36
∆
r AB ( m )=¿ 1081 76

Calculo de Errores
Y xB

ŕ BA (t)

ŕ A (t)

ŕ B (t)

A B X
Variable auxiliar
u=x2B +( v 0 B −v 0 A )2 t 2B

∂u ∂u ∂u ∂u
∆ u=
| | | | | | | |
∂ xB
∆ xB +
∂ v0 B
∆ v0 B +
∂ v0 A
∆ v0 A +
∂ tB
∆ tB

∂ u 1
|r⃗ BA (t B ) =
| √ | | | √ |
u ∆ ⃗r BA (t B ) =
∂u
∆ u= ∆u
√ 2 u

∂u ∂u
| |
∂ xB
∆ x B=2 x B ∆ x B
| |
∂ v0 B
∆ v 0 B=2(v 0 B −v0 A )t 2B ∆ v 0 B

8
 v0 B 1 v0B
tB= ∆ t B= ∆ v0B+ 2 ∆ g
g g g

∂u ∂u
| |
∂ v0 A
∆ v 0 A =2(v 0 B −v 0 A ) t 2B ∆ v 0 A
| |
∂ tB
∆ t B=2(v 0 B−v 0 A )2 t B ∆ t B

EJERCICIO 7: Se tiene una placa de forma rectangular, usando un instrumento de

medida con una precisión de 0.1 se midió el largo del rectángulo (25.6 m). Usando
otro instrumento de medida de precisión 0.2 se midió el ancho del rectángulo 112.4
m). Si se desea calcular el perímetro del rectángulo. ¿Cuál sería el error absoluto de
su perímetro?

a=25,6 ± 0,1
b
b=12,4 ± 0,2
a

PERIMETRO
P=2 a+ 2b

P=2 ( 25,6 cm ) +2(24,8 cm)

P=51,2 cm+24,8 cm

P=76 cm.

ERROR ABSOLUTO
Ea P=2( Ea a)+ 2( Ea b)

Ea P=2 ( 0,1 cm ) +2(0,2 cm)

Ea P=0,2cm+0,4 cm

Ea P=0,6 cm.

FINALMENTE
PERIMETRO=76 cm± 0,6 cm

EJERCICIO 8: Considere la siguiente función f ( x )=x 3−x −1.

a) Grafique o tabule la función para hallar sus raíces.

9
 b) Encontrar un intervalo de longitud 1 que aísle la raíz positiva.
c) Sin resolver, calcule cuantas iteraciones se debe realizar para satisfacer una
tolerancia de 0.0015.
d) Haga 3 iteraciones del método de Bisección e indique la solución luego de 4
iteraciones y el máximo error que se espera.
e) Haga 3 iteraciones del Método de falsa posición.
f) Haga 3 iteraciones del Método de punto fijo.
g) Haga 3 iteraciones del Método de Newton-Rapson.

Solución:

a) Tabule : f ( x )=x 3−x−1.

x F(x)
- -25
3
- -7
2
- -1
1
0 -1
1 -1
b) Raíz en (1,2) 2 5
3 23
+ +

-1 1.5 2

b−a 2−1 1
c) MEPn= = n = n ; tol=0.0015
2n 2 2

MEPn ≤tol n=¿ iteraciones

1
≤ 0.0015
2n
666.67 ≤ 2n
2.824 ≤ n log 2
n ≥ 9.381
n=9

Respuesta: Se deben realizar 9 iteraciones para satisfacer una tolerancia de 0.0015

d) Método de bisección:

10
 1° Iteración: (1,2)

c=1.5

f ( c ) =f ( 1.5 )=+ ¿

b−a 2−1
MEP= = =0.5
2 2

+ +

-1 1.5 2

 2° Iteración: (1,1.5)

c=1.25

f ( c ) =f ( 1.25 )=−¿

MEP=0.25

1 1.25 1.5

 3° Iteración: (1.25 ,1.5)

c=1.375

f ( c ) =+¿

MEP=0.125

+ +

1.25 1.375 1.5

 4° Iteración: (1.25 ,1.375)

c=1.3125

f ( c ) =−¿

MEP=0.0625

Solución:c=1.3125 (raíz aproximada)

MEP=0.0625

Raiz ∈ 1.3125 ±0.0625

11
 b−a
MEPn=
2n

2−1
MEP9 = =0.00195
29

e) MÉTODO DE FALSA POSICIÓN:

 1° Iteración: (1,2)

bf ( a ) −af (b)
c=
f ( a ) −f (b)

2 (−1 )−1(5)
c= =1.1667
−1−5

 Error :

MEP=b−a=2−1=1
Error: |f ( c )|=|f ( 1.167 )|=|−0.578|=0.5777
+

1 1.167 2
- -

 2° Iteración: (1.167,2)

2 (−0.578 )−1.167 ( 5 )
c= =1.2533
−0.578−5

 Error:

MEP=2−1.167=0.833

error =|f ( 1.253 )|=|−0.286|=0.2858

1.16 1.253 2
- -
 3° Iteración: (1.253, 2)

2 (−0.286 )−1.253 ( 5 )
c= =1.2937
−0.286−5

MEP=0.7467

Error: |f (1.293 )|=|−0.131|=0.1313

12
 Solución:

La Raíz Aproximada sería: 1.2937

MEP: 0.7467

Error: 0.1313

f) MÉTODO PUNTO FIJO

 Punto fijo: x 0=1.5
Hallamos un G(x) que converja en 1.5:

f ( x )=x 3−x −1 X0 EA

g ( x )=√3 x+1 0 1.5 -

1 1 1.3572 0.1428
g' ( x )= ¿
3
2 1.3309 0.0263
'
g (1.5 )=0.181
3 1.3259 0.005
 0.181<1 (Converge)

 Raiz=1.3259
 Error=0.005

g) MÉTODO DE NEWTON-RAPSON
 1° Iteración:

f ( x )=x 3−x −1 ; x 0=1.5

f ' ( x )=3 x 2−1

f ( x ¿¿ 0)
x 1=x 0− ¿
f '( x¿¿ 0) ¿

0.875
x 1=1.5−
5.75

x 1=1.34783

 2° Iteración:

f (x¿ ¿1)
x 2=x 1− ¿
f '(x ¿¿ 1)¿

0.1007
x 2=1.3478−
4.4499

x 2=1.3252

13
  3° Iteración:

f ( x ¿¿ 2)
x 3=x 2− ¿
f '( x ¿¿ 2)¿

0.00206
x 3=1.3252−
4.26847

x 3=1.3247

EJERCICIO 9: Una determinada sustancia se desintegra según la ecuación A= Pe -0.0248t

donde P es la cantidad inicial en el tiempo t=0 y A la cantidad resultante después de t
años. Si inicialmente se depositan 500 miligramos de dicha sustancia. ¿Cuánto
tiempo habrá de transcurrir para que quede el 1% de ésta? Usar el método que
estime conveniente de tal manera que la solución tenga tres dígitos significativos
hacer su programa en MATLAB de dichos métodos.

Solución:

Cantidad Resultante A = 1%(500) = 5

Cantidad Inicial P = 500

14
Tiempo Transcurrido = t =?

EJERCICIO 10: Probar que la ecuación x3 + cos(x) + 10x = 0 tiene una única solución
real. Aproximar la solución usando el método de Newton-Raphson con 3 iteraciones.

Solución:

15
60

40

20

-20

-40

-60
-3 -2 -1 0 1 2 3
X0 = -0.5

16
EJERCICIO 11: Calcular √24 con 6 decimales exactos, usando el método de falsa
posición.
f ( x )=x 2−24

Tabulando:

x f(x)
3 -15
4 -8 Cambio de signo
5 1
Raíz en (4,5)

 1° Iteración: (4,5)

f ( b )( b−a ) (1)(1)
c=b− =5−
f ( b )−f ( a ) 1+8

c=4.888889

- - +

4 4.888889 5

 2° Iteración: (4.888889;5)

c=4.898876

- - +

4.888889 4.898876 5

 3° Iteración: (4.898876;5)

c=4.898886

- - +

4.898876 4.898886 5

 3° Iteración: (4.898886;5)

f ( b )( b−a ) (1)(5−4.898886)
c=b− =5−
f ( b )−f ( a ) 1+(9.159590∗10−4)

c=4.898979

17
EJERCICIO 12: Aproxime con 10-4 de precisión las raíces de las siguientes ecuaciones
en los intervalos dados, usando el método indicado. haga 3 iteraciones y el programa
respectivo.
a) x 3−2 x 2−5=0 , [ 1,4 ] Bisección
b) x 3−3 x 2−1=0 , [ −4 , 0 ] Falsa posición
π
[ ]
c) x−cosx=0 , 0 , punto fijo
2
π
d) x−0.8−0.2 sinx=0 , 0 , [ ] 2
Newton Rapshon

SOLUCIÓN:

a) f ( x )=x 3−2 x 2−5=0 , [ 1,4 ] Bisección

x F(x)
 1° Iteración: (2, 3) 1 -6
2 -5
2+3 3 4
c= =2.5
2 4 27

3−2
MEP= =0.5
2

f (c) = -1.875
+

2 2.5 3
- -
 2° Iteración: (2.5, 3)

c=2.75

MEP=0.25

f ( c ) =0.671875

+ +

2.5 2.75 3
-

 3° Iteración: (2.5, 2.75)

c=2.625

MEP=0.125

f ( c ) =−0.693359

18
 2.5 0.625 2.75
- -
Sol:

c=2.625

MEP=0.125

b) f ( x )=x 3−3 x 2−1=0 , [−4 , 0 ] Falsa posición

x F(x)
-4 -17
3 -1
-2 3
-1 1
0 -1

 1° Iteración: (-3,-2)

3 (1 )
c=−2− =−2.75
4

MEP=1

f (−3 )=−1

f (−2 )=3

Error=|f (−2.75 )|=0.8906

+ +

-3 -2.75 -2
-

 2° Iteración: (-3, -2.75)

( 0.8906 )( 0.25 )
c=−2.75− ( 1.8906 )
=−2.8678

f (−3 )=−1

f (−2.75 )=0.8906

MEP=0.25

Error=|0.0872|=0.0872

+ +

19
 -3 -2.8678 -2.75
-

 3° Iteración: (-3, -2.8678)

f (−3 )=−1

f (−2.8678 )=0.0872

( 0.0872) ( 0.1322 )
c=−2.8678− ( 1.0872 ) =−2.8784

MEP=0.1322

Error=|0.0075|=0.0075

Sol 1:

c=−2.8784

MEP=0.1322

Error=0.0075

 1° Iteración: (-1,0)

f (−1 ) =1

f ( 0 )=−1

(−1 )( 1 )
c=0− ( −2 )=−0.5

MEP=1

Error=|−0.375|=0.375

-1 -0.5 0
- -

 2° Iteración: (-1, 0.5)

c=−0.6364

MEP=0.5

Error=|−0.0427|=0.0427

20
 -1 -0.6364 0.5
- -

 3° Iteración: (-1, -0.6364)

c=−0.6513

MEP=−0.3636

Error=|−0.0037|=0.0037

Sol 2:

c=−0.6513

MEP=−0.3636

Error=0.0037

c) Met. Punto fijo:

f ( x )=x−cosx=0 [0 ; π/2]

g ( x 1 ) → x=cosx

x f(x) x g(x) E.A E.R

0 -1 0 0 0.8660 - -
π/6 -0.3424 Raíz
1 0.8660 0.6479 0.2181 0.3366
π/4 0.0783 2 0.6479 0.7974 0.1495 0.1875
π/3 0.5472 3 0.7974 0.6986 0.0988 0.1414

Sol: C = 0.6986 E.A = 0.0988 E.R = 0.1414

π
[ ]
d) f ( x )=x−0.8−0.2 sinx=0 , 0 ,
2
Newton Raps h on

X f(x)
0 -0.8
π -0.3764
6
π -0.1560
4
π 0.0740
3

f ( x )=x−0.8−0.2 sinx=0

21
 f ' (x )=1+0.2 cosx

x 1=0

x 2=0− ( −0.8
1.2 )
=0.6667

−9.62 x 10−3
x 3=0.6667− ( 1.1572 )
=0.6750

−2 x 10−5
x 4 =0.6750− ( 1.156 )=0.6750

X=0.6750

EJERCICIO 13: Resolver:

a) Si en la siguiente función: f ( x )=x e x −2 x 3+ 4 x−15=0se despeja el primer x

del primer término. Si se aplica el método del punto fijo, indique si habrá
convergencia o divergencia y de qué tipo.

f ( x )=x e x −2 x 3+ 4 x−15=0
x f (x)
−3 −542
2 x 3−4 x+15 −2 51
G ( x) → x =
ex

|G ' ( x )|<1→ Convergencia

( 6 x 2−4 ) e x −(2 x 3−4 x+15)e x
G ' ( x) =
(e ¿¿ x)2 ¿

G' (−2 )=96 ≠|G ' (x )|< 1

Existe divergencia

G ' ( x)
X 0→ −2
X1 → 51.723 Existe divergencia Osciltoria
X2 → 9.52 x 10 −18

22
 b) Aplica el método de Newton Raphson (4 Iteraciones) e indique si hay o no
convergencia.
f =x e x +e x −6 x2 + 4
'(x )

x 0=−2.5

6.0448
x 1=−2.5− =−2.3202 ERROR = 0.1798
−33.6231

0.4720
x 2=−2.3202− =−2.3036 ERROR = 0.0166
−28.4297

0.0039
x 3=−2.3036− =−2.3035 ERROR = 0.0001
−27.9697

0.0011
x 4 =−2.3035− =−2.3035 ERROR = O
−27.9669

EJERCICIO 14: Si en el ejercicio 12 a, se desea aplicar el punto fijo, encuentre sin

iterar (aplicando el criterio de convergencia) un G(x) que permita hallar una raíz.

f ( x )=x 3−2 x 2−5 ;[1,4 ]→ f ( x )=x (x 2−2 x)−5

1+4
Determino el x o: x o= =2.5
2

Despejando en el 1° Término
5
 G ( x ) : x= 2 → x o=2.5
(x −2 x)

Determino la Función en Matlab

%Metodo del punto fijo con Criterio de Convergencia

% x0=2.5
Clc, clear
G= ‘5/X^2-2*x’
Syms x
dG=diff(G,x)
G=inline(dG)

23
 X0=2.5
if abs(dG(x0)) < 1
dG_X0= abs (dG (x0))
disp (‘Converge’)
else
disp (‘Diverge’)
end

EJERCICIO 15: Con los siguientes valores:

L/r (x) 140 180 220 240

p/a (y) 12800 7500 5000 3800

Donde p/a es la carga en lb/pulg 2 que causa la ruptura de una columna de hierro
dulce con extremos redondeados y L/r es la razón de la longitud de la columna al
mínimo radio de giro de la sección transversal.

a) Encuentre el polinomio de tercer grado que pasa por esos puntos

P3 ( x ) =a1 +a2 x+ a3 x2 +a 4 x3 (Aproximación Polinomial simple), Hallar P(200).

b) Usando polinomios de Lagrange hallar el valor de Y cuando x=200.

Solución:

a) P3 ( x ) =a1 +a2 x+ a3 x2 +a 4 x3

P3 ( 140 )=a 1+140 a 2+19600 a 3+2744000 a 4=12800

P3 ( 180 )=a 1+180 a 2+32400 a 3+5832000 a 4=7500

P3 ( 220 )=a 1+220 a 2+ 48400 a3 +10648000 a4 =5000

P3 ( 240 )=a 1+240 a 2+57600 a 3+ 13824000a 4=3800

Formamos nuestra matriz:

1 140 19600 2744000 12800

|1
1
1
180 32400 5832000
220 48400 5832000
240 57600 13824000
| |
7500
5000
3800

 1°: Multiplicamos fila 1 por -1 y le sumamos la fila 2.

 2°: Multiplicamos fila 1 por -1 y le sumamos la fila 3.
 3°: Multiplicamos fila 1 por -1 y le sumamos la fila 4.

Nos queda la siguiente matriz:

24
 1 140 19600 2744000 12800

|0 40 12800 3088000 −5300

0 80 28800 7904000 −7800
0 100 38000 11080000 −9000
| |
 1°: Multiplicamos fila 2 por -2 y le sumamos la fila 3.
 2°: Multiplicamos fila 3 por -2.5 y le sumamos la fila 4.

Nos queda la siguiente matriz:

1 140 19600 2744000 12800

|0
0
0
40 12800 3088000 −5300
0
0
3200 1728000 2800
6000 3360000 4250
| |
Y por último multiplicamos la fila 3 por -1.875 y le sumamos la fila 4.
1 140 19600 2744000 12800

| 0
0
0
40 12800 3088000 −5300
0
0
3200 1728000 2800
0 120000 −1000
| |
120000 a4 =−1000 → a4 =−0.00833

3200 a3 +1728000 (−0.00833 )=2800→ a3 =5.3732

40 a2+ 12800 ( 5.3732 ) +3088000 (−0.00833 )=−5300
a 2=−1208.848
a 1+140 (−1208.848 ) +19600 ( 5.3732 )+ 2744000 (−0.00833 )=12800

a 1=99581.52

Ahora hallamos P (200):

P3 ( x ) =a1 +a2 x+ a3 x2 +a 4 x3

Reemplazando, tenemos:

P3 ( x ) =99581.52−1208.848 x+5.3732 x 2−0.00833 x 3

P3 ( 200 )=99581.52−1208.848(200)+5.3732(200)2−0.00833(200)3

P3 ( 200 )=6099.92=6100
25
 b) P3 ( x ) =L1 y 1 + L2 y 2 + L3 y3 + L 4 y 4
Hallar y para x = 200

( x−180)(x−220)(x−240)

−1
L1 = → L1=
(140−180)(140−220)(140−240) 20

( x−140)(x−220)(x−240) 1
 L2 =
(180−140)(180−220)(180−240)
→ L2=
2

( x−140)(x −180)(x −240) 3

 L3 =
(220−140)(220−180)(220−240)
→ L3=
4

(x−140)( x−180)( x−220) −1

 L4 =
(240−140)(240−180)(240−220)
→ L4 =
5

−1 1 3 1
 P3 ( 200 )= ( 12800 ) + ( 7500 ) + ( 5000 )− ( 3800 )
20 2 4 5

 P3 ( 200 )=6100

EJERCICIO 16: En una reacción química, la concentración del producto CB cambia con
el tiempo como se indica en la tabla de abajo. Calcule la concentración CB cuando
t=0.82
Ca 0.00 0.30 0.55 0.80 1.10 1.15
t 0.00 0.10 0.40 0.60 0.80 1.00
Trabajar con un polinomio de segundo grado. Usar polinomios de Lagrange.

SOLUCIÓN:

P2 ( x ) =L1 y 1 + L2 y 2 + L3 y 3

x−x 2 x −0.8 0.82−0.8

∗x−x 3 ∗x−1 ∗0.82−1
x 1−x 2 0.6−0.8 −0.2
L1 = → = =−0.045
x 1−x 3 0.6−1 −0.4

x−x 1 x −0.6 0.82−0.6

∗x−x 3 ∗x−1 ∗0.82−1
x 2−x 1 0.8−0.6 0.2
L2 = → = =0.99
x 2−x 3 0.8−1 −0.2

x−x 1 x−0.6 0.82−0.6

∗x−x 2 ∗x −0.8 ∗0.82−0.8
x 3−x 1 1−0.6 0.4
L3 = → = =¿
x 3 −x2 1−0.8 0.2

26
P2 ( x ) =−0.045(0.80)+0.99(1.10)+0.055(1.15)

P2 ( x ) =−0.036+1.089+0.06325

P2 ( x ) =1.11625

EJERCICIO 17: De acuerdo con la taba siguiente se muestran los valores f(x,y):

Y X
0 40 80 100
4 10.881 10.376 10.063 9.531
12 12.160 11.013 10.786 9.663
20 13.997 12.801 11.570 10.451
28 15.846 13.652 12.418 11.301

a) Calcule f(40,15) (primer grado) aproximacion simple

Y X=40
x1 12 11.013 y 1 P1 ( x )=a1 +a2 x
x2 20 12.801 y 2

P12 ( 12 )=a1 +a2 ( 12 )=11.013

P20 ( 20 )=a1 + a2 ( 20 ) =12.801

[ 112012⋮⋮ 11.013
12.801 ] F (1 ) ×−1+ F (2) → [ 112 ⋮11.013 ]
0 8 ⋮ 1.788

 8 a 2=1.788 → a2=0.2235

 a 1+12 ( a 2) =11.013
a 1+12 ( 0.2235 )=11.013
a 1=8.331

ENTONCES:
P1 ( x )=8.331+ 0.2235 x
P15 ( 15 )=8.331+0.2235 ( 15 )

RESPUESTA

P15 ( 15 )=11.6835

b) Calcule f (50,28) (tercer grado) polinomios lagrange

27
 X =50
X Y=28

x1 0 15.846 y1 P3 ( x ) =L1 y 1 + L2 y 2 + L3 y3 + L 4 y 4
x 2 40 13.652 y2
x 3 80 12.418 y3
x 4 100 11.301 y4

x−x 2 x−x 3 x −x 4
L1 = × ×
x1 −x2 x 1−x 3 x 1−x 4

i=1
j= 1,2,3,4

50−40 50−80 50−100 −3

L1 = × × =
0−40 0−80 0−100 64

x−x 1 x−x 3 x−x 4

L2 = × ×
x 2−x1 x 2−x 3 x 2−x 4

i=2
j= 1,2,3,4
50−0 50−80 50−100 25
L2 = × × =
40−0 40−80 40−100 32
x−x 1 x−x 2 x−x 4
L3 = × ×
x 3−x 1 x3 −x2 x 3−x 4
i=3
j= 1,2,3,4
50−0 50−40 50−100 25
L3 = × × =
80−0 80−40 80−100 64

x−x 1 x−x 2 x−x3

L4 = × ×
x 4 −x 1 x 4−x 2 x 4−x 3
i=4
j= 1,2,3,4
50−0 50−40 50−80 −1
L4 = × × =
100−0 100−40 100−80 8

RESPUESTA:
−3 25 25 1
P3 ( 50 )= ( 15.846 ) + ( 13.652 ) + (12.418 )− ( 11.301 ) =13.361
64 32 64 8

c) Calcule f (90,25) (Segundo grado) polinomios de lagrange

28
 Y X
40 80 90 100
12 11.013 10.786 9.663
20 12.801 11.570 10.451
25 A b d c
28 13.652 12.418 11.301

a.
X = 25
X Y

x 1 12 11.013 y 1 P2 ( x ) =L1 y 1 + L2 y 2 + L3 y 3
x 2 20 12.801 y 2
x 3 28 13.652 y 3

x−x 2 x−x 3
L1 = ×
x1 −x2 x 1−x 3
i=1
j= 1,2,3
25−20 25−28 −15
L1 = × =
12−20 12−28 128

x−x 1 x−x 3
L2 = ×
x 2−x1 x 2−x 3
i=2
j= 1,2,3
25−12 25−28 39
L2 = × =
20−12 20−28 64

x−x 1 x−x 2
L3 = ×
x 3−x 1 x3 −x2
i=3
j= 1,2,3
25−12 25−20 65
L3 = × =
28−12 28−20 128

RESPUESTA:
−15 39 65
P2 ( 25 )= ( 11.013 )+ ( 12.801 ) + ( 13.652 )=13.443
128 64 128

29
 b.
X Y

x 1 12 10.786 y 1 P2 ( x ) =L1 y 1 + L2 y 2 + L3 y 3
x 2 20 11.570 y 2
x 3 28 12.418 y 3

x−x 2 x−x 3
L1 = ×
x1 −x2 x 1−x 3
i=1
j= 1,2,3
25−20 25−28 −15
L1 = × =
12−20 12−28 128

x−x 1 x−x 3
L2 = ×
x 2−x1 x 2−x 3
i=2
j= 1,2,3
25−12 25−28 39
L2 = × =
20−12 20−28 64

x−x 1 x−x 2
L3 = ×
x 3−x 1 x3 −x2
i=3
j= 1,2,3
25−12 25−20 65
L3 = × =
28−12 28−20 128

RESPUESTA:
−15 39 65
P2 ( 25 )= ( 10.786 ) + ( 11.570 ) + ( 12.418 )=12.0925
128 64 128

c.

X Y

x 1 12 9.663 y 1 P2 ( x ) =L1 y 1 + L2 y 2 + L3 y 3
x 2 20 10.451 y 2
x 3 28 11.301 y 3

30
 x−x 2 x−x 3
L1 = ×
x1 −x2 x 1−x 3
i=1
j= 1,2,3
25−20 25−28 −15
L1 = × =
12−20 12−28 128

x−x 1 x−x 3
L2 = ×
x 2−x1 x 2−x 3
i=2
j= 1,2,3
25−12 25−28 39
L2 = × =
20−12 20−28 64

x−x 1 x−x 2
L3 = ×
x 3−x 1 x3 −x2
i=3
j= 1,2,3
25−12 25−20 65
L3 = × =
28−12 28−20 128

RESPUESTA:
−15 39 65
P2 ( 25 )= ( 9.663 )+ ( 10.451 )+ (11.301 )=10.975
128 64 128

d.
X = 90

X Y

x 1 40 13.442 y 1 P2 ( x ) =L1 y 1 + L2 y 2 + L3 y 3
x 2 80 12.0925 y 2
x 3 100 10.975 y 3

x−x 2 x−x 3
L1 = ×
x1 −x2 x 1−x 3
i=1
j= 1,2,3
90−80 90−100 −1
L1 = × =
40−80 40−100 24

31
 x−x 1 x−x 3
L2 = ×
x 2−x1 x 2−x 3
i=2
j= 1,2,3
90−40 90−100 5
L2 = × =
80−40 80−100 8

x−x 1 x−x 2
L3 = ×
x 3−x 1 x3 −x2
i=3
j= 1,2,3
90−40 90−80 5
L3 = × =
100−40 100−80 12

RESPUESTA FINAL
−1 5 5
P2 ( 25 )= ( 13.443 ) + ( 12.0925 ) + ( 10.975 )=11.571
24 8 12

d) Calcule el valor de (y) para un x=60 cuyo f(x,y)=11.29(primer grado) El metodo

que prefiera

Y X
40 60 80 100
12 11.013 a 10.786 9.663
c 11.29
20 12.801 b 11.570 10.451
28 13.652 12.418 11.301

a.
X=60

X Y

x 1 40 11.013 y 1 P1 ( x )=L1 y 1+ L2 y 2
x 2 80 10.786 y 2

x−x 2
L1 =
x1 −x2

32
i=1
j= 1,2
60−80 1
L1 = =
40−80 2

x−x 1
L2 =
x 2−x1
i=2
j= 1,2
60−40 1
L2 = =
80−40 2
RESPUESTA:
1 1
P1 (25 )= ( 11.013 ) + ( 10.786 ) =10.8995
2 2

b.
X=60
X Y

x 1 40 12.801 y 1 P1 ( x )=L1 y 1+ L2 y 2
x 2 80 11.570 y 2

x−x 2
L1 =
x1 −x2
i=1
j= 1,2
60−80 1
L1 = =
40−80 2

x−x 1
L2 =
x 2−x1
i=2
j= 1,2
60−40 1
L2 = =
80−40 2
Rpta:
1 1
P1 (25 )= ( 12.801 )+ (11.570 )=12.1885
2 2

c.
X=11.29
33
 X Y

x 110.8995 12 y 1 P1 ( x )=L1 y 1+ L2 y 2
x 212.1885 20 y 2

x−x 2
L1 =
x1 −x2
i=1
j= 1,2
11.29−12.1885
L1 = =0.697
10.8995−12.1885

x−x 1
L2 =
x 2−x1
i=2
j= 1,2
11.29−10.8995
L2 = =0.303
12.1885−10.8995

RESPUESTA FINAL
P1 (11.29 )=0.697 ( 12 ) +0.303 ( 20 ) =14.424

EJERCICIO 18: Obtenga la segunda derivada evaluada en x=37 para la funcion que se
da enseguida

X 1 1.8 3 4.2 5 6.5

F(x) 3 4.34536 6.57735 8.88725 10.44721 13.39223
Trabajar con un polinomio de segundo grado. Usar aproximación simple
Resolución P2 ( x ) =a1 +a2 x+ a3 x 2
X Y

x1 1.8 4.34536 y 1 x2 3 6.57735 y 2

x3 4.2 8.88725 y 3
P2 (1.8 )=a1 + a2 ( 1.8 ) + a3 ( 1.8 )2 =4.34536
P2 (3 )=a 1+ a2 ( 3 )+ a3 ( 3 )2=6.57735
P2 ( 4.2 )=a1+ a2 ( 4.2 ) +a3 ( 4.2 )2=8.88725

34
 11.8 3.24 ⋮ 4.34536
[ 13 9 ⋮ 6.57735
1 4.217.64 ⋮ 8.88725 ]
F ( 1 ) ×−1+ F ( 2 )=¿

11.8 3.24 ⋮ 4.34536

[ ]
0 1.25.76 ⋮ 2.23199 F ( 1 ) ×−1+ F ( 3 )=¿
1 4.217.64 ⋮ 8.88725

11.8 3.24 ⋮ 4.34536

[ ]
0 1.2 5.76 ⋮2.23199 F ( 2 ) ×−2+ F ( 3 )=¿
0 2.4 14.4 ⋮ 4.54189

11.8 3.24 ⋮ 4.34536

[0 1.25.76 ⋮ 2.23199
0 0 2.88 ⋮0.07791 ]
a 3 ( 2.88 ) =0.07791

a 3=0.02705

a 2 ( 1.2 ) + ( 0.02705 ) ( 5.76 )=2.23199

a 2=1.73015

a 1+ ( 1.73015 ) ( 1.8 ) + ( 0.02705 ) ( 3.24 ) =4.34536

a 1=1.14345

Ecuacion:
P2 ( x ) =1.14345+ 1.73015 x +0.02705 x 2

Primera derivada:

P' 2 ( x )=1.73015+0.0541 x
Segunda derivada:
P' '2 ( x )=0.0541

Para X = 37
P' '2 ( 37 )=0.0541

35