Matemática Aplicada A La Estadística - Farmacia U Sevilla PDF
Matemática Aplicada A La Estadística - Farmacia U Sevilla PDF
Matemática Aplicada A La Estadística - Farmacia U Sevilla PDF
2. Interpolación polinómica 37
2.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
2.2. Análisis previo (existencia y unicidad) . . . . . . . . . . . . . . . . . . . . . . . 37
2.3. Método de las Diferencias Divididas . . . . . . . . . . . . . . . . . . . . . . . . 38
3
4 ÍNDICE GENERAL
4. Programación lineal 55
4.1. Introducción / motivación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
4.2. Aplicaciones básicas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
4.3. Resolución por el método geométrico . . . . . . . . . . . . . . . . . . . . . . . . 57
4.3.1. Análisis previo: Regiones factibles . . . . . . . . . . . . . . . . . . . . . 57
4.3.2. Método geométrico en regiones factibles acotadas . . . . . . . . . . . . . 59
4.3.3. Sobre la existencia de solución y regiones factibles no acotadas . . . . . 59
Tema 1
Reales R = “clausura/completado” de Q
(N ⊂ Z ⊂ Q ⊂ R) recta real; axioma del supremo = Teor. de Cantor = Teor. de Bolzano
imposibilidad de resolver ecuaciones como x2 + 1 = 0.
Complejos C = {a + bi : a, b ∈ R, i2 = −1}
(N ⊂ Z ⊂ Q ⊂ R ⊂ C) Toda ecuación algebraica tiene solución en C
5
6 TEMA 1. FUNCIONES, LÍMITES, CONTINUIDAD Y DERIVABILIDAD
Propiedades:
Reflexiva: ∀a ∈ Q : a ≤ a.
Antisimétrica: a ≤ b, b ≤ a ⇒ a = b.
Transitiva: a ≤ b, b ≤ c ⇒ a ≤ c.
Además, dicha relación es total: ∀a, b ∈ Q, a ≤ b ó b ≤ a.
La relación de orden es compatible con la suma y el producto:
Para “completar” la recta hasta llenarla debemos ampliar hasta los números reales. [más
adelante veremos axiomas sobre su construcción: del supremo, y el Teorema de Bolzano.]
R es un cuerpo ordenado
1.1. NÚMEROS REALES Y NOCIONES BÁSICAS SOBRE FUNCIONES 7
En todo R se pueden definir operaciones suma (+) y producto (·) que extienden las
definidas antes, y una relación de orden (≤) que también es total.
Dados a, b ∈ R, diremos que
a ≥ b si b ≤ a.
a < b si a ≤ b y a 6= b.
a > b si b < a.
Propiedades de R:
∀x ∈ R, x > 0, ∃n ∈ N : x < n.
Entre dos números reales distintos siempre existe un racional, y en consecuencia, infinitos
racionales.
8 TEMA 1. FUNCIONES, LÍMITES, CONTINUIDAD Y DERIVABILIDAD
Subconjuntos de R. Intervalos
Definición 1.4. Se llama dominio de una función f (lo denotaremos por Dom f ) al conjunto
de valores para los que está bien definida f (x) :
f : Domf ⊂ R → R : x 7→ f (x).
Habrá ocasiones en que nos den una fórmula y debamos hallar el dominio en que es válida
aplicarla:
√
f (x) = x Domf = R+ .
1
f (x) = Domf = R\{0}.
x
Imagen o rango o recorrido de f ≡ valores que toma f :
y
6 y = x2
¡µ y=x
¡
¡
¡
¡ √
¡ y= x
¡
¡
¡
¡
¡ -
x
Simetrı́a respecto la bisectriz
-
• •
x
A • B
¡ b b
¢
V V = − 2a , p(− 2a )
1.2. FUNCIONES ELEMENTALES 11
√
−b ± b2 − 4ac
A, B =
2a
Para grados superior a 2 debemos esperar a tener más herramientas para conocer su rep-
resentación gráfica.
Para el estudio del signo de un polinomio, resultará muy útil conocer procedimientos de
factorización. Recordamos brevemente con ejemplos cómo se aplica la Regla de Ruffini para
dividir (y eventualmente factorizar) un polinomio de grado cualquiera entre otro de grado
uno.
Se toman los coeficientes del polinomio del numerador y el elemento opuesto (cambio de
signo) del coeficiente de grado cero del denominador y opera como en el ejemplo (se baja el
primer coeficiente, multiplica, se sube el resultado y se suma):
x3 − 3x2 − 7x − 8
Ejemplo 1.8. Queremos calcular el cociente
x−5
1 -3 -7 8
5 5 10 15
1 2 3 7
Observa los números obtenidos a través de la operación anterior. Salvo el último (recuadra-
do), que es el resto, los demás, denotan los coeficientes de un polinomio un grado menor, o
sea, 2:
Esto nos dice que
x3 − 3x2 − 7x − 8 = (x − 5)(x2 + 2x + 3) + 7
El caso interesante se produce cuando el resto es cero (en lugar de 7), eso dice que cierto
número (en este caso habrı́a sido 5) es un cero o raı́z del polinomio y, por tanto, que éste
puede factorizarse como (x − 5) por otro polinomio de un grado menos.
Ejemplo 1.9. Sea q(x) = x4 + x3 − 7x2 + 4. Puedes comprobar que q(2) = 0, eso indica
que q(x) = (x − 2)r(x) con r(x) otro polinomio de grado 3. Para hallarlo desarrollamos por
Ruffini (ojo, hay que poner los coeficientes de todos los monomios, incluidos 0 por aquéllos
que no aparecen):
1 1 -7 0 4
2 2 6 -2 -4
1 3 -1 -2 0
Comprueba, desarrollando el producto, que se tiene la siguiente igualdad:
Funciones racionales
p(x)
f (x) = con p y q funciones polinómicas.
q(x)
Domf = {x ∈ R : q(x) 6= 0} = R\{x ∈ R : q(x) = 0}.
Función exponencial
f (x) = ax , a ∈ R, a > 0, x ∈ R. Domf = R.
Propiedades:
ar+s = ar as , (ar )s = ar·s , (a · b)r = ar · br
Si x, y ∈ R, x < y :
ax < ay cuando a > 1 (creciente).
ax > ay cuando a < 1 (decreciente).
2x
50
45
40
35
30
25
20
15
10
−6 −4 −2 0 2 4 6
x
Propiedades:
x
loga 1 = 0, loga a = 1, loga (x · y) = loga x + loga y, loga (xm ) = m · loga x, loga y = loga x −
loga y.
1.2. FUNCIONES ELEMENTALES 13
log(x)
2
1.5
0.5
−0.5
−1
−1.5
−2
0 1 2 3 4 5 6
x
Funciones trigonométricas
sin(x)
0.5
−0.5
−1
−6 −4 −2 0 2 4 6
x
cos(x)
0.5
−0.5
−1
−6 −4 −2 0 2 4 6
x
Función tangente:
sen x
f (x) = tg x = .
π
cos x
Domf = R\{ 2 + kπ : k ∈ Z}, Imf = R, π−periódica.
tan(x)
−2
−4
−6
−6 −4 −2 0 2 4 6
x
1 1 1
cosec x = , sec x = , cotg x = .
sen x cos x tg x
1.3. LÍMITES DE FUNCIONES 15
Ejemplo 1.12. f : R\{0} → R : x 7→ xsen x. Se tiene que lı́m f (x) = 0 (basta tomar δ = ε).
x→0
Comprobar
½ que ocurre lo mismo con la función
x si x ∈ Q,
g(x) =
0 si x 6∈ Q.
1
Ejemplo 1.13. lı́m = +∞.
x→0 |x|
1 1
Ejemplo 1.14. lı́m = lı́m = 0.
x→+∞ x x→−∞ x
16 TEMA 1. FUNCIONES, LÍMITES, CONTINUIDAD Y DERIVABILIDAD
A veces no podemos asegurar la existencia de lı́mite en torno a un punto, pero sı́ estudiar
los lı́mites laterales (por la derecha e izquierda respectivamente):
Ejemplo 1.15. La función parte entera E : R → R verifica lı́m E(x) = n pero lı́m E(x) =
x→n+ x→n−
n − 1.
Teorema 1.16. Dada f : S ⊂ R → R se cumple que lı́m f (x) existe y vale l ∈ R si y solo
x→x0
si existen los lı́mites laterales lı́m f (x) y lı́m f (x) y tienen el mismo valor.
x→x+
0 x→x−
0
El resultado, leı́do de forma negativa, dice que si los lı́mites laterales no coinciden, no
existe lı́mite de la función en ese punto.
Teorema 1.19. Sean f, g : S ⊂ R → R tales que existen l, l0 ∈ R, y lı́m f (x) = l, lı́m g(x) =
x→x0 x→x0
l0 . Entonces existen los siguientes lı́mites:
lı́m [f (x) ± g(x)] = l ± l0 ,
x→x0
lı́m f (x) · g(x) = l · l0 ,
x→x0
f (x) l
Si l0 6= 0 lı́m = 0.
x→x0 g(x) l
Observación 1.20 (Indeterminaciones). Cuando los valores l y l0 no son reales, sino +∞
ó −∞, nos encontraremos con problemas. Se trata de expresiones que no tienen a priori
un valor concreto (por ello las llamaremos indeterminaciones), sino que hay que hallarlo
explı́citamente en cada caso. Expresiones indeterminadas que pueden aparecer son:
∞ 0
1∞ , 0 · ∞, ∞ − ∞, , , ∞0 , 00 .
∞ 0
Algunas que podremos resolver en este mismo tema usando las propiedades precedentes son
∞−∞ (por ejemplo sacando factor común o multiplicando por el conjugado) o ∞ ∞ (será simple
si se trata de funciones racionales). Un dominio total de las indeterminaciones requerirá her-
ramientas del próximo tema de derivación, como por ejemplo la Regla de L’Hôpital.
x0 ∈6 Domf.
6 ∃ lı́m f (x).
x→x0
∃ lı́m f (x) 6= f (x0 ).
x→x0
c) Discontinuidad de salto infinito: cuando alguno de los lı́mites laterales (o ambos) valen
±∞.
1 1
Ejemplos: f (x) = , g(x) = .
x |x|
f (S) = {f (x) : x ∈ S}
Teorema 1.28.
• (continuidad implica acotación local): f : (a, b) → R, f continua en x0 ⇒ ∃δ > 0 : f
acotada en (x0 − δ, x0 + δ).
• Si a ∈Domf, y f es continua en a por la derecha, entonces ∃δ > 0 : f acotada en
[a, a + δ). (Resultado análogo para el extremo b).
• (continuidad en intervalo cerrado y acotado sı́ implica acotación) f continua en [a, b] ⇒
f acotada en [a, b].
Máximos y mı́nimos
Definición 1.29. Sea f : S ⊂ R → R. Se dice que f tiene un máximo absoluto en c ∈ S
(resp. mı́nimo) si f (x) ≤ f (c) para todo x ∈ S (resp. f (x) ≥ f (c)).
Se dice que tiene un máximo (resp. mı́nimo) relativo o local en c ∈ S cuando existe δ > 0
tal que f (x) ≤ f (c) (resp. f (x) ≥ f (c)) para todo x ∈ (c − δ, c + δ).
Teorema 1.30.
•[Weierstrass] Sea f : [a, b] → R continua. Entonces f alcanza máximo y mı́nimo abso-
lutos. (Nota: el resultado es falso si el intervalo de partida no es cerrado o no es acotado,
incluso aunque f sea acotada.)
• [Bolzano] f : [a, b] → R continua y con signo distinto en los extremos f (a) y f (b),
entonces ∃c ∈ (a, b) tal que f (c) = 0.
• [Darboux, valores intermedios] f : [a, b] → R continua. Entonces f alcanza todos los
valores comprendidos entre su máximo y su mı́nimo.
20 TEMA 1. FUNCIONES, LÍMITES, CONTINUIDAD Y DERIVABILIDAD
Probar que todo polinomio de grado impar tiene al menos una raı́z real, y que su imagen
es todo R
Probar que la ecuación x−senx = 1 tiene al menos una raı́z real.
Funciones monótonas
Definición 1.31. Sea f : S ⊂ R → R.
Teorema 1.32. Sea f : [a, b] → R una función creciente. Entonces para todo c ∈ (a, b) se
verifica:
∃ lı́m f (x) = f (c− ), ∃ lı́m f (x) = f (c+ ).
x→c− x→c+
Observación 1.33.
Para las funciones continuas y monótonas estrictas, es posible construir la función in-
versa respecto de la composición, y es también continua y estrictamente monótona.
Como ejemplo podemos considerar la función f (x) : R+ → R+ : x 7→ f (x) = xn con
n ∈ N.
Observación 1.35. Para que la derivada exista tiene que existir el lı́mite, es decir, deben
existir los lı́mites laterales y coincidir.
f (2 + h) − f (2) (2 + h)2 − 22 h2 + 4h
f 0 (2) = lı́m = lı́m = lı́m = lı́m (h + 4) = 4.
h→0 h h→0 h h→0 h h→0
|h| − |0| h
f+0 (0) = lı́m = lı́m =1
h→0+ h h→0 h+
pero
|h| − |0| −h
f−0 (0) = lı́m = lı́m = −1.
h
h→0− h→0 − h
Luego existen las derivadas laterales, pero los lı́mites no coinciden. Entonces, la función
valor absoluto no es derivable en a = 0.
h1/3 − 01/3 1
f 0 (0) = lı́m = lı́m 2/3 = +∞
h→0 h h→0 h
luego no se trata de un número real. En este caso, se dice que la función tiene derivada
+∞ en a = 0.
0 0.5 1 1.5 2
x
–1
− − − − es la función y = x2
¦ ¦ ¦ ¦ es la tangente en el punto (1, 1), y = 2x − 1
× × × × es la recta normal en el punto (1, 1), y = (3 − x)/2
Definición 1.39. Se define la recta de pendiente m que pasa por el punto (x0 , y0 ) como:
y − y0 = m (x − x0 )
Dos rectas de pendientes m y m,
e respectivamente, se dice que son perpendiculares cuando
o
forman un ángulo de 90 . Entonces, se puede comprobar que la relación entre sus pendientes
es m
e = −1/m.
Definición 1.40. Si f es derivable en a y f 0 (a) 6= 0, entonces la pendiente de la recta tangente
en el punto (a, f (a)) es f 0 (a), y la pendiente de la recta normal es − f 01(a) .
1. f ± g es derivable en a, siendo:
(λ · f )0 (a) = λ · f 0 (a)
3. f · g es derivable en a, siendo:
f
4. Si g 0 (a) 6= 0, es derivable en a, siendo:
g
µ ¶0
f f 0 (a) · g(a) − f (a) · g 0 (a)
(a) =
g (g(a))2
24 TEMA 1. FUNCIONES, LÍMITES, CONTINUIDAD Y DERIVABILIDAD
Potencia
(xn )0 = nxn−1 (f (x)n )0 = nf (x)n−1 f 0 (x)
Exponenciales
¡ ¢0
(ex )0 = ex ef (x) = ef (x) f 0 (x)
¡ ¢0
(ax )0 = ax · (ln a) af (x) = (ln a)af (x) f 0 (x)
Logarı́tmicas
1 1
(ln x)0 = , x > 0 (ln f (x))0 = f 0 (x)
x f (x)
1 1 1 1 0
(loga (x))0 = (loga f (x))0 = f (x)
ln a x ln a f (x)
Trigonométricas
(senx)0 = cos x (senf (x))0 = f 0 (x) cos f (x)
(cos x)0 = −senx (cos f (x))0 = −f 0 (x) senf (x)
1
(tan x)0 = 1 + (tan x)2 = (tan f (x))0 = [1 + (tan f (x))2 ] f 0 (x)
(cos x)2
−1
(cotanx)0 = −(1 + (cotanx)2 ) = (cotanf (x))0 = − [1 + (cotanf (x))2 ] f 0 (x)
(senx)2
Inversas trigonométricas
1 f 0 (x)
(arcsenx)0 = √ , si |x| < 1 (arcsenf (x))0 = p
1 − x2 1 − f (x)2
−1 −f 0 (x)
(arccosx)0 = √ , si |x| < 1 (arc cos f (x))0 = p
1 − x2 1 − f (x)2
1 f 0 (x)
(arctanx)0 = (arctan f (x))0 =
1 + x2 1 + f (x)2
−1 −f 0 (x)
(arccotanx)0 = (arccotanf (x))0 =
1 + x2 1 + (f (x))2
1.6. APLICACIONES DE LAS DERIVADAS 25
y = (x3 + 2x + 3)4 .
Teorema 1.47. Sea f : S ⊂ R → R, derivable en a ∈ S con f 0 (a) > 0 (ó +∞) (respectiva-
mente, f 0 (a) < 0 (ó −∞)). Entonces, existe un intervalo (a−δ, a+δ) tal que ∀x ∈ (a−δ, a+δ),
x 6= a se tiene:
½
f (x) < f (a) si x < a (respectivamente, f (x) > f (a))
f (x) > f (a) si x > a (respectivamente, f (x) < f (a))
Demostración: Si f 0 (a) > 0, entonces por el Teorema anterior f serı́a localmente estric-
tamente creciente en un intervalo (a − δ, a + δ). Luego no tendrı́a ni máximo ni mı́nimo en
a.
Si f 0 (a) < 0, entonces por el Teorema anterior f serı́a localmente estrictamente decreciente
en un intervalo (a − δ, a + δ). Luego no tendrı́a ni máximo ni mı́nimo en a.
Luego, f 0 (a) = 0.
1
0.8
0.6
0.4
0.2
Se calculan las imágenes de estos puntos y en el (los) punto (puntos) con imagen mayor,
f alcanza el valor máximo absoluto. En el (los) punto (puntos) con imagen menor, f
alcanza el valor mı́nimo absoluto.
2
1.8
1.6
1.4
1.2
1
0.8
0.6
0.4
0.2
0 1 2 3 4
x
Puede comprobarse que f es continua en [0, 4]. Estudiamos cada uno de los puntos:
Comparando las imágenes de los puntos calculados, obtenemos que el valor máximo
absoluto de f en [0, 4] es 2 y se alcanza en x = 4, y el valor mı́nimo absoluto de
f en [0, 4] es 0 y se alcanza en los puntos x = 0 y x = 2.
28 TEMA 1. FUNCIONES, LÍMITES, CONTINUIDAD Y DERIVABILIDAD
1.6
1.4
1.2
0.8
Teorema 1.54. Sea f : [a, b] → R continua en [a, b]. Supongamos que f es derivable en (a, b)
salvo quizás un punto c ∈ (a, b).
Observación 1.55. Este resultado se puede usar para determinar la existencia de extremos
relativos en puntos en los que la función no es derivable.
1.6. APLICACIONES DE LAS DERIVADAS 29
Observación 1.57.
f (x) g(x)
f (x) · g(x) = 1 o f (x) · g(x) = 1 .
g(x) f (x)
5. Si lı́m f (x) = ∞ = lı́m g(x), entonces se puede usar el Teorema para hallar
x→a,∞ x→a,∞
lı́m (f (x) − g(x)) = ”∞ − ∞”. Para ello, se escribe:
x→a,∞
1 1
g(x) − f (x)
f (x) − g(x) = 1 .
f (x)·g(x)
x3 − 3x2 + 4
1. lı́m
x→2 x2 − 4x + 4
µ ¶
x3 − 3x2 + 4 0 3x2 − 6x 6x − 6
lı́m 2 = “ ” = lı́m = lı́m =3
x→2 x − 4x + 4 0 x→2 2x − 4 x→2 2
π
2. lı́m x(arctgx − )
x→+∞ 2
π
π arctgx − 2
lı́m x(arctgx −
) = “∞ · 0” = lı́m 1
x→+∞ 2 x→+∞
x
1
“0” 2 −x2
= = lı́m 1+x1 = lı́m = −1
0 x→+∞ − 2 x→+∞ 1 + x2
x
f 0 (x) f (x)
Observación 1.59. De la no existencia de lı́m 0
no se implica nada sobre lı́m .
g (x) g(x)
–3 –2 –1 0 1 x 2 3
f es convexa en (−1, 2)
–3 –2 –1 0 1 x 2 3
–2
–4
f es cóncava en (−2, 1)
Definición 1.62. Una función derivable en un punto a se dice que es convexa (respectiva-
mente cóncava) en a si existe un entorno de dicho punto en el que la curva se mantiene por
encima (respectivamente, por debajo) de la recta tangente a ella en el punto x = a.
8
6
4
2
–3 –2 –1 0 1 x 2 3
–2
–4
–6
f es convexa en (−2, 2)
32 TEMA 1. FUNCIONES, LÍMITES, CONTINUIDAD Y DERIVABILIDAD
10
8
6
4
2
–3 –2 –1 0 1 x 2 3
–2
–4
f es cóncava en (−2, 2)
Definición 1.63. Decimos que f tiene un punto de inflexión en (a, f (a)) si f cambia en
x = a de cóncava a convexa o de convexa a cóncava.
2
1.8
1.6
1.4
1.2
1
0.8
0.6
0.4
0.2
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
x
f pasa de cóncava a convexa en (1, 1)
1.6. APLICACIONES DE LAS DERIVADAS 33
2
1.8
1.6
1.4
1.2
1
0.8
0.6
0.4
0.2
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
x
f pasa de convexa a cóncava en (1, 1)
Para saber si la curva corta a la ası́ntota horizontal o no, tenemos que resolver la
ecuación:
f (x) = c
f (x)
m = lı́m ∈R
x→±∞ x
los valores de x ∈ Dom(f ) en los que f 0 se anula (si f 0 (x) > 0 la función es
creciente, y si f 0 (x) < 0 la función es decreciente),
los puntos x ∈ Dom(f ) para los que la derivada no existe.
Recordemos que si x = k ∈
/ Dom(f ), entonces puede existir una ası́ntota vertical en
x = k.
8. Máximos y mı́nimos: Con las herramientas estudiadas se buscan los puntos en los
que se alcanzan los máximos y mı́nimos relativos de la función.
1.6. APLICACIONES DE LAS DERIVADAS 35
10. Puntos de inflexión: Son los puntos en los que la función pasa de cóncava a convexa
o de convexa a cóncava. Si f es 2 veces derivable en su dominio, entonces se encuentran
entre los puntos que verifican que f 00 (x) = 0, PERO no todos los que verifican dicha
condición son puntos de inflexión.
36 TEMA 1. FUNCIONES, LÍMITES, CONTINUIDAD Y DERIVABILIDAD
Tema 2
Interpolación polinómica
2.1. Introducción
Teorema 2.1 (Unicidad del polinomio de interpolación). Si p(x) y q(x) son dos polinomios
de grado ≤ n tales que p(xi ) = q(xi ) para i = 0, 1, . . . , n con xi distintos entre si, entonces
p(x) = q(x) ∀x ∈ R.
37
38 TEMA 2. INTERPOLACIÓN POLINÓMICA
donde
x0 y0
y1 −y0
y01 = x1 −x0
y12 −y01
x1 y1 y012 = x2 −x0
y2 −y1 y123 −y012
y12 = x2 −x1 y0123 = x3 −x0 ···
y23 −y12
x2 y2 y123 = x3 −x1
y3 −y2
y23 = x3 −x2
x3 y3
.. .. .. .. ..
. . . . .
xi -2 -1 1 2
Ejemplo 2.3. Hallar el polinomio de grado ≤ 3 que interpola
yi -5 1 1 7
Comenzamos calculando los coeficientes:
xi yi yij yijk yijkl
−2 −5
1−(−5)
−1−(−2) =6
0−6
−1 1 1−(−2) = −2
1−1 2−(−2)
1−(−1) =0 2−(−2) =1
6−0
1 1 2−(−1) =2
7−1
2−1 =6
2 7
ti 0 2 4 6 8
yi 2.64 ×105 1.84×107 2.53×108 2.95×109 4.35×1010
Hallar el polinomio de interpolación que da una aproximación de la ecuación que sigue el
crecimiento celular.
40 TEMA 2. INTERPOLACIÓN POLINÓMICA
Tema 3
Propiedades:
Z Z
1. Si k ∈ R, entonces k f (x) dx = k f (x) dx.
Z Z Z
2. (f (x) ± g(x)) dx = f (x) dx ± g(x) dx
41
42 TEMA 3. LA INTEGRAL. INTEGRACIÓN NUMÉRICA.
Tipo
Z trigonométricas inversas Z
dx f 0 (x)
= arctg x + C dx = arctg(f (x)) + C
1 + x2 1 + f (x)2
Z Z
dx f 0 (x)
√ = arcsen x + C p dx = arcsen (f (x)) + C
1 − x2 1 − f (x)2
1.
Z Z Z
5 5/2 5 1
dx = 3 2 dx = 2
q dx
2 + 3x2 1 + 2x 1 + ( 32 x)2
q Ãr !
Z 3
5 2 5 3
=√ ³q ´2 dx = √ arctg x + C.
6 3 6 2
1+ 2x
2.
Z Z ex µ ¶
ex 2 ex
√ dx = 2 q dx = arcsen + C.
4 − e2x 2
x
1 − ( e2 )2 2
44 TEMA 3. LA INTEGRAL. INTEGRACIÓN NUMÉRICA.
= (x − 1) ex + C.
Z ½ dx
u = arctgx ⇒ du = 1+x2
arctgx dx =
dv = dx ⇒ v = x
Z
x
= x arctgx − dx
1 + x2
1
= x arctgx − 2 ln |x2 + 1| + C.
Ejemplo 3.8.
Z
2x − 3
dx =
2x3 − x2 − x
£ ¤
dado que 2x3 − x2 − x = x (x − 1) (2x + 1)
Z µ ¶
A B C
= + +
x x − 1 2x + 1
1 8
= 3 ln |x| − ln |x − 1| − ln |2x + 1| + C.
3 3
2x − 3 A B C
= + +
2x3− x2 − x x x − 1 2x + 1
A(x − 1)(2x + 1) + Bx(2x + 1) + Cx(x − 1)
=
x(x − 1)(2x + 1)
x2 (2A + 2B + C) + x(−A + B − C) − A
=
2x3 − x2 − x
A = 3
0 = 2A + 2B + C
1
B = −
2 = −A + B − C ⇒ 3
−3 = −A C = − 16
3
También podemos dar valores convenientes a x para obtener un sistema de ecuaciones más
simple para A, B y C.
A1 A2 A3 Ak
+ 2
+ 3
+ ··· +
x − x0 (x − x0 ) (x − x0 ) (x − x0 )k
46 TEMA 3. LA INTEGRAL. INTEGRACIÓN NUMÉRICA.
Ejemplo 3.9.
Z
2x − 3
dx =
x3 − 3x2 + 3x − 1
£ ¤
dado que x3 − 3x2 + 3x − 1 = (x − 1)3
Z µ ¶
A B C
= + + dx
x − 1 (x − 1)2 (x − 1)3
Z µ ¶
0 2 −1
= + + dx
x − 1 (x − 1)2 (x − 1)3
1 1 1
= −2 + +C
(x − 1) 2 (x − 1)2
donde se ha resuelto:
2x − 3 A B C
= + +
x3 − 3x2 + 3x − 1 x − 1 (x − 1)2 (x − 1)3
Ax2 + x(−2A + B) + (A − B + C)
=
(x − 1)3
Hay otras muchas combinaciones, como mezcla de raı́ces reales y complejas (simples y/o
múltiples). Nosotros sólo trataremos aquı́ el caso anterior, y el caso en que la raı́z compleja
1
es imaginaria pura, es decir, del tipo 2 , que ya sabemos es de tipo arcotangente.
a + x2
Consiste en hacer un cambio de variable que transforme la integral en otra que sepamos
calcular. No hay que olvidar, una vez resuelta, deshacer el cambio. Veamos algunos ejemplos:
Z
ex
√ dx = (ex = t ⇒ ex dx = dt)
4 − e2x
Z Z
dt dt
= √ q
=
4 − t22 1 − ( 2t )2
¡ ¢ ¡ x¢
= arcsen 2t + C = arcsen e2 + C.
3.3. LA INTEGRAL DEFINIDA 47
Z √
x
√ dx = (x = t6 ⇒ dx = 6 t5 dt)
3
x+1
Z
t8
= 6 dt (dividimos)
t2 + 1
Z Z
6 4 2 dt
= 6 (t − t + t − 1) dt + 6
1 + t2
µ ¶
t7 t5 t3
= 6 − + − t + 6 arctg(t) + C
7 5 3
à 7 5 1
!
x6 x6 x2 1
= 6 − + − x 6 + 6 arctg(x1/6 ) + C.
7 5 3
Por ejemplo,
Z Z Z
sen3 x t3 dt t3
dx = √ √ = dt
cosx 1−t 2 1 − t2 1 − t2
que se ha transformado en una integral de tipo racional.
También podemos intentar el cambio cosx = t o tgx = t.
Observación 3.10. Estos tres cambios no siempre funcionan, pero son fáciles de hacer.
¡ ¢
El cambio que más habitualmente funciona es tg x2 = t, pero los cálculos son más
complicados.
2. x = tn
Entonces dx = ntn−1 dt. Ver ejemplo anterior.
Observemos que sólo se diferencia de las primitivas o integrales indefinidas en que aparecen
lı́mites de integración a, b.
Si dada la función f (x) conocemos una primitiva de ésta, F (x), entonces se verifica la
Regla de Barrow:
Z b
f (x) dx = F (x)|ba = F (b) − F (a).
a
Z
Observemos que f (x) dx = F (x) + C, es decir, el valor de C no afecta a la aplicación
de la Regla de Barrow.
3.4. APLICACIONES DE LA INTEGRAL DEFINIDA 49
Z ¯2
2
2 x3 ¯¯ 8 1 7
Ejemplo 3.12. x dx = ¯ = − = .
1 3 1 3 3 3
Z b Z b
2. Si k ∈ R, entonces k f (x) dx = k f (x) dx.
a a
Z b Z b Z b
3. (f (x) ± g(x)) dx = f (x) dx ± g(x) dx.
a a a
Z b
5. Si f (x) ≥ 0, ∀x ∈ [a, b], entonces f (x) dx ≥ 0.
a
Z b Z b
6. Si f (x) ≤ g(x) ∀x ∈ [a, b], entonces f (x) dx ≤ g(x) dx.
a a
Si la función tiene cambios de signo en [a, b], hay que separar los intervalos donde f (x)
tiene signo constante y aplicar lo anterior. Por ejemplo, si f (x) ≥ 0 en [a, c] y f (x) ≤ 0
en [c, b], entonces:
Z c Z b
A= f (x) dx − f (x) dx.
a c
50 TEMA 3. LA INTEGRAL. INTEGRACIÓN NUMÉRICA.
Ejemplo 3.13. Calcular el área encerrada por f (x) = x, g(x) = x2 en el intervalo (0, 1) y
en el intervalo (0, 2).
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4
x
El área en el intervalo (0, 2), donde hay dos zonas diferenciadas, una en la que f está sobre
g y otra donde g está sobre f , es:
Z 2 Z 1 Z 2
2
A = |f (x) − g(x)| dx = (x − x ) dx + (x2 − x) dx
0 0 1
µ ¶¯2 µ ¶
1 x3 x2 ¯¯ 1 8 1 1
= + − = + −2− −
6 3 2 ¯1 6 3 3 2
1 8 1 1 1 + 16 − 12 − 2 + 3
= + −2− + = = 1.
6 3 3 2 6
3.4. APLICACIONES DE LA INTEGRAL DEFINIDA 51
Ejemplo: Área del cı́rculo: La curva √ que define el contorno de un cı́rculo de centro
(0, 0) y radio r es x2 + y 2 = r2 , luego y = r2 − x2 .
1
0.8
0.6
y
0.4
0.2
Gracias a la simetrı́a de la figura, para calcular el área lo hacemos para x ∈ [0, r] y multipli-
camos por 4:
Z ³ x ´2 Z r
r p r
A = 4 r2 − x2 dx1− = 4r dx
0 0 r
( x
= sent
haciendo el cambio de variable r
dx = r cost
Z π/2
= 4r2 cos2 t dt
0
1 + cos(2t)
usando la fórmula trigonométrica cos2 (t) =
2
Z π/2 µ ¶ µ ¶¯π/2
1 + cos(2t) t sen(2t) ¯¯
= 4r2 dt = 4r2 + ¯ = π r2 .
0 2 2 4 0
La integral definida nos va a servir para calcular el valor medio o el promedio de una
magnitud. Por ejemplo, pensemos que se ha medido la concentración de nitrógeno en el suelo
cada metro, en una sección transversal de tundra húmeda y se han obtenido los datos
52 TEMA 3. LA INTEGRAL. INTEGRACIÓN NUMÉRICA.
Distancia (m) 1 2 3 4 5
Concentración (g/m3 ) 589.3 602.7 618.5 667.2 641.2
Distancia (m) 6 7 8 9 10
Concentración (g/m3 ) 658.3 672.8 661.2 652.3 669.8
y cada parte tiene una longitud de 10−0 n que denotaremos por h. Por tanto, h = 10
n y de
aquı́ n1 = 10
1
h, que sustituyéndolo en (3.1) da
n
1 X
c̄ = c(xk )h.
10
k=1
Cuando n tiende a ser muy grande, matemáticamente decimos que n tiende a infinito, la suma
Xn Z 10
c(xk )h se parece mucho a c(x)dx. Por tanto,
k=1 0
Z 10
1
c̄ = c(x)dx,
10 0
es decir, la integral de c(x) en el intervalo [0, 10] dividida por la longitud de este intervalo. A
esto se le llama valor medio de c.
La concentración media entre a y b es
Z b
1
c̄ = c(x)dx.
b−a a
Ejercicio 3.14. La temperatura de una cámara varı́a en un periodo de 24 horas de la siguiente
forma:
π
T (t) = 68 + sen( t)
12
para 0 ≤ t ≤ 24.
2
porque no se puede obtener de forma explı́cita una primitiva de e−x . Otras veces, no se
tiene una expresión de la función que se quiere integrar sino sque ólo se conocen valores de la
función en ciertos puntos.
En estas situaciones son necesarios métodos numéricos que proporcionen una aproximación
de la integral.
La idea que subyace en el método que explicaremos a continuación consiste en aproximar
el integrando f (x) por una función más sencilla, haciendo interpolación, de manera que sólo
sea necesario conocer los valores de f en determinados puntos. Además, la nueva función
(interpoladora) se integrará de una manera muy simple porque será un polinomio, y en nuestro
caso, un polinomio de grado 1. Por consiguiente, en este apartado deduciremos una fórmula
Z b
que nos da de manera aproximada el valor de f (x)dx.
a
En todo lo que sigue supondremos que el intervalo [a, b] se divide en n partes iguales, de modo
que la longitud de cada subintervalo es h = b−a n , y que f es una función continua en [a, b] de
la que conocemos sus valores en los puntos x0 , x1 , . . . , xn , que se han obtenido de dividir el
intervalo [a, b] en n partes, es decir, x0 = a, x1 = x0 + h, x2 = x1 + h,...,xn = b. Y conocemos
f (x0 ), . . . , f (xn ).
En cada subintervalo [xi , xi+1 ] vamos a aproximar f (x) por el polinomio de interpolación de
grado menor o igual que 1 en [xi , xi+1 ], es decir,
f (xi+1 ) − f (xi )
pi (x) = f (xi ) + (x − xi ).
h
Z xi+1 Z xi+1
Entonces, f (x)dx se aproxima por pi (x)dx.
xi xi
Z xi+1 Z xi+1
f (xi+1 ) − f (xi ) f (xi+1 ) + f (xi )
f (x)dx ∼ pi (x)dx = f (xi )h + h= h.
xi xi 2 2
Ahora expresamos la integral en [a, b] de f (x) como suma de integrales en los subintervalos
[xi , xi+1 ] (por la propiedad que ya vimos de las integrales definidas) y sustituimos cada una
de estas integrales por su aproximación:
Z b X Z xi+1
n−1
f (x)dx = f (x)dx ∼
a i=0 xi
n−1
X f (xi+1 ) + f (xi )
∼ h=
2
i=0
54 TEMA 3. LA INTEGRAL. INTEGRACIÓN NUMÉRICA.
h
= [f (a) + 2f (x1 ) + · · · + 2f (xn−1 ) + f (b)].
2
Es decir, h/2 por la suma de f en los dos extremos a y b y la suma del doble de f en los
puntos intermedios x1 , · · · , xn−1 . A esta fórmula se la conoce como fórmula de los trapecios
Z b
porque, si f fuese positiva, f (x)dx representarı́a el área limitada por f entre a y b, y la
a
fórmula de los trapecios no es más que la suma de las áreas de todos los trapecios que se
forman uniendo por rectas los puntos (xi , f (xi )), como muestra la figura.
El error cometido al aproximar la integral por la fórmula de los trapecios, denotada por Tn ,
viene dado por la estimación siguiente
¯Z b ¯
¯ ¯ (b − a)3 b−a 2
¯ f (x)dx − T ¯ ≤ máx |f 00
(x)| = máx |f 00 (x)| h .
¯ n¯
x∈[a,b] 12n 2 x∈[a,b] 12
a
si llamamos h = b−a2 , y la fórmula de Simpson generalizada para n puntos (n debe ser par)
Z b
h
f (x)dx ∼ [f (a) + 4(f (x1 ) + f (x3 ) + · · · + f (xn−1 )) + 2(f (x2 ) + · · · + f (xn−2 )) + f (b)]
a 3
El error cometido por la fórmula generalizada de Simpson, denotada por Sn , es ahora del
orden de h4 ¯Z b ¯
¯ ¯ 5
¯ f (x)dx − S ¯ ≤ máx |f IV ) (x)| (b − a)
¯ n ¯ x∈[a,b] 2880( n )4
a 4
y la fórmula es exacta hasta polinomios de grado tres.
Tema 4
Programación lineal
Mientras que para funciones reales de variable real la derivación ha permitido resolver el
problema de optimalidad en su conjunto, en este tema, la programación lineal resuelve prob-
lemas de optimización en varias variables para funciones objetivos de un tipo muy particular:
funciones lineales.
Existen multitud de aplicaciones relativas a la programación lineal, y métodos de resolu-
ción mucho más complejos del único que aquı́ describimos, que será el método geométrico.
Ilustraremos con ejemplos el tipo de marco que se puede recoger y tratar a través de este
método.
Definición 4.1. Los problemas de óptimos para expresiones lineales donde las re-
stricciones son desigualdades dadas a partir de más expresiones lineales (inecuaciones)
conforman la PROGRAMACIÓN LINEAL.
Ejemplo 4.2. Un ejemplo de tal tipo de problemas es el siguiente:
55
56 TEMA 4. PROGRAMACIÓN LINEAL
En una granja se da una dieta “para engordar” con una composición mı́nima de 15
unidades de una sustancia A y otras 15 de una sustancia B. En el mercado sólo se encuentran
dos clases de compuestos: el tipo X con una composición de una unidad de A y cinco de B, y
el tipo Y, con una composición de cinco unidades de A y una de B. El precio del tipo X es de
10 euros y el del tipo Y es de 30 euros. Se pregunta: ¿Qué cantidades se han de comprar de
cada tipo para cubrir las necesidades con un coste mı́nimo?
El planteamiento matemático es: Hallar x (latas del tipo X) e y (latas del tipo Y) que
resuelven el problema
mı́n(10x + 30y)
D
siendo
¯ x ≥ 0, y ≥ 0,
¯
D = (x, y) ∈ R2 ¯ x + 5y ≥ 15, .
5x + y ≥ 15
-Problema del transporte: organizar reparto de mercancı́as con coste mı́nimo de tiempo,
dinero o riesgo
Para atender el suministro diario de gas a tres ciudades C1 , C2 y C3 una empresa tiene
destinadas dos fábricas F1 y F2 que producen 20 y 30 m3 respectivamente. Las necesidades de
las tres ciudades son: 20, 18 y 12 m3 respectivamente. Si los costes de transporte por tonelada
de las industrias a las ciudades son, en cientos de euros, los indicados en la tabla adjunta,
planificar el reparto óptimo para que dicho coste sea mı́nimo.
F1 4 3 1
F2 2 2 1
C1 C2 C3
-Problema de la ruta más corta (o del viajante): ordenar etapas de un viaje con el propósito
de minimizar el recorrido
4.3. RESOLUCIÓN POR EL MÉTODO GEOMÉTRICO 57
Usaremos el método geométrico, que implica dominar los conceptos de función objetivo,
restricciones, región factible, rectas de nivel ası́ como la resolución de sistemas de
ecuaciones e inecuaciones lineales.
D = {(x, y) ∈ R2 | 0 ≤ x, 0 ≤ y, 3x + y ≤ 4, 2x + 3y ≤ 6}.
4.5
3.5
2.5
1.5
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Ahora marcamos la intersección de los semiplanos (región coloreada) que verifican todas las
desigualdades que definen D.
4.5
3.5
2.5
1.5
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
son rectas, y todas son paralelas entre si cuando la constante cambia (dicha constante es el
valor de la función objetivo a lo largo de todos los puntos (x, y) de esa recta),
Dicha función, “trasladada paralelamente” sobre la región genera los valores posibles de la
constante. Ası́, los valores mayores y menores se obtienen en los extremos, por ejemplo, el
máximo saldrı́a de...
4.5
3.5
2.5
1.5
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Paso 1: calculamos los vértices Esto es, los puntos de la intersección dos a dos de las
“ecuaciones asociadas” d1 x + d2 y = d3 a las inecuaciones d1 x + d2 y < (>, ≤, ≥) d3 que definen
el dominio (ojo: asegurarse que están en la región factible).
En el caso de la región D del dibujo anterior se trata de los puntos (0, 0), (0, 2), (4/3, 0),
(6/7, 10/7).
E = {(x, y) ∈ R2 | 0 ≤ x, 0 ≤ y, 3x + y ≥ 4, 2x + 3y ≥ 6}.
el primero tiene solución mientras que el segundo no posee solución: f tiene mı́nimo sobre E
pero no máximo, de hecho veremos que sup f (x, y) = +∞.
E
Ahora se puede comprobar que, si bien lı́neas que minimizan el valor constante= x + y
tienen un “tope” pues van descendiendo y el último punto factible es el (6/7, 10/7), en sentido
ascendente (generando valores cada vez mayores constante= x + y) no tiene máximo finito:
5
4.5
3.5
2.5
1.5
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Existen puntos (x, y) ∈ E tan grandes como queramos, haciendo, como anunciamos, que
sup f (x, y) = +∞.
E