UNIVERSIDAD NACIONAL AUTONOMA DE HONDURAS

Facultad de Ciencias - Escuela de Matemática

Guía de Ejercicios #2 de Análisis Numérico

Problema 1 Construya los polinomios interpolantes de Lagrange para las siguientes funciones y
obtenga una cota del error absoluto en el intervalo [ x0 , xn ].

a) f ( x ) = e2x cos 3x, x0 = 0, x1 = 0.3, x2 = 0.6, n = 2

b) f ( x ) = sen(ln x ), x0 = 2.0, x1 = 2.4, x2 = 2.6, n = 2

c) f ( x ) = ln x, x0 = 1, x1 = 1.1, x2 = 1.3, x3 = 1.4, n = 3

Problema 2 Sea P3 ( x ) el polinomio interpolante para los datos (0, 0), (0.5, y), (1, 3), (2, 2). Encuen-
tre y si el coeficiente de x3 en P3 ( x ) es 5.46

Problema 3 Sea P3 ( x ) el polinomio interpolante de los datos

x 8 8.5 9 9.5
g( x ) 16.63553 18.19056 19.77502 α

Determine el valor de α si el coeficiente de x en P3 ( x ) es 1.66546.

Problema 4 Considere la función S definida como:

 1 + x − x2

0≤x<1
S= −4 + 16x − 15x2 + 4x3 1≤x<2
52 − 68x + 27x2 − 3x3 2≤x<3


¿S es un spline cúbico con nodos x0 = 0, x1 = 1, x2 = 2 y x3 = 3?

Problema 5 Determine los trazadores cúbicos sujetos en los nodos dados y aproxime f ( x ) y f 0 ( x )
a las siguientes funciones para el valor dado de x. Además determine el error real.

a) f ( x ) = sen(e x − 2); x0 = 0.8, x1 = 1.0; x = 0.9

b) f ( x ) = x3 + 4.001x2 + 4.002x + 1.101; x0 = −0.5, x1 = −0.25, x2 = 0; x = −1/3

Primer período - 2020 Lic. Angel Rivera

Problema 6 Un trazador cúbico natural S de la función f está definido por

 S0 ( x ) = 1 + 2x − x3

0≤x<1
S( x ) =
S1 ( x ) = 2 + b ( x − 1 ) + c ( x − 1 ) 2 + d ( x − 1 ) 3 1≤x≤2


Obtenga b, c y d.

Problema 7 Obtenga los polinomios de mínimos cuadrados de primero, segundo y tercer grado
para los datos de la tabla anexa. En cada caso calcule el error real E.

xi 1.0 1.1 1.3 1.5 1.9 2.1

yi 1.84 1.96 2.21 2.45 2.94 3.18

Problema 8 Sea f ( x ) = tan( x ) y usando h = 0.1, 0.0001

f ( x +h)−2 f ( x )+ f ( x −h)
a) Calcule aproximaciones a f 00 (0.8) usando f 00 ( x ) ≈ h2
, use 6 cifras decima-
les significativas o más.
2 f ( x )−5 f ( x +h)+4 f ( x +2h)− f ( x +3h)
b) Calcule aproximaciones a f 00 (0.8) usando f 00 ( x ) ≈ h2
, use 6 ci-
fras decimales significativas o más.
2 f ( x )−5 f ( x −h)+4 f ( x −2h)− f ( x −3h)
c) Calcule aproximaciones a f 00 (0.8) usando f 00 ( x ) ≈ h2
, use 6 ci-
fras decimales significativas o más.

d) Determine el error real para los resultados obtenidos en a), b) y c).

Problema 9 Deduzca las fórmulas presentadas en el ejercicio anterior.

Problema 10 Adapte la fórmula

f ( x + h) − f ( x − h)
f 0 (x) ≈
2h
Para calcular las derivadas parciales f x ( x, y) y f y ( x, y). (Sugerencia: recuerde la definición de una
derivada parcial).

Problema 11 Usando la fórmula de Taylor para f ( x + h), f ( x − h), f ( x ), f ( x + 2h), f ( x − 2h), de-
duzca la fórmula de diferencias centradas:
f ( x + 2h) − 4 f ( x + h) + 6 f ( x ) − 4 f ( x − h) + f ( x − 2h)
f (4) ( x ) ≈
h4

Primer período - 2020 Lic. Angel Rivera

Problema 12 Aproxime las siguientes integrales aplicando la regla del trapecio, simpson y punto
medio.
Z 1
a) x4 dx
0
Z π/4
b) x sen xdx
0

Obtenga una cota del error para cada fórmula de cuadratura y compárela con el error real.

Problema 13 Aproxime las siguientes integrales usando la cuadratura gaussiana con n = 2 y

obtenga el error real.
Z 1.5
a) x2 ln xdx
1
Z 1
b) x4 dx
0
Z π/4
c) x sen xdx
0

Teóricos y aplicaciones

x−y
Problema 14 Sea f ( x, y) = . Usando h = 0.1 y h = 0.01 calcule aproximaciones a f x (2, 3) y
x+y
f y (2, 3) usando las fórmulas (problema 10).

Problema 15 Se definen las medias artiméticas x̄ y ȳ para los puntos {( xk , yk )}kN=1 como

1 N 1 N
N k∑ N k∑
x̄ = xk ȳ = yk
=1 =1

Pruebe que el punto ( x̄, ȳ) está en la recta de regresión de minimos cuadrados determinada por
el conjunto de puntos dado.

Problema 16 Obtenga el grado de precisión de la fórmula de cuadratura

Z 1 √ ! √ !
− 3 3
f ( x )dx = f +f
−1 3 3

Primer período - 2020 Lic. Angel Rivera

Problema 17 Determine las constantes a, b, c y d que producirán una fórmula de cuadratura
Z 1
f ( x )dx = a f (−1) + b f (1) + c f 0 (−1) + d f 0 (1)
−1

cuyo grado de precisión es 3.

Problema 18 Al medir la velocidad (con un tubo de Pitot) en una tubería circular de diámetro
interior de 20 cm, se encontró la siguiente información:

v(cm/s) 600 550 450 312 240

r (cm) 1 3 5 7 8

donde r es la distancia en cm medida a partir del centro del tubo.

a) Obtenga la curva v = f (r ) que mejor se ajusta a estos datos experimentales (se sabe que el
perfil de velocidades en una tubería generalmente es de tipo parabólico).

b) Calcule la velocidad en el punto r = 4 cm.

c) Contruya la aproximación por mínimos cuadrados de la forma be ax (linealice) y calcule el

error.

Problema 19 Un investigador obtiene de un experimento los siguientes datos que se muestran a

continuación:
x 1.0 2.0 3.0 4.0 5.0
y 2.0 5.0 10.0 17.0 26.0

D
Se sabe que dichos datos pueden modelarse por la curva de ajuste y = , donde D y C
x+C
son parámetros.Determine la linealización de la curva de ajuste y la recta de regresión lineal de
mínimos cuadrados para estimar D y C.

Problema 20 Un investigador obtiene de un experimento los siguientes datos que se muestran a

continuación:
x 1.2 1.4 1.6 1.8 2.0 2.2
y 0.5073 0.5616 0.6098 0.6452 0.6824 0.7110

x
Se sabe que dichos datos pueden modelarse por la curva de ajuste y = , donde c
ln(c/d + de x )
y d son parámetros. Determine la linealización de la curva de ajuste y la recta de regresión de
mínimos cuadrados para estimar c y d.

Primer período - 2020 Lic. Angel Rivera

Problema 21 La longitud de arco de una curva continua y derivable y = f ( x ) se define por
b
Z q
L= 1 + [ f 0 ( x )]2 dx
a

Sea f ( x ) = x1/3 + x2/3 aproxime la longitud de curva de f ( x ) de 0 ≤ x ≤ 2 mediante cuadratura

gaussiana con dos puntos y utilice para aproximar la derivada la fórmula de 5 puntos,

1
f 0 (x) = [ f ( x − 2h) − 8 f ( x − h) + 8 f ( x + h) − f ( x + 2h)]
2h
con tamaño de paso h = 0.10

Problema 22 Mediante los interpolantes de trazadores cúbicos se va a aproximar la parte superior

de una camioneta. Se escalo la imagen de la camioneta y se obtuvo la información recolectada que
se muestra en la figura.

Con la información proporcionada, construya

a) el trazador cúbico sujeto del punto A al punto B.

b) el trazador cúbico sujeto del punto B al punto C.

c) el trazador cúbico sujeto del punto A al punto C (use la versión matricial para determinar
el trazador).

Primer período - 2020 Lic. Angel Rivera

Problema 23 Utilice la transformación t = x −1 y use la cuadratura gaussiana para aproximar las
siguientes integrales con n = 3 y n = 4.
∞
Z
1
a) dx
1
x2 + 9
∞
Z !
1
b) x −4 sen dx.
1
x

Problema 24 Use cuadratura gaussiana de dos puntos para aproximar la siguiente integral doble
Z 2 Z 2x
( x2 + y3 )dydx
1 x

Problema 25 Las áreas de sección transversal (A) son requeridas para diferentes tareas en la inge-
niería de abastecimiento de aguas; entre otras, pronósticos de inundación y diseño de reservorios.
Un ejemplo de una corriente común con su sección transversal se muestra en la figura. Los puntos
representan ubicaciones donde se ancló un bote y tomó lecturas a diferentes profundidades.

Estime el área de la sección transversal a partir de estos datos utlizando:

a) Regla compuesta del trapecio con h = 2 m y h = 4 m.

R∞
Problema 26 La integral impropia 0 f ( x )dx no puede convertirse en una integral con límites
finitos por medio de la sustitución t = x −1 porque el límite en cero se vuelve infinito. El problema
R∞ R1 R∞
se resuelve escribiendo primero 0 f ( x )dx = 0 f ( x )dx + 1 f ( x )dx. Aplique este método para
aproximar las siguientes integrales usando la cuadratura gaussiana con n = 3:
Z ∞
1
a) dx
0 1 + x4

Primer período - 2020 Lic. Angel Rivera

 Z ∞
1
b) dx
0 (1 + x 2 )3
∞
Z
2
c) e− x dx
−∞

Problema 27 La velocidad hacia arriba de un cohete se puede calcular con la siguiente fórmula:
!
m0
v = u ln − gt
m0 − qt

donde v = velocidad hacia arriba, u = velocidad a la cual se expulsa el combustible relativo al

cohete, m0 = masa inicial del cohete en el tiempo t = 0, q = razón de consumo de combustible
y g = aceleración hacia abajo debido a la gravedad (se supone constante e igual a 9.8 m/s2 ). Si
u = 2000 m/s, m0 = 150000 kg y q = 2600 kg/s. Use la cuadratura gaussiana con 3 puntos para
determinar que tan alto volará el cohete en 30 segundos.

Problema 28 La interpolación iterada inversa consiste en aproximar la raíz utilizando la inter-

polación de Lagrange para x = f −1 (y), es decir, suponga que f ∈ C1 [ a, b], f 0 ( x ) 6= 0 en [ a, b] y que
f tiene un cero p en [ a, b]. Sean x0 , · · · , xn , n + 1 números distintos en [ a, b] con f ( xk ) = yk para
cada k = 0, 1, · · · , n. Si quiere aproximar p, construya el polinomio interpolante de grado n en los
nodos y0 , y1 , · · · , yn para f −1 . Puesto que yk = f ( xk ) y 0 = f ( p), se deduce que f −1 (yk ) = xk y
p = f −1 (0 ).
Use la interpolación iterada inversa para obtener una aproximación a la solución de x − e− x = 0,
por medio de los datos

x 0.3 0.4 0.5

y = x − e− x -0.44 -0.27 -0.11

Bibliografía:

Burden, R. Faires, D. (2011) Análisis Numérico (9na. Edición), Canadá: Cengage Learning
Editores.

Mathews, J.H., & Fink, K.D. (2000). Métodos numéricos con Matlab. Madrid: Prentice Hall
Madrid.

Steven, C., & Raymond, C. (2006). Métodos Numéricos para ingenieros. Inc. México. McGraw-
Hill.

Primer período - 2020 Lic. Angel Rivera