Matemáticas I

José Leandro de Marı́a y José Luis Estévez Balea

1
Índice general

Índice general 2

1 Fundamentos 5
1.1. Álgebra básica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.2. Intervalos y valores absolutos . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.3. Leyes de los exponentes . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.4. Cı́rculos y parábolas. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14

2 Álgebra 17
2.1. Sistemas de ecuaciones lineales . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.2. Matrices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.3. Sistemas en forma matricial . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.4. Determinantes. Regla de Cramer . . . . . . . . . . . . . . . . . . . . . . . 26
2.5. Teorema de Rouché-Frobenius . . . . . . . . . . . . . . . . . . . . . . . . . 31
2.6. Espacios vectoriales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
2.7. Producto escalar. Ortogonalidad . . . . . . . . . . . . . . . . . . . . . . . . 39
2.8. Aplicaciones lineales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
2.9. Isometrı́as en Rn . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
2.10. Ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
2.11. Diagonalización de aplicaciones lineales simétricas . . . . . . . . . . . . . . 63

3 Geometrı́a 71
3.1. Espacio afı́n . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
3.2. Problemas de incidencia . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
3.3. El espacio euclidiano tridimensional . . . . . . . . . . . . . . . . . . . . . . 81
3.4. Producto escalar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
3.5. Espacio métrico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
3.6. El espacio euclidiano . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
3.7. Ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88

4 Análisis 95
4.1. Funciones elementales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
4.2. Sucesiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
4.3. Seno y coseno de un ángulo . . . . . . . . . . . . . . . . . . . . . . . . . . 101
4.4. El concepto de lı́mite de una función en un punto . . . . . . . . . . . . . . 102
4.5. Continuidad de una función . . . . . . . . . . . . . . . . . . . . . . . . . . 105
4.6. Lı́mites infinitos y lı́mites en el infinito . . . . . . . . . . . . . . . . . . . . 105
4.7. La derivada . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107

2
ÍNDICE GENERAL 3

4.8. Teoremas de Rolle, Cauchy y del Valor Medio . . . . . . . . . . . . . . . . 110

4.9. La derivada segunda. Concavidad y convexidad. . . . . . . . . . . . . . . . 113
4.10. Teorema de Taylor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 113
4.11. Aplicaciones de la fórmula de Taylor . . . . . . . . . . . . . . . . . . . . . 118
4.12. Cálculo de primitivas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 120
4.13. La integral de Riemann . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 125

Índice 129
Capı́tulo 1

Fundamentos

Los números reales están ordenados como los puntos de una recta.

1.1. Álgebra básica

A los números 1, 2, 3, · · · que aparecen de forma natural en el proceso de contar, les
llamamos números naturales. La introducción del 0, ası́ como la de los números negativos
−1, −2, −3, · · · aparece como consecuencia de la operación de substracción. Del mismo
modo la operación de división nos lleva a considerar los números racionales pq donde p y
q son números enteros.

Ejemplo 1 Determinar si los siguientes números son naturales, enteros o racionales.

4+5 4+5
a) 0 b) 3 − 2
c) 7 − 4 d) −3

Solución: a) 0 no es un número natural. b) 3 − (4 + 5)/2 = 3 − 9/2 = −3/2 no es

un número natural ni un número entero , pero es un número racional . c) 7 − 4 = 3 es
un número natural, un entero y un número racional. d) (4+5)/(-3)=-3 no es un número
natural, pero es un entero y un número racional.

Los Griegos ya sabı́an que algunas lı́neas de figuras geométricas muy simples tenı́an lon-
gitudes que no se correspondı́an con cocientes √ de números enteros. La longitud de la
diagonal de un cuadrado cuyos lados miden 1 es 2, que no puede expresarse en la forma
p/q, con p y q números enteros. Lo mismo ocurre con la circunferencia cuyo diámetro es
1: su longitud π no es un número racional. Los números que no se pueden expresar como
racionales se llaman irracionales. Estos números junto con los racionales constituyen el
conjunto de los números reales.

5
6 CAPÍTULO 1. FUNDAMENTOS

Dos identidades algebraicas que nos serán de utilidad son

(a + b)(a − b) = a2 − b2 y (a + b)2 = (a + b)2 = a2 + 2ab + b2

Ejemplo 2 Simplificar la expresión (a + b)(a − b) + b2 .

Solución: Como (a + b)(a − b) = a2 − b2 , tenemos (a + b)(a − b) = a2 − b2 + b2 = a2 .

Ejemplo 3 Desarrollar la expresión (a + b)3 .

Solución: Tenemos (a + b)2 = a2 + 2ab + b2 . Entonces

(a + b)3 = (a + b)2 (a + b)
= (a2 + 2ab + b2 )(a + b)
= (a2 + 2ab + b2 )a + (a2 + 2ab + b2 )b
= a3 + 2a2 b + ab2 + a2 b + 2ab2 + b3
= a3 + 3a2 b + 3ab2 + b3

Ejemplo 4 Resolver la ecuación x2 − 5x + 3 = 0 completando el cuadrado.

Solución: Transformamos la ecuación añadiendo y restando ( 25 )2 :

x2 − 5x + ( 25 )2 − ( 52 )2 + 3 = 0
(x − 52 )2 − 134
= 0
(x − 25 )2 = 13
4√
x − 52 = ± 213√
x = 52 ± 213

Otra forma de completar el cuadrado es haciendo

x2 − 5x + 3 = (x + p)2 + q

y entonces
x2 − 5x + 3 = x2 + 2px + p2 + q.
5
Igualando coeficientes, tenemos que p = 2
y q = 3 − p2 = − 13
4
.

Ecuación de segundo grado.

1.2. INTERVALOS Y VALORES ABSOLUTOS 7

Aplicando el método de completar el cuadrado a la ecuación general de segundo grado

ax2 + bx + c = 0 con a 6= 0 tenemos la fórmula conocida
√
−b ± b2 − 4ac
x=
2a

y dependiendo de si el discriminante (b2 − 4ac) es mayor, igual o menor que 0, tendremos

dos, una o ninguna solución respectivamente.

Ejercicio Comprobar las siguientes igualdades, que pueden considerarse una generaliza-
ción de que la suma por diferencia es la diferencia de cuadrados, y que pueden ser de
utilidad en el cálculo algebraico:

a3 − b3 = (a − b)(a2 + ab + b2 )
a4 − b4 = (a − b)(a3 + a2 b + ab2 + b3 )
en general se tiene:

an − bn = (a − b)(an−1 + an−2 b + · · · + abn−2 + bn−1 )

1.2. Intervalos y valores absolutos

La notación para los distintos tipos de intervalos es la siguiente

(a, b) es el conjunto de todos los x tales que a < x < b.

[a, b] es el conjunto de todos los x tales que a 6 x 6 b.

[a, b) es el conjunto de todos los x tales que a 6 x < b.

(a, b] es el conjunto de todos los x tales que a < x 6 b.

[a, ∞) es el conjunto de todos los x tales que a 6 x.

(a, ∞) es el conjunto de todos los x tales que a < x.

(−∞, b] es el conjunto de todos los x tales que x 6 b.

(−∞, b) es el conjunto de todos los x tales que x < b.

(−∞, ∞) es el conjunto de todos los números reales.

Valor absoluto El valor absoluto |x| de un número real x es igual a x si x > 0 y −x si

x < 0.
8 CAPÍTULO 1. FUNDAMENTOS

Propiedades del valor absoluto

Para dos números reales x e y se verifica:

1. |x + y| 6 |x| + |y|

2. |xy| = |x||y|
√
3. |x| = x2

Ejemplo 5 Mostrar que no siempre se tiene |x + y| = |x| + |y|.

Solución: Sean x = 5 e y = −2. Entonces se tiene que |x + y| = |5 − 2| = 3, y por otra

parte |x| + |y| = 5 + 2 = 7. De hecho la igualdad nunca se tiene cuando x e y tienen
distinto signo.

√
Ejemplo 6 Probar que |x| = x2 .

Solución: Para cualquier número x se tiene que (−x)2 = x2 ,pentonces |x|2 = x2 inde-
pendientemente del signo de x. Como |x| > 0 se tiene que |x|2 = |x| y de aquı́ la
igualdad.

Ejercicio Encontrar todos los números reales tales que |x| = 5.

Ejercicio Describir, utilizando el lenguaje de intervalos, el conjunto de los x que verifican

las siguientes desigualdades:

1. x2 + 3x > 0

2. x2 − 2x < 0

3. x2 − x − 2 > 0

Solución:

1. La condición la podemos escribir

x(x + 3) > 0
1.2. INTERVALOS Y VALORES ABSOLUTOS 9

y esta desigualdad se tiene cuando los dos factores tienen el mismo signo. Entonces tiene
que cumplirse una de las dos condiciones siguientes: a) x > 0 y x > −3 o bien b) x < 0 y
x < −3.

La condición a) se cumple en la intersección de los intervalos (0, +∞) ∩ (−3, ∞) = (0, ∞)

y la segunda en (−∞, −3) ∩ (−∞, 0) = (−∞, −3). Resulta entonces que la inecuación del
enunciado se cumple en la unión de estos intervalos:

I = (−∞, −3) ∪ (0, ∞)

2. De manera análoga se tiene en este caso que el conjunto de valores de x que verifican
la inecuación del enunciado es el conjunto

J = (0, 2)

3. En este caso, completando el cuadrado tenemos:

x2 − x − 2 > 0
1 1
x2 − x + − − 2 > 0
4 4
2 1 9
x −x+ − > 0
4 4
 2
1 9
x− >
2 4
1 3
x− >
2 2

y esta desigualdad se tiene, bien cuando x − 1/2 > 3/2, que es decir x > 2 o bien cuando
−x + 1/2 > 3/2, que es cuando −x > 1 ⇔ x < −1. En definitiva tenemos el conjunto

K = (−∞, −1) ∪ (2, ∞)

Ejercicio Expresar las siguientes desigualdades en términos de: “x pertenece al interva-

lo...”:

1. 2 < x > 6

2. |x| < 3

3. |x − 3| > 5

4. x2 − 3x + 2 > 0
10 CAPÍTULO 1. FUNDAMENTOS

Ejercicio Escribir las siguientes expresiones en terminos de desigualdades que implican

valores absolutos:

1. x ∈ (−3, 3)

2. −x ∈ (−3, 3)

3. x ∈ [−6, 14]

1.3. Leyes de los exponentes

La expresión an , donde a es un número real llamado base y n es un número natural
llamado exponente, se define como el producto de a por si mismo n veces.

an = a · a · · · a (n veces)

Leyes de los exponentes: potencias enteras

1. bn bm = bn+m

2. (bn )m = bnm

3. (bc)n = bn cn

Ejemplo 7 Simplificar las expresiones:

(3 · 5)10 + 39
a) 210 · 510 b)
39

Solución:

a) 210 · 510 = (2 · 5)10 = 1010 .

b)
(3 · 5)10 + 39 310 · 510 + 39 39 · 3 · 510 + 39
= = = 3 · 510 + 1.
39 39 39

La primera de las leyes anteriores es básica para extender la exponenciación para permitir
exponentes negativos. Si b 6= 0, entonces bn 6= 0 y de b0 bn = bn se deduce que b0 = 1. Hay
que notar que la expresión 00 no tiene sentido.
1.3. LEYES DE LOS EXPONENTES 11

Exponentes negativos
Si b es un número real y n un entero positivo entonces se define

1
b−n = .
bn
Las leyes de los exponentes permanecen válidas para n y m positivos, negativos o cero.

Ejemplo 7 Simplificar las expresiones:

(2 · 3)−2 · 4
a) .
(1/3)2

(2 · 3−2 ) · 4
b) .
(1/3)2

c) [(8/3)2 − (3/8)3 ][(8/3)−2 − (3/8)−3 ].

Solución:

(2 · 3)−2 · 4 4 4 · 32
a) = = =1
(1/3)2 (2 · 3)2 · (1/3)2 22 · 32

(2 · 3−2 ) · 4 2·4 2·4

b) 2
= 2 2
= =8
(1/3) 3 · (1/3) 1

c)

 2  −2  2  −3  3  −2  3  −3

8 8 8 3 3 8 3 3
− − + =
3 3 3 8 8 3 8 8
 2  3  3  2
8 8 3 3
=1− − +1
3 3 8 8
 5  5
8 3
=2− − .
3 8
12 CAPÍTULO 1. FUNDAMENTOS

Potencias racionales
Las potencias racionales están definidas del siguiente modo:

bn = b · · · b (n veces); b0 = 1
b−n = 1/b
√
n
n
b1/n = b si b > 0 y n es un número natural.
bm/n = m 1/n
(b )
Si b, c > 0 y p, q son racionales, entonces:

1. bp+q = bp + bq
2. bpq = (bp )q
3. (bc)p = bp cp
4. bp < bq si b > 1 y p < q; bp > bq si b < 1 y p < q

Ejemplo 8 Hallar el valor de: 27−4/3 .

√
Solución: 27−4/3 = 1/274/3 = 1/( 3 27)4 = 1/34 = 1/81.

Ejemplo 9 Si x > 0, simplificar las expresiones: [x3/4 x−4/3 ]5/3 y (x4/3 )3/2 /x1/3

Solución:

√
[x3/4 x−4/3 ]5/3 = (x3/4−4/3 )5/3 = (x−7/12 )5/3 = x−35/36 = 1/( 36 x)35 .

√
(x4/3 )3/2 /x1/3 = x4/3·3/2−1/3 = x5/3 = ( 3 x)5 .

Ejemplo 10 Se definió am/n como (am )1/n , para a > 0. Veamos que am/n es también
(a1/n )m .

Solución: Para ello hay que demostar que (a1/n )m es la raı́z n−sima de am . Tenemos que
[(a1/n )m ]n = (a1/n )mn = (a1/n )nm = [(a1/n )n ]m = am .

Ejercicios Simplificar las expresiones:

x7/2 (x−5/2 + 3x1/2 + 3x3/2 )

sp
5 x5 y 3
p3
x4 y 4
1.3. LEYES DE LOS EXPONENTES 13

Lı́neas rectas en el plano

Consideramos la representación del plano como dos lı́neas rectas perpendiculares, ejes,
y en cada una de ellas situamos los números reales, ası́ cada punto P del plano queda
determinado por un par de números reales (p1 , p2 ),coordenadas , que resultan al trazar
dos perpendiculares a los ejes.

Distancia entre dos puntos del plano

Sean los puntos P1 de coordemadas (x1 , y1 ) y P2 de coordenadas (x2 , y2 ). La distancia de

P1 a P2 que se denota d(P1 , P2 ) es:

p
(x1 − x2 )2 + (y1 − y2 )2

Pendiente determinada por dos puntos del plano Sea r la recta que pasa por los
puntos (x1 , y1 ) y (x2 , y2 ) donde x1 6= x2 , la pendiente de r es:
y2 − y1
m=
x2 − x1

Ecuación de la recta La ecuación de la recta que pasa por los puntos (x1 , y1 ) y (x2 , y2 )
donde x1 6= x2 , es

y = y1 + m(x − x1 )

Observación: Si los puntos considerados antes verifican x1 = x2 e y1 6= y2 entonces la

ecuación de la recta que pasa por ellos es y = x1 .

Ejemplo 11 Encontrar la ecuación de la recta que pasa por (1,2) con inclinación 3.

Solución: Sustituyendo los valores en la ecuación tenemos

y = 2 + 3(x − 1).

Ejemplo 12 Sea r la recta que pasa por los puntos (3, 1) y (2, −2). Hallar el punto donde
r corta al eje x.

Solución: La ecuación de la recta es:

 
−2 − 1
y =1+ (x − 3) = 3x − 8
2−3

La recta interseca al eje x en el punto de ella donde la coordenda y vale 0, es decir

0 = 3x − 8 de donde tenemos x = 8/3.
14 CAPÍTULO 1. FUNDAMENTOS

Ejemplo 13 Hallar la inclinación de la recta 2y − 5x + 8 = 0.

Solución: Transformamos la ecuación en otra equivalente pero de forma más familiar:

2y − 5x + 8 = 0
2y = 5x − 8
5
y = x−4
2

donde se ve que la inclinación es 52 .

1.4. Cı́rculos y parábolas.

Consideraremos en esta sección dos figuras geométricas del plano que se pueden describir
por medio de ecuaciones muy sencillas. El cı́rculo de radio r y centro (a, b) es el lugar de
los puntos (x, y) del plano cuya distancia a (a, b) es r. En ecuaciones se escribirı́a:
p
(x − a)2 + (y − b)2 = r

o lo que es lo mismo

(x − a)2 + (y − b)2 = r2

Ejemplo 14 Hallar la ecuación del cı́rculo de centro (2, 1) y radio 7.

Solución: De la ecuación tenemos (x − 2)2 + (y − 1)2 = 49 o simplificando x2 + y 2 − 4x −

2y − 44 = 0.

Ejemplo 15 Mostrar que la ecuación x2 + y 2 − 4x − 14y + 41 = 0 es de un cı́rculo.

Encontrar su centro y radio.

Solución: Completando cuadrados tenemos:

0 = x2 + y 2 − 4x − 14y + 41 = 0
= (x2 − 4x + 4) + (y 2 − 14y + 49) − 4 − 49 + 41
= (x − 2)2 + (y − 7)2 − 12

√
donde se ve que es la ecuación de un cı́rculo de centro (2, 7) y radio 12.

La otra figura geométrica que analizaremos es la parábola. La ecuación más sencilla de

esta figura es y = x2 , y en general cualquier ecuación de la forma y = ax2 donde a 6= 0
1.4. CÍRCULOS Y PARÁBOLAS. 15

describe todas las posibles parábolas con vértice en (0, 0) y eje de simetrı́a perpendicular
al eje x.

La construcción geométrica de la parábola es la siguiente: sea una recta d contenida en

el plano y un punto F no contenido en esa recta. El conjunto de los puntos que distan lo
mismo de d que de F constituyen una parábola de foco F y directriz d. Tomemos como
directriz la recta y = −a/2 y foco el punto F (0, a/2). La distancia entre F y d es a y expre-
sando la condición de equidistancia analı́ticamente dist ((x, y), d) = dist ((x, y), (0, a/2)),
tenemos:

r
a 2
 a 2
y+ = x + y−
2 2
 a 2  a 2
y+ = x2 + y −
2 2
x2
y =
2a

Una propiedad geométrica de la parábola es la siguiente: si consideramos, por ejemplo, la

parábola y = x2 , todas las semirectas perpendiculares al eje x y contenidas en el semiplano
y > 0 se reflejan en la parábola y las rectas reflejadas pasan todas por el foco F .

Ejercicio En la figura se muestra la parábola de ecuación y = x2 . Cuál es la ecuación de

la directriz que aparece dibujada en azul?.
Capı́tulo 2

Álgebra

2.1. Sistemas de ecuaciones lineales

Ecuaciones lineales
Una ecuación lineal con n incógnitas es una expresión de la forma t

a1 x 1 + a2 x 2 + · · · + an x n = b

donde x1 , x2 , · · · , xn son las incógnitas y a1 , a2 , · · · , an y b son números complejos (no-

sotros trabajaremos con números reales). Los números a1 , a2 , · · · , an se llaman los coefi-
cientes de la ecuación y b se llama constante de la ecuación.

Soluciones
Una solución de la ecuación anterior es un conjunto de valores reales, también llamados
escalares k1 , k2 , · · · , kn para las incógnitas x1 , x2 , · · · , xn , respectivamente, que satisfacen
la ecuación; es decir, que se verifica la igualdad

a1 k1 + a2 k2 + · · · + an kn = b

Ejemplo Consideremos la ecuación lineal ax = b. entonces tenemos que

1. Si a 6= 0, x = b/a es la única solución de ax = b.

2. Si a = 0 y b 6= 0, la ecuación ax = b no tiene solución.

3. Si a = 0 y b = 0 todo número k es solución de la ecuación ax = b.

Ecuaciones lineales con dos incógnitas

Una ecuación lineal con dos incógnitas la podemos escribir de forma general como

17
18 CAPÍTULO 2. ÁLGEBRA

ax + by = c
donde x e y son las incógnitas. Vamos a suponer que la ecuación es no degenerada, es
decir, que a y b no son ambos nulos. Entonces cada solución de la ecuación es un par de
números reales, s = (u1 , u2 ), que pueden hallarse despejando una de las incógnitas.

Si representamos las soluciones de la ecuación anterior como puntos en el plano de coor-

dendas cartesianas R2 , tenemos que el conjunto de soluciones es una lı́nea recta.

Sistemas de ecuaciones lineales

Un sistema de m ecuaciones lineales en n incógnitas es una expresión de la forma

L1 : a11 x1 + a12 x2 + · · · + a1n xn = b1

L2 : a21 x1 + a22 x2 + · · · + a2n xn = b1
··· ··· ···
Lm : am1 x1 + am2 x2 + · · · + amn xn = bm

Solución de un sistema de ecuaciones

Un conjunto de números x1 = λ1 , · · · , xn = λn que verifican las ecuaciones L1 , · · · , Lm
del sistema anterior, de llama solución particular del sistema. El conjunto S de todas
las soluciones particulares del sistema se llama solución general o conjunto solución. Si
S = ∅ se dice que el sistema es incompatible y en caso contrario se dice que el sistema es
compatible.

Sistemas equivalentes
Dos sistemas de ecuaciones lineales en las mismas incógnitas son equivalentes cuando el
conjunto de solución es el mismo.

Operaciones elementales
Llamamos operaciones elementales en un sistema de ecuaciones a las siguientes:

E1 Intercambiar las ecuaciones Li y Lj :

Li ↔ Lj

E2 Multiplicar la ecuación Li por un escalar distinto de cero:

Li ↔ kLi , k 6= 0

E3 Sustituir la ecuación Li por kLj + Li , siendo k un escalar:

Li ↔ kLj + Li
2.1. SISTEMAS DE ECUACIONES LINEALES 19

La razón de hablar de operaciones elementales en este momento está en el siguiente

resultado

Teorema 2.1. Sean S1 y S2 dos sistemas de ecuaciones de tal modo que uno se obtiene
del otro tras un número finito de operaciones elementales. Entonces S1 y S2 tienen la
misma solución general.

Este teorema nos proporciona un método de resolución de un sistema: si a partir de un

sistema, obtenemos por medio de operaciones elementales otro de forma más sencilla,
resolvemos esta última forma y ya tenemos la solución del sistema.

Ejemplo Considerermos el sistema

L1 : x1 + 2x2 − 2x3 = 3
L2 : 3x1 + 9x2 − 12x3 = −15
L3 : 10x1 − x2 + 4x3 = 12

Sustituimos primero L2 por −3L1 + L2 y tenemos

L1 : x1 + 2x2 − 2x3 = 3
L2 : 3x2 − 6x3 = −24
L3 : 10x1 − x2 + 4x3 = 12

donde seguimos llamando a las nuevas ecuaciones L1 , L2 y L3 . Ahora sustituimos L3 por

−10L1 + L3 y resulta

L1 : x1 + 2x2 − 2x3 = 3
L2 : 3x2 − 6x3 = −24
L3 : −21x2 + 24x3 = −18

y finalmente se sustituye L3 por 7L2 + L3 para obtener

L1 : x1 + 2x2 − 2x3 = 3
L2 : 3x2 − 6x3 = −24
L3 : −18x3 = −186
20 CAPÍTULO 2. ÁLGEBRA

Sistemas en forma triangular y escalonada

Un sistema de ecuaciones como el anterior se llama sistema en forma triangular. En
general, un sistema con el mismo número de ecuaciones que de incógnitas donde la primera
incógnita de la k-ésima ecuación es xk , se llama sistema en forma triangular.

Los sistemas en forma triangular tienen solución única. Ésta se obtiene partiendo de la
última ecuación, que nos da el valor de xn , luego sustituyendo este valor en la ecuación
Ln−1 de donde obtendremos el valor de xn−1 y ası́ sucesivamente. Por ejemplo, para
resolver el sistema anterior tendrı́amos enn primer lugar

186 31
x3 = =
18 3
de aquı́, sustituyendo en L2 se tiene

31
3x2 − 6 · = −24
3
y entonces

64 38
3x2 = − 24 = 38 ⇒ x2 =
3 3
ya sólo queda sustituir en L1 los valores de x3 y x2 y tendremos

38 31 76 62 5
x1 + 2 · −2· = 3 ⇒ x1 = 3 − + =−
3 3 3 3 3
y ya tenemos la solución del sistema.

De forma más general podemos considerar sistemas en forma escalonada. Un sistema en

forma escalonada es aquel donde no existen ecuaciones degeneradas (identicamente nulas)
y la primera incógnita de la ecuación Li se encuentra a la derecha de la primera incógnita
de la ecuación Li−1 .

La forma general de un sistema en forma escalonada es

L1 : a11 x1 + a12 x2 + a13 x3 + · · · + a1n xn = b1

L2 : a2j2 xj2 + a2j2 +1 xj2 +1 + · · · + a2n xn = b2
··· ··· ··· ···
Lm : amjm xjm + amjm +1 xjm +1 + · · · + amn xn = bm
donde 1 < j2 < · · · < jm y a11 6= 0, a2j2 6= 0, · · · , amjm 6= 0

Ejercicio 2.1. Hallar la solución del sistema

x1 + 5x2 − x3 = 2
−12x3 = −1
2.2. MATRICES 21

2.2. Matrices
Una matriz A = (aij ) es una expresión de la forma

 
a11 a12 ··· a1n
 a21 a22 ··· a2n 
A= 
 ··· ··· ··· 
am1 am2 ··· amn
donde los números aij pertenecen a un cuerpo, generalmente R ó C. Se dice que la matriz
A tiene m filas y n columnas.

En el espacio de matrices Mm×n (K) está definida la operación suma:

A + B = C, donde cij = aij + bij
además podemos multiplicar por un escalar del cuerpo
λA = B donde bij = λaij
El producto de matrices AB está definido cuando el número de columnas de A es igual
al número de filas de B.

    
a11 · · · a1n b11 · · · b1r c11 · · · c1r
 .. . .. ..   ... . . . ...
.   .. . . . 
= . . .. 
 
 .
am1 · · · amn bn1 · · · bnr cm1 · · · cmr
| {z }| {z } | {z }
∈Mm×n ∈Mn×r ∈Mm×r

donde n
X
cij = aik bkj
k=1

Se llama matriz traspuesta AT de la matriz A a la que resulta de cambiar filas por

columnas, esto es
aTij = aji

Un tipo de matrices particular es el de matrices cuadradas que son aquellas que tienen el
mismo número de filas que de columnas. Se dice que una matriz cuadrada A es invertible
cuando existe una matriz A−1 tal que AA−1 = A−1 A = I donde I es la matriz identidad,
todas sus entradas son cero salvo las de la diagonal principal, que son 1.
 
1 0 ··· 0
 0 1 ··· 0 
In =  .. .. . . .. 
 
 . . . . 
0 0 ··· 1
22 CAPÍTULO 2. ÁLGEBRA

2.3. Sistemas en forma matricial

En este apartado estudiaremos la resolución de sistemas de ecuaciones lineales utilizando
el cálculo matricial.

Observemos en primer lugar que el sistema



 a11 x1 + a12 x2 + · · · + a1n xn = b1
a21 x1 + a22 x2 + · · · + a2n xn = b1

(2.1)

 ··· ··· ···
am1 x1 + am2 x2 + · · · + amn xn = bm


puede escribirse del siguiente modo

A·X =B

donde

     
a11 a12 ··· a1n x1 b1
 a21 a22 ··· a2n   x2
 , B =  b2  ,
  
A= , X = 
 ··· ··· ···   ···   ··· 
am1 am2 ··· amn xn bm

Operaciones elementales entre filas

Las operaciones elementales empleadas en un sistema para llevarlo a otro equivalente,
normalmente en forma más sencilla se traducen ahora en las operaciones elementales
entre filas de la matriz de coeficientes A. Ası́ tenemos

1. Intercambiar dos filas cualesquiera.

2. Multiplicar por una constante no nula una fila de la matriz.

3. Reemplazar cualquier fila por la suma de ésta y otra ecuación multiplicada por una
constante.

Equivalencia de matrices
Se dice que dos matrices A y B son equivalentes por filas si una se obtiene de la otra
por medio de las operaciones elementales entre filas descritas en el apartado anterior y se
escribe A ∼ B.
2.3. SISTEMAS EN FORMA MATRICIAL 23

Definición 2.1. Dada una matriz A, diremos que una fila de A es no nula si tiene por lo
menos un elemento distinto de cero. Además llamamos entrada principal de una fila no
nula al primer elemento por la izquierda que es no nulo. Decimos que la matriz A está en
forma escalonada si:

i) Todas las filas no nulas están por encima de las nulas.

ii) Debajo de cada entrada principal sólo hay ceros.

Ejemplo 2.1. Las siguientes matrices están en forma escalonada

   
2 4 0 0 0 0 0 2 2 6 5 8 0
 0 2 3 2 1 1 9   0 −1 1 1 3 3 
   
A=
 0 0 0 0 3 0 −3  , B =  0
  0 5 9 0 3 

 0 0 0 0 0 1 −2   0 0 0 −2 3 8 
0 0 0 0 0 0 0 0 0 0 0 1 3

Definición 2.2. Si una matriz A que está en forma escalonada verifica que sus entradas
principles son todas iguales a 1 y por encima de éstos sólamente hay ceros, entonces se
dirá que A está en escalonada reducida o forma canónica por filas.

Ejemplo 2.2. Las siguientes matrices están en forma canónica por filas:

   
1 0 3 1 0 0 4 1 0 0 0 0 2

 0 1 3 2 0 0 9 

 0 1
 0 0 0 3 

A=
 0 0 0 0 1 0 −3  , B = 

 0 0 1 0 0 4 

 0 0 0 0 0 1 −2   0 0 0 1 0 5 
0 0 0 0 0 0 0 0 0 0 0 1 6

Método de Gauss

Teorema 2.2. Toda matriz A tiene una única forma canónica.

Para trasformar una matriz en otra equivalente por filas y que tenga forma canónica se
sigue el siguiente algoritmo.
24 CAPÍTULO 2. ÁLGEBRA

1. Comezamos por la primera columna por la izquierda y elegimos el elemento distinto

de cero que esté más arriba.

2. Si el elemento anterior no está en la primera fila, lo llevamos a ella por una permu-
tación de filas.

3. Creamos ceros debajo de esta primera entrada principal.

4. Si al realizar los pasos anteriores, aparecen filas nulas, se situarán por debajo de las
filas no nulas.

5. Pasamos a la siguiente columna que contenga una entrada principal y se aplican los
pasos 1 a 4 a la submatriz que resulta al suprimir la primera fila y las columnas a
la izquierda de la que estamos tratando.

6. Se repite el punto anterior con todas las columnas que contengan entradas princi-
pales.

7. Empezando por la primera columna a la derecha con entrada principal, se crean por
medio de operaciones elementales ceros encima de ella. Si esta entrada no es igual
a 1 dividimos por su valor todos los elementos de la su fila para que ésta entrada se
conviera en 1.

8. Realizamos el paso 7 con las columnas que contienen entradas principales, siguiendo
el orden de derecha a izquierda.

Observación Los pasos 1 a 6 del proceso anterior nos proporcionan una forma escalonada
de A. Este proceso se conoce como eliminación gaussiana que en términos de sistemas de
ecuaciones lineales se conoce como método de eliminación de Gauss. Además los pasos 7
a 8 nos serán útiles para calcular la inversa A−1 de A, como veremos más adelante.
 
Si (A|B) es la matriz ampliada correspondiente al sistema (1) y A|eB
e su forma esca-
lonada,
  podemos establecer el siguiente criterio: El sistema (1) es compatible si y sólo si
B
e no contiene entradas principales.

 
En caso de que el sistema (1) sea compatible, la forma escalonada A|
eB e nos da más
información.

 
Definición 2.3. Sea A|B es la forma escalonada de la matriz ampliada correspondiente
e e
al sistema (1). Las variables correspondientes a columnas con alguna entrada principal de
A
e se llaman variables principales y las otras se llaman variables libres.

Tenemos ahora el siguiente

2.3. SISTEMAS EN FORMA MATRICIAL 25

 
Criterio 2.1. Si A|B es la forma escalonada de la matriz del sistema (1) y el sistema
e e
es compatible, se tiene:

i) El sistema es compatible determinado si y sólo si no hay ninguna variable libre.

ii) El sistema es compatible indeterminado si existe, al menos, una variable libre.

Ejemplo Consideremos el sistema homogéneo


 2x1 + 5x2 − 33x3 = 0
−2x1 + 3x2 − 23x3 = 0
15x1 + 2x2 + x3 = 0


cuya matriz ampliada es

 
2 5 33 0
 −2 3 23 0 
15 2 1 0

Hallamos la forma escalonada

   
2 5 −33 0 2 5 −33 0
 −2 3 −23 0  →  0 8 −56 0 →
15 2 1 0 15 2 1 0
   
2 5 −33 0 2 5 −33 0
 0 8 −56 0  →  0 8 −56 0 
0 −71
2
497
2
0 0 0 0 0

Se tiene entonces que es sistema es compatible porque la última fila de la matriz ampliada
no contiene entradas principales y además indeterminado porque una de las columnas
del bloque Ae no contiene entradas principales, concretamente la columna 3; es decir, la
variable x3 correspondiente a esta columna es libre. El sistema es equivalente al siguiente


2x1 + 5x2 − 33x3 = 0
8x2 − 56x3 = 0

de donde se tiene la solución por sustitución hacia arriba

x1 = −x3 , x2 = 7x3 , x3 ∈ R
26 CAPÍTULO 2. ÁLGEBRA

Ejemplo Sea el sistema


 2x1 + 5x2 − 3x3 = 0
6x1 − 3x2 + 2x3 = 0
8x1 − 2x2 + x3 = 0


Hallamos la forma escalonada de la matriz ampliada

   
2 5 −3 0 2 5 −3 0
 6 −3 2 0  →  0 −18 11 0 →
8 −2 1 0 8 −2 1 0
 
2 5 −3 0
 0 −18 11 0 
0 0 − 49 0
y el sistema es equivalente a


 2x1 + 5x2 − 3x3 = 0
−18x2 + 11x3 = 0
− 49 x3 = 0


que es compatible determinado y cuya solución es

x1 = 0, x2 = 0, x3 = 0

2.4. Determinantes. Regla de Cramer

Sea A una matriz cuadrada n × n; es decir, de n filas y n columnas.

 
a11 a12 · · · a1n
 a21 a22 · · · a2n 
A=
 
.. .. ... .. 
 . . . 
an1 an2 · · · ann

Si llamamos Ai,j a la matriz que resulta de suprimir en A la fila i y la columna j, tenemos

que

det A = (−1)1+1 a1,1 det A1,1 + (−1)1+2 a1,2 det A1,2 + · · · + (−1)1+n a1,n det A1,n
que es lo que se llama desarrollo del determinante de A, también denotado |A|, por la
primera fila.
2.4. DETERMINANTES. REGLA DE CRAMER 27

Menores y cofactores
Consideremos una matriz cuadrada A = (aij ) de orden n, (n × n) y sea Aij la submatriz
cuadrada de orden n − 1 que resulta al suprimir la fila i y la columna j. Entonces al
determinante |Aij | le llamamos menor del elemento aij de A, y llamamos cofactor o
adjunto de aij al número
Cij = (−1)i+j |Aij |

El siguiente teorema es de gran utilidad para el cálculo de determinantes

Teorema 2.3. El determinante de la matriz A = (aij ) es igual a la suma de los productos

de los elementos de una fila o columna por sus respectivos cofactores. Si desarrollamos
por la fila i resultarı́a
n
X
|A| = ai1 Ci1 + ai2 Ci2 + · · · + ain Cin = aij Cij
j=1

y si lo hacemos por la columna j serı́a

n
X
|A| = a1j C1j + a2j C2j + · · · + anj Cnj = aij Cij
i=1

Propiedades de los determinantes

Las siguientes propiedades pueden deducirse del teorema anterior

1. El determinante de una matriz A y su traspuesta AT son iguales.

2. El determinante del producto es el producto de los determinantes, |AB| = |A||B|.

3. Si A tiene una fila (o columna) de ceros, entonces |A| = 0.

4. Si A tiene dos filas iguales (o columnas), entonces |A| = 0.

5. Si A es triangular, entonces |A| = a11 a22 · · · ann .

6. Si B resulta de A por el intercambio de dos filas (o columnas), entonces |B| = −|A|.

7. Si B resulta al multiplicar una fila (o columna) de A por un escalar λ, entonces

|B| = λ|A|.

8. Si B resulta al sumar el múltiplo de una fila (o columna) a otra de A, entonces

|B| = |A|.
28 CAPÍTULO 2. ÁLGEBRA

9. Si A, B y C son matrices con todas sus filas iguales, salvo la i-ésima que es respecti-
vamente (αbi1 +βci1 , αbi2 +βci2 , · · · , αbin +βcin ), (bi1 , bi2 , · · · , bin ) y (ci1 , ci2 , · · · , cin ),
entonces det A = α det B + β det C

Ejemplo 2.3. Calcular el determinante de la matriz cuadrada n × n

 
a+1 a a ··· a

 a a+1 a ··· a 

A=
 a a a+1 ··· a 

 .. .. .. ... .. 
 . . . . 
a a a ··· a + 1

Si restamos la última fila a todas las anteriores tenemos

1 0 0 ··· −1
0 1 0 ··· −1
det A = 0 0 1 ··· −1
.. .. .. .. ..
. . . . .
a a a ··· a + 1
y ahora restamos a la última fila cada una de las anteriores multiplicada por a, y tenemos

1 0 0 ··· −1
0 1 0 ··· −1
det A = 0 0 1 ··· −1 = na + 1
.. .. .. .. ..
. . . . .
0 0 0 · · · na + 1

Ejemplo 2.4. Calcular el determinante

1 x1 x21
1 x2 x22
1 x3 x23

Solución: Haciendo las transformaciones de columnas

Cj → Cj − x1 Cj−1
2.4. DETERMINANTES. REGLA DE CRAMER 29

tenemos
1 x1 x21 1 0 0
x2 − x1 x2 (x2 − x1 )
1 x2 x22 = 1 x2 − x1 x2 (x2 − x1 ) = =
x3 − x1 x3 (x3 − x1 )
1 x3 x23 1 x3 − x1 x3 (x3 − x1 )
1 x2
= (x3 − x1 )(x2 − x1 ) = (x3 − x1 )(x2 − x1 )(x3 − x2 )
1 x3
Ejemplo 2.5. Probar que
1 x1 · · · xn−1
1
1 x2 · · · xn−1
2
Y
Vn = .. .. .. .. = (xj − xi )
. . . . i<j
n−1
1 x3 · · · x3

Solución: Iigual que en el ejercicio anterior hacemos las transformaciones de columnas

Cj → Cj − x1 Cj−1
y tenemos entonces
1 0 0 0 ··· 0
2 n−2
1 x2 − x1 x2 (x2 − x1 ) x2 (x2 − x1 ) · · · x2 (x2 − x1 )
Vn = .. .. .. ... .. =
. . . .
1 xn − x1 xn (xn − x1 ) xn (xn − x1 ) · · · xn−2
2
n (xn − x1 )
= (xn − x1 ) · · · (x2 − x1 )Vn−1
y tenemos el resultado por recurrencia

Rango de una matriz

Un concepto muy relacionado con el de determinante de una matriz es el de rango, aunque,
ası́ como el de determinante era aplicable únicamente a las matrices cuadradas, el rango
se estudia en cualquier tipo de matriz.

Definición 2.4. El rango de una matriz A es mayor entero n tal que A contiene una
submatriz n × n con determinante distinto de cero.

Ejemplo 2.6. Calculemos el rango de la matriz

 
1 −2 3 5
 3 −8 2 7 
−2 6 1 −2
30 CAPÍTULO 2. ÁLGEBRA

El determinante de la submatriz de orden 2 que se obtiene al eliminar la fila 3 y las

columnas 3 y 4 es

 
1 −2
det = 14 6= 0
3 −8
y por lo tanto el rango de A es mayor o igual a 2. Por otra parte las submatrices de A de
orden 3 son

       
1 −2 3 1 3 5 1 3 5 −2 3 5
 3 −8 2  ,  3 2 7  ,  3 2 7  y  −8 2 7 
−2 6 1 −2 1 −2 −2 1 2 6 1 −2
pero sus determinantes vales todos 0, por lo tanto, se tiene que el rango de A es 2.

Regla de Cramer

Definición 2.5. Se dice que un sistema es de Cramer si la matriz A del sistema es una
matriz cuadrada y |A| =
6 0.

El siguiente teorema se conoce como la Regla de Cramer y nos proporciona un método

para la resolución de sistemas de ecuaciones que son de Cramer. Como veremos más
adelante este método se podrá extender a sistemas más generales.

Teorema 2.4 (Regla de Cramer). Sea A una matriz cuadrada de orden n y no singular.
Para cualquier B = (b1 , · · · , bn )T , el sistema de ecuaciones lineales

a11 · · · a1j · · · a1n

    
x1 b1
 .. . . . .. .. .
. ..   .   .. 
 . .   ..   . 
 ai1 · · · aij · · · ain   xi  =  bi 
    
 . .
. . ... . . . ..   .   .. 
 .. .   ..   . 
an1 · · · an2 · · · ann xn bn
| {z } | {z } | {z }
A X B

es compatible y tiene como única solución

det A1 (B) det A2 (B) det An (B)

x1 = , x2 = , · · · , xn =
det A det A det A
donde Ai (B) es la matriz que se obtiene al sustituir en A la columna j-ésima por B,
2.5. TEOREMA DE ROUCHÉ-FROBENIUS 31

a11 · · · b1 · · · a1n
 
 .. . . . .. . . . ..
 . . .


Ai (B) =  ai1 · · · bi · · · ain
 

 . .
 .. . . ... .. .
. ..


an1 · · · bn · · · ann

Matriz inversa
Supongamos que A es una matriz n × n y det A 6= 0. Nos planteamos el problema de
calcular la matriz inversa de A, tenemos entonces

AA−1 = I (2.2)

Si llamamos Xi a la columna i-ésima de A−1 , podemos expresar (2) como n sistemas

lineales

AXi = Ei i = 1, · · · , n
donde Ei es la columna i-ésima de la matriz identidad. Resolviendo el sistema i por la
regla de Cramer, es decir, calculando la columna i-ésima de A−1 tenemos

det A1 (Ei ) det A2 (Ei ) det An (Ei )

x1i = , x2i = , · · · , xni =
det A det A det A
y de esta forma podemos escribir

T
C11 C12 · · · C1n

.
1  C21 C22 · · · ..

A−1 =
 
det A 
 .. .. .. .. 
. . . . 
Cn1 Cn2 · · · Cnn
donde Cij es el cofactor correspondiente al elemento aij .

2.5. Teorema de Rouché-Frobenius

Teorema 2.5 (Rouché-Frobenius). Un sistema de ecuaciones lineales AX = B es com-

patible si y sólo si rg A = rg( A|B). Explı́citamente, el sistema de ecuaciones


 a11 x1 + a12 x2 + · · · + a1n xn = b1
a21 x1 + a22 x2 + · · · + a2n xn = b2

S:

 ··· ··· ··· ··· ···
am1 x1 + am2 x2 + · · · + amn xn = bm

32 CAPÍTULO 2. ÁLGEBRA

tiene solución si las matrices

   
a11 a12 · · · a1n a11 a12 · · · a1n b1
 a21 a22 · · · a2n   a21 a22 · · · a2n b2 
A=  ···
 y A|B =  
··· ··· ···   ··· ··· ··· ··· ··· 
am1 am2 · · · amn am1 am2 · · · amn bm

tienen el mismo rango.

Demostración. El sistema S podemos escribirlo considerando los vectores columna de A

de la siguiente forma
       
a11 a12 a1n b1
 a21   a22   a2n   b1 
 ..  x1 +  ..  x2 + · · · +  ..  xn =  ..  (2.3)
       
 .   .   .   . 
am1 am2 amn b1

Supongamos primero que el sistema tiene solución, entonces existen ciertos números
x1 , x2 , . . . , xn de forma que el vector columna B es combinación lineal de los vectores
columna de A. Resulta entonces que si añadimos la columna B a la matriz A, ésta no
aumenta su rango, por lo que rg A = rg( A|B).

Reciprocamente, si rg A = rg( A|B), esto se puede interpretar como que la columna B es

combinación lineal de las columnas de A. Entonces existen ciertos escalares, que llamamos
x1 , x2 . . . . , xn que hacen que el vector B sea combinación lineal de los vectores columna de
la matriz A. Es decir tendrı́amos una expresión como (2.3) y esto significa que el sistema
tiene a los escalares xi como solución.
Definición 2.6. Se dice que un sistema de ecuaciones lineales es:

Compatible determinado si tiene una única solución.

Compatible indeterminado si tiene infinitas soluciones.
Incompatible si tiene no tiene solución

El teorema de Rouché-Frobenius se enuncia en ocasiones de manera más completa, ha-

ciendo alusión a los casos en los que el sistema es indeterminado. Aquı́ enunciamos esta
segunda parte como corolario que permite distinguir los casos en los que hay una o infinitas
soluciones.
Corolario 2.1. Sea r el rango de la matriz A del sistema del teorema anterior. Si el
sistema tiene solución (rg A = rg( A|B)), puede ocurrir:

r = n y en este caso el sistema es compatible determinado.

r < n y el sistema es compatible indeterminado.
2.6. ESPACIOS VECTORIALES 33

Demostración. Si el rango de la matriz A (y de su ampliada) es r entonces existe un

menor de orden r en la matriz A. Podemos suponer que se trata del menor formado por
las r primeras filas y r primeras columnas y escribimos el sistema S de la siguiente forma:

 a11 x1 + · · · + a1r xr = b1 − a1r+1 xr+1 − · · · − a1n xn
S: ··· ··· ··· ··· ··· ··· ···
ar1 x1 + · · · + arr xr = br − arr+1 xr+1 − · · · − arn xn


Si consideramos la matriz  
a11 · · · a1r
Ar =  ... . . . ... 
 
ar1 · · · arr
tenemos que |Ar | = 6 0 y podemos aplicar la Regla de Cramer para expresar las incógnitas
x1 , . . . , xr en función de las xr+1 , . . . , xn . Esto quiere decir que el conjunto de soluciones
del sistema depende de n − r parámetros. En resumen pueden darse los siguientes casos:

r = n y |A| =
6 0 y el sistema tiene solución única.

r < n y el sistema tiene infinitas soluciones. El conjunto de soluciones depende de

n − r parámetros.

2.6. Espacios vectoriales

Un espacio vectorial es una estructura algebraica compuesta por dos conjuntos: uno V
cuyos elementos se llaman vectores y otro K de escalares. El conjunto V está dotado de
una ley de composición interna (que denotaremos con el signo ”+”) para la que V es un
grupo conmutativo. Por otra parte K tiene otra operación interna (denotada con el signo
”·”) con la que este conjunto es un cuerpo y además esta operación es ley de composición
externa que asocia a cada λ ∈ K y v ∈ V un elemento λv ∈ V . Estas dos operaciones
están relacionadas por la propiedad distributiba. Concretamente podemos establecer

Definición 2.7. Un espacio vectorial es una estructura algebraica formada por el par
(V, K) donde V es el conjunto de vectores y K el de los escalares. Los elementos de V
se llaman vectores y nosotros los notaremos en negrita, y los elementos de K se llaman
escalares. Existen dos operaciones, la primera es la suma de vectores
+
V × V −→ V

y, la segunda que es una ley de composición externa

·
K × V −→ V

y se verifican las siguientes propiedades

34 CAPÍTULO 2. ÁLGEBRA

1. Propiedad asociativa para la suma: u + (v + w) = (u + v) + w

2. Propiedad conmutativa para la suma: u + v = v + u
3. Esistencia de elemento neutro para la suma: 0 + u = u
4. Existencia de elemento simétrico para la suma: Para cada u ∈ V existe u0 ∈ V tal
que u + u0 = 0
5. Propiedad asociativa para el producto por escalares: r(su) = (rs)u
6. Propiedad distributiva con respecto de la suma de vectores: r(u + v) = ru + rv
7. Propiedad distributiva con respecto de la suma de escalares: (r + s)u = ru + su
8. Existencia de elemento neutro para el producto por escalares: 1u = u

Ejemplos Los siguientes conjuntos son algunos ejemplos clásicos de espacios vectoriales.

1. El conjunto de las n-uplas (x1 , · · · , xn ) de números reales, o elementos del producto

cartesiano Rn , forman un espacio vectorial sobre el cuerpo R.
2. El conjunto de las matrices n × m con coeficientes reales respecto de las operaciones
de suma de matrices y producto de un escalar por una matriz habituales.
3. El espacio P6n de polinomios con coeficientes en R de grado menor o igual que
n, donde la suma y el producto por escalares son los siguientes: si p(x) = an xn +
an−1 xn−1 + · · · + a1 x + a0 y q(x) = bn xn + bn−1 xn−1 + · · · + b1 x + b0 son dos polinomios
de grados n y m respectivamente y λ ∈ R, entonces:

p(x) + q(x) = (an + bn )xn + (an−1 + bn−1 )xn−1 + · · · + (a1 + b1 )x + a0 + b0

y
λp(x) = λan xn + λan−1 xn−1 + · · · + λa1 x + λa0

4. El espacio P de polinomios con coeficientes en R.

5. El espacio de funciones f : R −→ R.

Proposición 2.1. En un espacio vectorial se verifican las siguientes propiedades:

1. 0u = 0
2. r0 = 0
3. (−r)u = −(ru) = r(−u)
2.6. ESPACIOS VECTORIALES 35

4. Si ru = 0, se tiene r = 0 o u = 0
5. Si ru = su, y u 6= 0, entonces r = s.
6. r(u − v) = ru − rv

Definición 2.8. Dado un espacio vectorial V y un subconjunto W ⊂ V , se dice que W

es un subespacio vectorial cuando también es un espacio vectorial.

El siguiente resultado proporciona una caracterización de los subespacios vectoriales.

Teorema 2.6. Dado un espacio vectorial V , el subconjunto W ∈ V es un subespacio
vectorial de V si ∀w1 , w2 ∈ W y λ, µ ∈ R se tiene que

λw1 + µw2 ∈ W

Demostración. Si el conjunto W verifica la propiedad del enunciado, entonces las opera-

ciones definidas en V están definidas en W y su imagen está en W . Todas las propiedades
de la definición de espacio vectorial las hereda W de V . Lo único que tenemos que probar
es que el elemento neutro está en W y el simétrico de cada elemento de W también está
en W . Para ver que 0 ∈ W hacemos λ = µ = 0 y la condición del enunciado implica que
0 ∈ W . Si v ∈ W entonces para λ = −1 y µ = 0 tenemos

−v ∈ W

Definición 2.9. Sean U y V dos espacios vectoriales. Llamamos suma de U y V y lo

denotamos U + V al espacio vectorial que resulta al considerar todos los vectores u + v
donde u ∈ U y v ∈ V .
Observación 2.1. A partir de este punto, salvo que se indique otra cosa, trabajaremos
con el espacio vectorial Rn con R.
Definición 2.10. Sean v1 , · · · , vn ∈ V . Dados n números reales λ1 , · · · , λn , el vector

v = λ1 v1 +, · · · , +λn vn

se llama combinación lineal de los vectores v1 , · · · , vn .

Dependencia e independencia lineal

Definición 2.11. Dado un conjunto de vectores C = {v1 , . . . , vn ∈ W }, diremos que es

un conjunto ligado o linealmente dependiente si alguno de ellos es combinación lineal
de los otros. En caso contrario; es decir, que ninguno de ellos sea combinación lineal de
los otros diremos que el conjunto es linealmente independiente.
36 CAPÍTULO 2. ÁLGEBRA

Definición 2.12. Dada una familia de vectores A = {v1 , · · · , vn }, llamamos clausura

lineal de A al conjunto
 
X 
L[B] = λj vj : B ⊂ A finito, λi ∈ R
 
vj ∈B

es decir, todas las posibles combinaciones lineales de elementos de A. La clausura lineal

de A también suele notarse hAi.

Es fácil demostrar que L[A] es un subespacio vectorial de V . además, si W es un subespacio

vectorial que contiene a A, entonces L[A] ⊂ W . Ası́ podemos decir que L[A] es el menor
espacio vectorial que contiene a A.
Definición 2.13. Sea A una familia de vectores de V . Diremos que A es un sistema de
generadores de V si L[A] = V .

El siguiente resultado proporciona un criterio muy útil para el estudio de la dependencia

o independencia lineal de un conjunto de vectores.
Proposición 2.2. Una familia A = {v1 , · · · , vn } de vectores no nulos de V es linealmente
independiente si dada una combinación lineal de elementos de A que dé como resultado
el vector cero: X
λi vi = 0V ,
vi ∈A

se tiene que todos los λi son iguales a 0.

Observación 2.2. En la situación de la proposición anterior, si de la combinación lineal
no resultan todos los λi iguales a cero, se puede despejar un vector vi como combinación
lineal de los restantes, y el conjunto no serı́a linealmente independiente.

Ejemplo Sea v ∈ V un vector distinto de 0V , entonces v es linealmente independiente.

Ejemplo Sean v1 = (1, −2) y v2 = (2, 0) dos vectores de R2 . Veamos que son linealmente
independientes. Una combinación lineal de ambos igualada a cero, nos proporciona la
ecuación

λ1 (1, −2) + λ2 (2, 0) = (0, 0)

de donde se tiene el sistema de ecuaciones lineales


λ1 + 2λ2 = 0
−2λ1 = 0
2.6. ESPACIOS VECTORIALES 37

y tiene como solución única λ1 = 0 y λ2 = 0.

Ejemplo El conjunto P = {3x − 1, x + 2} es linealmente independiente dentro del espacio

vectorial de polinomios P en la variable x con coeficientes reales. En efecto, si hacemos
λ1 (3x − 1) + λ2 (x + 2) = 0 tenemos


3λ1 + λ2 = 0
−λ1 + 2λ2 = 0
y tiene como solución única λ1 = 0 y λ2 = 0.

Definición 2.14. Un sistema de generadores B de un espacio vectorial V es una base,

cuando sus elementos son linealmente independientes.

Lema 2.1. Sea V un espacio vectorial y S = {v1 , . . . , vm } una base de V y consideremos

los vectores w1 , . . . , wn . Si n > m, entonces los vectores w1 , . . . , wn son linealmente
dependientes.

Demostración. Supongamos lo contrario; esto es, que los vectores w1 , . . . , wn son lineal-
mente independientes. Bajo este supuesto, junto con la hipótesis de que S constituye una
base de V vamos a ir sustituyendo los vectores de esta base por vectores wi . Comenzamos
con el primero, como S es base podemos escribir
w1 = α1 v1 + α2 v2 + · · · + αn vn

En la combinación lineal anterior no puede ser α1 = α2 = · · · = αn = 0, porque entonces

serı́a w1 = 0 y hemos supuesto que los wi son linealmente independientes. Entonces, salvo
reordenación de los vi suponemos que α1 6= 0 y entonces resulta
v1 = α1−1 w1 − α1−1 α2 w2 − · · · − α1−1 αm wm
esto demuestra que el subespacio que generan w1 , v2 , . . . , vm contiene a los v1 , v2 , . . . , vm ,
y esto implica que contiene a V .

La demostración sigue, reemplazando paso a paso los vectores, ya que si tenemos que
w1 , . . . , wr , vr+1 , . . . , vm generan V entonces existirá una combinación lineal
wr+1 = β1 w1 + · · · + βr wr + γr+1 vr+1 + · · · + γm vm
y por la misma razón que antes no puede ser γr+1 = · · · = γm = 0, ası́ que si suponemos
que γr+1 6= 0 obtendrı́amos de la combinación lineal
γr+1 vr+1 = wr+1 − β1 w1 − · · · − βr wr − γr+2 vr+2 − · · · − γm vm
38 CAPÍTULO 2. ÁLGEBRA

que vr+1 es combinación lineal de w1 , . . . , wr+1 , vr+2 , . . . , vm , y entonces estos vectores

generan V .

Siguiendo este proceso llegamos finalmente a que el conjunto T = {w1 , . . . , wm } gene-

ra V , pero como n > m tendrı́amos que wn podrá escribirse como combinación lineal
del conjunto T lo que es una contradicción ya que habı́amos supuesto que el conjunto
w1 , . . . , wn era linealmente independiente.

Teorema 2.7. Sea V un espacio vectorial y B y B 0 dos bases de V . Entonces B y B 0

tienen el mismo número de elementos.

Demostración. Sean B = {v1 , . . . , vm } y B 0 = {w1 , . . . , wn }. Por el teorema anterior se

tienen las desigualdades n ≤ m y m ≤ n, y entonces ha de ser m = n.

Definición 2.15. El número de elementos de una base de el espacio vectorial V se llama

dimensión del espacio vectorial.

Espacio vectorial cociente

Sea V un espacio vectorial sobre un cuerpo K y E ⊂ V un subespacio. En este apartado
vamos a definir el concepto de espacio vectorial cociente y demostraremos resultados
relativos a las dimensiones de los espacios tratados. En esta situación la estructura de E
como subgrupo de V permite definir una relación de equivalencia ∼E , de manera que

v ∼E v0 si y sólo si v − v0 ∈ E

Ejercicio 2.2. Comprobar que con la relación ∼E definida antes el conjunto V /E se dota
de estructura de espacio vectorial sobre K.

Definición 2.16. Sea V un espacio vectorial y E ⊂ V un subespacio. El espacio V /E

formado por las clases v + E, donde v ∈ V , se llama espacio vectorial cociente, (de V
módulo E).

Teorema 2.8. Sea E ⊂ V un subespacio vectorial de dimensión finita. Entonces se tiene

que:

1. dimE ≤ dimV
2.7. PRODUCTO ESCALAR. ORTOGONALIDAD 39

2. dim(V /E) = dimV − dim(E)

Demostración. Sea e1 , . . . , er una base de E aplicando el lema 2.1 podemos ampliar esta
base de forma que
e1 , . . . , er , v1 , . . . , vs
es una base de V . Esto demuestra la primera afirmación.

Vamos a ver que [v1 ], . . . , [vs ] es una base de V /E. Para ver que son un sistema de
generadores basta tener en cuenta que e1 , . . . , er , v1 , . . . , vs genera V y [e1 ] = · · · = [er ] =
[0].

Para ver que son linealmente independientes igualamos una combinación lineal arbitraria
a cero, y tenemos " s #
s
X X
[0] = λi [vi ] = λi vi
i=1 i=1

pero esto significa que

s
X
λi vi ∈ E
i=1

y entonces podemos expresar esta combinación como combinación de la base de E

λ1 v1 + · · · + λs vs = µ1 e1 + · · · + µr er

esto es s s
X X
λi vi − µj e j = 0
i=1 i=j

pero como e1 , . . . , er , v1 , . . . , vs son linealmente independientes tenemos

λ1 = λ2 = · · · = λs = µ1 = · · · = µr = 0

y en particular todos los λi = 0 que es lo que querı́amos probar.

2.7. Producto escalar. Ortogonalidad

Definición 2.17. Sea E un espacio vectorial sobre R. Se llama producto escalar a una
aplicación
h·, ·i : E × E → R
que cumple las siguientes propiedades:

i) hu, vi = hv, ui, ∀u, v ∈ E

40 CAPÍTULO 2. ÁLGEBRA

ii) hλu + µv, wi = λhu, wi + µhv, wi ∀λ, µ ∈ R, ∀u, v, w ∈ E

Si además se cumple

iii) hv, vi > 0, ∀v ∈ E\ {0E }

entonces el producto escalar es definido positivo.

Ejemplo 2.7. Si x = (x1 , x2 ) y y = (y1 , y2 ) son vectores de R2 , la aplicación

hx, yi = 2x1 y1 + 4x2 y2

define un producto escalar. Sin embargo

hx, yi = x1 y1 + x1 y2 + x2 y1 + x2 y2

no, ya que hx, xi = (x1 + x2 )2 y se tiene que para x = (1, −1) es hx, xi = 0 y no se cumple
la propiedad i) de la definición.

Ejercicio 2.3. Comprobar que en el espacio P2 de polinomios en t con coeficientes reales,

de grado menor o igual que 2, la aplicación
3
X
hp, qi = p(i)q(1)
i=1

define un producto escalar.

Definición 2.18. Sea E un espacio vectorial sobre R. Se llama noma a una aplicación

k·k:E →R

que satisface las siguientes propiedades

i) kvk > 0 ∀v ∈ E\0E y kvk = 0 si y sólo si v = 0E

ii) kλvk = |λ|kvk, ∀λ ∈ R y ∀v ∈ E
iii) kv + wk 6 kvk + kwk ∀v, w ∈ E
2.7. PRODUCTO ESCALAR. ORTOGONALIDAD 41

Proposición 2.3. Dado un producto escalar h·, ·i definido en un espacio vectorial E sobre
R, la aplicación p
kvk = + hv, vi
define una norma en E.

Decimos que un vector u ∈ V es unitario si kuk = 1. Si u 6= 0, entonces u/kuk es

unitario.

Definición 2.19. Sea h·, ·i un producto escalar definido en un espacio vectorial E sobre
R. Dados dos vectores u y v, se define el ángulo entre ellos como el número θ ∈ [0, π) tal
que
hu, ui
cos θ =
kukkvk

Definición 2.20. Dos vectores u y v del espacio E son ortogonales respecto del pro-
ducto escalar h·, ·i si hu, vi = 0. Una familia de vectores {u1 , · · · , un } se dice ortogonal
si hui , uj i = 0 para todos los pares (i, j) con i 6= j. La familia {u1 , · · · , un } se llama
ortonormal si es ortogonal y además todos sus vectores tienen norma 1.

El producto escalar definido en Rn por

hx, yi = x1 y1 + x2 y2 + · · · + xn yn

se llama producto escalar usual o estándard.

Ejemplo La familia (1, 2, 0), (−2, 1, 3), (6, −3, 5) es ortogonal en R3 con el producto es-
calar usual. Para obtener una familia ortonormal, dividimos cada vector por su norma,
ası́ obtenemos
1 1 1
√ (1, 2, 0), √ (−2, 1, 3), √ (6, −3, 5)
5 14 70

Definición 2.21. Sea A un conjunto de un espacio vectorial E con un producto escalar.

Se llama conjunto ortogonal a A y se denota A⊥ , al conunto de los vectores de E que son
ortogonales a todos los vectores de A.
42 CAPÍTULO 2. ÁLGEBRA

Proposición 2.4. El conjunto A⊥ es subespacio vectorial de E. Si A es subespacio vecto-

rial de E y éste es de dimensión finita, se tiene E = A ⊕ A⊥ y A⊥ se llama complemento
ortogonal de A.

Ejemplo Calcular el ángulo formado por las matrices

   
1 −2 0 1
A= B=
3 0 1 −1

respecto del producto escalar hA, Bi = trazaAB T en el espacio vectorial M2×2 de matrices
reales 2 × 2.

Solución: Tenemos
 
T −2 3
hA, Bi = trazaAB = traza =1
0 3

además
p √
kAk = + hA, Ai = 14
p √
kBk = + hB, Bi = 3

y finalmente se tiene que el ángulo θ es tal que

1
cos θ = √ √
14 3
es decir θ = arc cos √141√3 .

Proyección ortogonal

Definición 2.22. Sea E un espacio vectorial dotado de una norma y F un subespacio de

E y sea v ∈ E\F . Llamamos proyección ortogonal de v sobre F y lo denotamos ΠF (v)
al vector v0 de F que verifica

kv − v0 k = inf {kv − uk : u ∈ F } (2.4)

Proposición 2.5. Si la dimensión de F es n, existe un único vector v0 ∈ F que verifica

(3). Además, dada una base ortonormal {e1 , · · · , en } de F se tiene

n
X
v0 = hv, ei i ei
i=1
2.7. PRODUCTO ESCALAR. ORTOGONALIDAD 43

Teorema 2.9. Para cualesqiera u y v ∈ V se tiene

|hu, vi| ≤ kukkvk

Teorema 2.10. Sean u, v ∈ V , entonces:

1. (Teorema de Pitágoras) Si u, v son perpendiculares entre sı́, entonces

ku + vk2 = kuk2 + kvk2 .

2. (Ley del palalelogramo) Para cualesquiera u, v se tiene

ku + vk2 + ku − vk2 = 2kvk2 + 2kvk2 .

Demostración. Para la primera parte tenemos

ku + vk2 = hu + vi = hu, ui + 2hu, vi + hv, vi

= kuk2 + kvk2 .

La segunda se deja como ejercicio.

Teorema 2.11. Para cualesquiera v, w ∈ V tenemos

|hv, wi| ≤ kvkkwk

Teorema 2.12. Sean v, w ∈ V , entonces

kv + wk ≤ kvk + kwk
44 CAPÍTULO 2. ÁLGEBRA

Bases ortogonales
Se dice que una base {v1 , . . . , vn } de V es ortogonal si sus elementos son mutuamente
perpendiculares. Si además, cada uno de los vectores de la base es unitario, la base es
ortonormal.

Teorema 2.13. Sea V un espacio vectorial de dimensión finita con un producto escalar
definido positivo. Sea W ⊂ V un subespacio vectorial. Supongamos que W 6= V y sea
{w1 , . . . , wm } una base ortonormal de W . Entonces existen elementos wm+1 , . . . , wn de
forma que {w1 , . . . , wn } es una base ortonormal de V .

Demostración. (Ortogonalización de Gram-Schmidt). Podemos ampliar la base {w1 , . . . , wm }

a una base
{w1 , . . . , wm , vm+1 , . . . , vn }
que no es necesariamente ortogonal. Si llamamos Wm+1 al espacio generado por w1 , . . . , wm , vm+1
y con esta construimos una base ortogonal. Para ello le restamos a vm+1 su proyección
sobre el subespacio generado por {w1 , . . . , wm }. Entonces si llamamos

hvm+1 , w1 i hvm+1 , wm i
c1 = , . . . , c1 = .
hw1 , w1 i hwm , wm i

Sea
wm+1 = vm+1 − c1 w1 − · · · − cm wm
Resulta que wm+1 es perpendicular a w1 , . . . , wm . Por otra parte wm+1 6= 0 y vm+1
pertenece al subespacio generado por w1 , . . . , wm+1 , ya que

vm+1 = vm+1 + c1 w1 + · · · + cm wm

Tenemos entonces que {w1 , . . . , wm+1 } es una base ortogonal de Wm+1 . Se obtiene por
induccoón una base ortogonal para el espacio Wm+s , s = 1, . . . , n − m.

Corolario 2.2. Sea V un espacio vectorial de dimensión finita con un producto escalar
definido positivo. Entonces V tiene una base ortogonal.
2.8. APLICACIONES LINEALES 45

2.8. Aplicaciones lineales

Definición Dados dos espacios vectoriales E y F sobre R. La aplicación f : E → F es
lineal si
f (λu + µv) = λf (u) + µf (v), ∀λ, µ ∈ R, ∀u, v ∈ E

Definición Se llama imagen de la aplicación lineal f : U → V al conjunto de vectores

de V que son imagen por f de algún vector de U . La imagen de la aplicación f se suele
denotar Imf .

Por otra parte, se llama núcleo de la aplicación f al conjunto de vectores de U cuya

imagen es el vector cero, 0V de V . El núcleo de f suele notarse Kerf .

Sea f : U → V , una aplicación lineal como al principio de esta sección y sean {u1 , u2 , · · · , un }
y {v1 , v2 , · · · , vm } bases de U y V respectivamente. Tenemos entonces

f (u1 ) = a1,1 v1 + a2,1 v1 · · · + am,1 vm

f (u2 ) = a1,2 v1 + a2,2 v1 · · · + am,2 vm
.. .. .. ..
. . . .
f (un ) = a1,n v1 + a2,n v1 · · · + am,n vm

Entonces si w = ni=1 xi vi es la expresión de un vector cualquiera de U comoo combina-

P
ción lineal de los elementos de la base anterior, tenemos

n n m
! m n
!
X X X X X
f (w) = xi f (ui ) = xi aj,i vj = aj,i xi vj
i=1 i=1 j=1 j=1 i=1

De esta forma, si escribimos f (w) = y1 v1 + · · · + vm , tenemos la representación matricial

de la imagen del vector v por la aplicación f :

    
y1 a1,1 a1,2 · · · a1,n x1
 y2   a2,1 a2,2 · · · a2,n  x2 
=
    
 .. .. .. .. ..  .. 
 .   . . . .  . 
ym am,1 am,2 · · · am,n xn
| {z } | {z } | {z }
Y A X

Definición Llamamos a la matriz A anterior, matriz de la aplicación lineal f respecto de

las bases {u1 , u2 , · · · , un } y {v1 , v2 , · · · , vm }.
46 CAPÍTULO 2. ÁLGEBRA

Definición Llamamos rango de una aplicación lineal f al rango de la matriz asocida a f

con respecto de las bases canónicas.

Ejercicio 2.4. Demostrar que el rango de una aplicación lineal f no depende de las bases
que se elijan para hallar la matriz A.

Teorema 2.14. [Isomorfı́a] Sea f : E → F una aplicación lineal. Entonces la aplicación

f : E/Ker f → Im f

definida como f [e] = f (e) es un isomorfismo. Se tiene entonces

dim E = dim Imf + dim Kerf

Corolario 2.3. Si A es la matriz m × n de una aplicación lineal f : U → V entonces

dim Kerf = n − rg A

Corolario 2.4. Sean V y W dos subespacios vectoriales de un espacio vectorial E de

dimensión finita. Entonces

dim (V + W ) = dim V + dim W − dim (V ∩ W )

Demostración. Consideremos la aplicación

f : V −→ (V + W )/W

definida de manera natural f (v) = [v]. Fácilmente se ve que f es suprayectiva y que

Ker f = V ∩ W . Tenemos entonces por el teorema de isomorfı́a

dim V = dim (V ∩ W ) + dim (V + W )/W

= dim (V ∩ W ) + dim (V + W ) − dim W
2.8. APLICACIONES LINEALES 47

Ejemplo Sea f : R → R la aplicación lineal que envı́a los vectores (2, 1, 0), (1, 0, 1) y
(−2, 3, 1) en (1, 2, 1), (3, 0, 2) y (1, −1, 0) respectivamente. Hallar la matriz de f respecto
de la base canónica.

Solución: Expresando las condiciones del enunciado en forma matricial, si llamamos A

a la matriz pedida,

   
2 1 −2 1 3 1
A 1 0 3  =  2 0 −1 
0 1 1 1 2 0
y entonces

  −1
1 3 1 2 1 −2
A =  2 0 −1   1 0 3 
1 2 0 0 1 1
  
1 3 1 1/3 1/3 −1/3
=  2 0 −1   1/9 −2/9 8/9 
1 2 0 −1/9 2/9 1/9
 
5/9 −1/9 22/9
=  7/9 4/9 −7/9 
5/9 −1/9 13/9

Ejemplo Sea f : R → R la aplicación lineal cuya expresión en coordenadas es f (x1 , x2 , x3 ) =

(2x1 − 3x2 , x2 − x3 , x1 − 2x3 ). Hallar la matriz de la aplicación f repecto de la base
B = {(1, 0, 2), (−1, 0, 0), (2, 3, 0)}

Solución: De la expresión en coordenadas de f se tiene que la matriz de f respecto de

la base canónica en R3 es

 
2 −3 0
M = 0 1 −1 
1 0 −2

Observemos ahora que la matriz

 
1 −1 2
A= 0 0 3 
2 0 0
48 CAPÍTULO 2. ÁLGEBRA

es la que efectúa el cambio de base de la canónica con B, es decir, si v tiene coordenadas

(v1 , v2 , v3 ) con respecto de la base B, entonces AvT son las coordenadas de este vector en
la base canónica.

La matriz pedida en el enunciado es

MB = A−1 M A
 −1   
1 −1 2 2 −3 0 1 −1 2
=  0 0 3   0 1 −1   0 0 3 
2 0 0 1 0 −2 2 0 0
 
−3/2 −1/2 1
=  −29/6 3/2 8 
−2/3 0 1

Ejemplo Sea P2 el espacio vectorial de los polinomios de grado menor o igual que 2.
Consideramos la aplicación lineal f : P2 → P2 , definida por f (1) = 1, f (t − 1) = t y
f (t2 − t) = t2 − 1.

a) Calcular Imf y ker f y sus dimensiones.

b) Hallar la matriz de f respecto de la base B = {1, t, t2 }

Solución: a) Es fácil comprobar que el conjunto {1, t − 1, t2 − t} es una base de P2 .

Primero veamos que es un sistema linealmente independiente, para esto supongamos que
una combinación lineal de ellos es igual al polinomio cero, 0P2 :

α1 1 + α2 (t − 1) + α3 (t2 − t) = 0P2
entonces reagrupando tenemos

(α1 − α2 ) + t(α2 − α3 ) + t2 α3 = 0P2

de donde se tiene necesariamente que α3 = α2 = α1 = 0, lo que prueba la independencia
lineal.

Además si p(t) = a1 + a2 t + a3 t2 es un elemento cualquiera de P2 y hacemos

a1 + a2 t + a3 t2 = α1 + α2 (t − 1) + α3 (t2 − t)
= (α1 − α2 ) + t(α2 − α3 ) + t2 α3
2.8. APLICACIONES LINEALES 49

e igualando coeficientes tenemos un sistema en las incógnitas α1 , α2 , α3 , que tiene una

única solución, ya que la matriz del sistema
 
1 −1 0
 0 1 −1 
0 0 1

tiene determinante distinto de cero y por lo tanto tiene rango 3.

Como cualquier polinomio p(t) puede expresarse en la base anterior tenemos que p(t) =
a1 + a2 (t − 1) + a3 (t2 − t), entonces para calcular ker f hacemos

f (p(t)) = f a1 + a2 (t − 1) + a3 (t2 − t)

= a1 f (1) + a2 f (t − 1) + a3 f t2 − t

= a1 1 + a2 t + a3 t2 − 1 = 0P2

y tiene como solución a1 = a2 = a3 = 0, por lo que resulta ker f = {0P2 }.

Por otra parte, de la igualdad

dim P2 = dim Imf + dim ker f

se tiene que dim Imf = 3, ya que dim P2 = 3.

b) De los datos del enunciado tenemos que

f (1) = 1 = (1, 0, 0)B

f (t − 1) = f (t) − f (1) = f (t) − 1 = t ⇒ f (t) = t + 1 = (1, 1, 0)B
f (t2 − t) = f (t2 ) − f (t) = f (t2 ) − (t + 1) = t2 − 1 ⇒ f (t2 ) = t2 + t = (0, 1, 1)B

y entonces la matriz de la aplicación f respecto a la base B es

 
1 1 0
 0 1 1 
0 0 1

Ejemplo Sea f : R3 −→ R4 la aplicación lineal definida por

f (x, y, z) = (x − y, y + z, x + y, x + y − 2z)

Hallar ecuaciones de Kerf y de Imf .

50 CAPÍTULO 2. ÁLGEBRA

Solución: Si {e1 , e2 , e3 } es la base canónica de R3 , los vectores

u1 = f (e1 ) = (1, 0, 1, 1)
u2 = f (e2 ) = (−1, 1, 1, 1)
u3 = f (e3 ) = (0, 1, 1, −2)

forman un sistema de generadores de Imf . Como estos vectores son linealmente indepen-
dientes, forman una base de Imf . Entonces podemos escribir

Imf = hu1 , u2 , u3 i = {(α − β, β + γ, α + β + γ, α + β − 2γ) : α, β, γ ∈ R} ⊂ R4

El Kerf está formado por los (x, y, z) ∈ R3 tales que

x − y = y + z = x + y = x + y − 2z = 0

y la solución de este sistema de ecuaciones es la trivial. Se tiene entonces

Kerf = {(0, 0, 0)}

y observamos que
dim R3 = 3; dim Imf = 3; dim Kerf = 0

Ejemplo. Sea R3 [x] el espacio vectorial de los polinomios de grado 3. Consideremos la

aplicación lineal F : R3 [x] −→ R3 [x] definida por

F (p(x)) = p0 (x)

Calcular las dimensiones de Imf y de Kerf .

Solución: Explı́citamene, la aplicación F es

F (a + bx + cx2 + dx3 ) = b + 2cx + 3dx3

El espacio vectorial R3 [x] tiene dimensión 4 (una base de este espacio está formada por
{1, x, x2 , x3 }) y la imagen de F está formada por los polinomios de grado 2, que tiene
dimensión 3. Por otra parte, Kerf está formada por los polinomios constantes (ya que su
derivada es el polinomio cero). Una base de Kerf estarı́a formada por un elemento, {1}.
En definitiva tenemos

dim R3 [x] = 4; dim Imf = 3; dim Kerf = 1

2.9. ISOMETRÍAS EN Rn 51

2.9. Isometrı́as en Rn
En esta sección estudiaremos las trasnformaciones lineales que preservan las distancias en
Rn . El concepto de distancia se tienen cuando está definida una norma, que a su vez es
consecuencia del producto escalar. Ver definiciones 2.17 y 2.18.

Definición 2.23. Una aplicación

f : Rn → Rn

es una isometrı́a si para cualquier par de vectores u, v de Rn se tiene

kf (u) − f (v)k = ku − vk (2.5)

√
donde kuk = + u1 , . . . , un .

Ejercicio 2.5. Probar que las siguientes aplicaciones son isometrı́as:

1. f : R2 → R2 , (x, y) 7→ (−x, y)

2. g : Rn → Rn , u 7→ −u

3. t : Rn → Rn , u 7→ a + u donde a es un vector fijo.

En lo que sigue nos vamos a fijar en las isometrı́as que fijan el vector cero. El siguiente
teorema caracteriza estas aplicaciones.

Teorema 2.15. Sea f : Rn → Rn una aplicación. Entonces son equivalentes:

1. f es una isometrı́a y f (0) = 0.

2. f preserva el producto escalar; esto es: f (u) · f (v) = u · v, ∀ u, v ∈ Rn .

Demostración. 1. (1 ⇒ 2) Suponemos que f cumple 1 en el enunciado del teorema.

De la propiedad

kf (u) − f (v)k = ku − vk (2.6)

tenemos, por una parte

kf (u)k = kuk (2.7)

52 CAPÍTULO 2. ÁLGEBRA

Porque f fija el vector 0.

Si ahora elevamos al cuadrado la ecuación (2.6) tenemos

(f (u) − f (v)) · (f (u) − f (v)) = (u − v) · (u − v)

f (u) · f (u) − 2f (u) · f (v) + f (v) · f (v) = u · u − 2u · v + v · v
kf (u)k2 − 2f (u) · f (v) + kf (v)k2 = kuk2 − 2u · v + kvk2

y teniendo en cuenta 2.7 se tiene

f (u) · f (v) = u · v (2.8)

2. (2 ⇒ 1) Supongamos ahora que f preserva el producto escalar; esto es, que se

verifica 2.8 para cualquier par de vectores. Entonces

kf (u) − f (v)k2 = (f (u) − f (v)) · (f (u) − f (v))

= f (u) · f (u) − 2f (u) · f (v) + f (v) · f (v)
= u · u − 2u · v + v · v
= (u − v) · (u − v)
= ku − vk2

Con lo que hemos demostrado que f es una isometrı́a. Queda por probar que f fija
el vector cero.
Si sustituimos en (2.8) los vectores u y v por el vector 0 tenemos que

0 = f (0) · f (0) = kf (0))k2

y tiene que ser f (0) = 0, (Definición 2.18, i)).

El siguiente teorema demuestra que las isometrı́as que fijan el origen son aplicaciones
lineales, además demuestra una propiedad importante que verifican las matrices de estas
aplicaciones. Estas matrices, y como consecuencia las aplicaciones que representan, forman
un grupo con la composición. Este es el Grupo Ortogonal.

Para probar el teorema necesitamos un lema previo.

Lema 2.2. Sea M una matriz n × n y u, v ∈ Rn , entonces

u · Mv = MT u · v

Demostración. Se deja como ejercicio.

2.9. ISOMETRÍAS EN Rn 53

Teorema 2.16. Sea f : Rn → Rn una aplicación. Entonces son equivalentes:

1. f es una isometrı́a y f (0) = 0.

2. f es lineal y la matriz A de f verifica AAT = I, (I representa la matriz identidad).

Demostración. 1. (1 ⇒ 2) Supongamos que se verifica 1, entonces por el teorema

anterior f preserva el producto escalar; es decir, para cualquier par de vectores se
tiene f (u) · f (v) = u · v.
Sea ahora un vector arbitrario w ∈ Rn , entonces

f (u + v) · f (w) = (u + v) · w

(f (u) + f (v)) · f (w) = f (u) · f (w) + f (v) · f (w)

= u · w + v · w = (u + v) · w

De los miembros de la derecha de las ecuaciones primera y tercera tenemos que para
cualquier terna de vectores u, v, w

f (u + v) · f (w) = (f (u) + f (v)) · f (w) (2.9)

Ahora por ser f una isometrı́a, si e1 , . . . , en es una base ortonormal de Rn entonces

f (e1 ), . . . , f (en ) también es una base ortonormal. De la igualdad 2.9 tenemos para
los vectores f (ei )

f (u + v) · f (ei ) = (f (u) + f (v)) · f (ei ), i = 1, . . . , n (2.10)

Cada miembro de esta igualdad expresa la proyección de cada uno de los vectores
f (u + v) y f (u) + f (v) sobre los vectores f (ei ); es decir, las coordenadas de esos
vectores con respecto a esa base.
Resula entonces de 2.10 que para cualesquiera u, v ∈ Rn

f (u + v) = f (u) + f (v)

De manera análoga se demuestra

f (λu) = λf (u)

que se deja como ejercicio.

2. (2 ⇒ 1) Sea A la matriz de f . Entonces, aplicando el lema 2.2 tenemos para cualquier
par de vectores de Rn

Au · Av = AT (Au) · v = (AT A)u · v

54 CAPÍTULO 2. ÁLGEBRA

en particular para los vectores de la base e1 , . . . , en se tiene

(AT A)ej · ei = δij

lo que quiere decir que

AT A = I.

Se deja como ejercicio la comprobación de este último paso de la demostración.

2.10. EJERCICIOS 55

2.10. Ejercicios
Ejercicio Comprobar si los siguientes subconjuntos son subespacios vectoriales de R3 .

i) V1 = {(x1 , x2 , x3 ) : x21 + x2 + x3 = 0}

ii) V2 = {(x1 , x2 , x3 ) : x1 + λx2 + x3 = 0, λ ∈ R}

iii) V3 = {(x1 , x2 , x3 ) : x1 + (λ + µ)x2 + x3 = 0, λ, µ ∈ R}

iv) V4 = {(x1 , x2 , x3 ) : x1 + x2 + x3 = 1}

Ejercicio Sean U y V dos subespacios vectoriales de un espacio vectorial W . Probar que

a) U ∩ V es subespacio vectorial de W , y por lo tanto de U y de V . Este resultado se

extiende claramente al caso de la intersección de n subespacios de W .

b) U + V es subespacio vectorial de W .

c) En general, U ∪ V no es subespacio vectorial de W . ¿En qué condiciones es U ∪ V

subespacio vectorial de W ?.

Ejercicio Sea F el espacio vectorial de las funciones de R en R, y sean

P = {f ∈ F : f (−x) = f (x), ∀x ∈ R}

I = {f ∈ F : f (−x) = −f (x), ∀x ∈ R}
los subconuntos formados por las funciones pares e impares. Probar que P e I son subes-
pacios vectoriales de F y que además se tiene P ⊕ I = F.

Ejercicio Sean U = L[B1 ] y V = L[B2 ] dos subespacios vectoriales de W donde B1 y B2

son bases.

i) Probar que B1 ∪ B2 es un sistema de genradores de U + V .

ii) Si U ∩ V 6= {0W } y B3 es una base de U ∩ V , completarla hasta obtener bases C1

y C2 de U y V respectivamente. Probar que

B3 ∪ (C1 \B3 ) ∪ (C2 \B3 )

es base de U + V .
56 CAPÍTULO 2. ÁLGEBRA

iii) Si U y V son espacios vectoriales de dimensión finita, demostrar, aplicando el apar-

tado anterior que

dim (U + V ) + dim (U ∩ V ) = dim U + dim V

Ejemplo Consideremos el subespacio W de R4 formado por los vectores (x1 , x2 , x3 , x4 )

tales que 
3x1 − x3 + 2x4 = 0
x1 − x2 + 2x4 = 0

1. Hallar las ecuaciones paramétricas de W .

2. Calcular una base de W .
3. Determinar si el vector (1, 2, −1, −3) pertenece a W .

Solución: El sistema de ecuaciones que define W es indeterminado y equivalente a


3x1 + 2x4 = x3
x1 + 2x4 = x2

cuya solución nos da la forma general de los vectores de W

(x1 , x1 + 2x4 , 3x1 + 2x4 , x4 )

donde puede verse como espacio generado por los vectores (1, 1, 3, 0) y (0, 2, 2, 1). Las
ecuaciones vectoriales de W serán

(x1 , x2 , x3 , x4 ) = λ(1, 1, 3, 0) + µ(0, 2, 2, 1)

y las paramétricas



 x1 = λ
x2 = λ + 2µ


 x3 = 3λ + 2µ
x4 = µ


Hemos visto que el espacio W tiene dimensión 2 y una base estarı́a formada por los
vectores e1 = (1, 1, 3, 0) y e2 = (0, 2, 2, 1). Para expresar el vector (1, 2, −1, −3) en esta
base sustituimos en las ecuaciones paramétricas y tenemos



 1 = λ
2 = λ + 2µ


 −1 = 3λ + 2µ
−3 = µ

2.10. EJERCICIOS 57

donde se ve claramente que resulta un sistema incompatible, esto quiere decir que (1, 2, −1, −3) ∈
/
W-

Ejemplo Consideremos en R5 los subespacios U y V dados por las siguientes ecuaciones


x1 − 2x4 = 0
U
x1 − x3 − x5 = 0



 x1 = λ
 x2 = λ + 2µ


V x3 = λ+µ
x4 = λ−µ




x5 = µ


i) Hallar sendas bases de U y de V .

ii) Expresar U ∩ V en ecuaciones paramétricas y determinar su dimensión.

iii) Calcular la dimensión de U + V .

Solución: El sistema que define U puede escribirse


x1 = 2x4
x1 = x3 + x5

ası́ podemos considerar x1 y x4 como variables dependientes y x2 , x3 y x5 como variables

independientes. de esta forma un vector genérico del subespacio U puede escribirse
 
1
x3 + x5 , x2 , x3 , (x3 + x5 ), x5
2

y el conjunto formado por los vectores u1 = (0, 1, 0, 0, 0), u2 = (1, 0, 1, 1/2, 0) y u3 =

(1, 0, 0, 1/2, 1) forman una base de U que tiene, de esta forma dimensión 3.

Por otra parte en las ecuaciones que definen V vemos que el conjunto formado por los
vectores v1 = (1, 1, 1, 1, 0) y v2 = (0, 2, 1, −1, 1) forman una base de este espacio, que por
lo tanto tiene dimensión 2.

Para encontrar las ecuaciones de U ∩V sustituimos las ecuaciones paramétricas que definen
V en las implı́citas de U , ası́ tendrı́amos el sistema

λ − 2(λ − µ) = 0
λ − (λ + µ) − µ = 0
58 CAPÍTULO 2. ÁLGEBRA

que es lo mismo que

λ − 2(λ − µ) = 0
µ = 0

de donde se tiene que λ = µ = 0 y sustituyendo en las ecuaciones paramétricas tenemos

que U ∩ V = 0.

Tenemos entonces que dim U ∩ V = 0 y estudiamos ahora la dimensión de U +V . Tenemos

que U + V = L[u1 , u2 , u3 , v1 , v2 ], para conocer la dimensión de este espacio calculamos el
rango de la matriz

 
0 1 0 0 0

 1 0 1 1/2 0 


 1 0 0 1/2 1 

 1 1 1 1 0 
0 2 1 −1 1

Para calcular el rango de esta matriz empezamos por su determinante, resolviendo primero
por la primera fila tenemos:

0 1 0 0 0
1 1 1/2 0
1 0 1 1/2 0
1 0 1/2 1
1 0 0 1/2 1 = −
1 1 1 0
1 1 1 1 0
0 1 −1 1
0 2 1 −1 1
0 0 1/2 0
1 0 1/2 1
= −
1 1 1 0
0 1 −1 1
1 0 1
1
= − 1 1 0
2
0 1 1
= 1 6= 0

donde en el segunda paso hemos sumado a la primera fila, la tercera multiplicada por −1
y después hemos resuelto de nuevo por la primera fila. Ası́ tenemos que la matriz tiene
rango 5. Ésta es la dimensión del espacio U + V .

También podemos aplicar la fórmula que relaciona las dimensiones de dos subespacios
vectoriales con las de su suma e intersección y tendrı́amos:

dim (U + V ) + dim (U ∩ V ) = dim U + dim V

dim (U + V ) + 0 = 3 + 2.
2.10. EJERCICIOS 59

Sistemas de ecuaciones

1. Hallar la solución del sistema


 3x + 2y + z = 1
5x + 3y + 4z = 2
x+y−z = 1


Solución: (Método de Gauss) Escribimos la matriz ampliada del sistema permu-

tando la primera y tercera filas
 
3 2 1 1
 5 3 4 2 
1 1 −1 1

Mediante el proceso de hacer ceros en el triángulo inferior izquierdo de la matriz

tenemos
     
3 2 1 1 3 2 1 1 3 2 1 1
 5 3 4 2  →  0 −1 4 −2  →  0 −1 4 −2 
1 1 −1 1 0 −2 9 −3 0 0 1 1

y el sistema inicial es equivalente al sistema


 x+y−z = 1
−y + 4z = −2
z = 1


que se resuelve escalonadamente y tiene como solución

(−4, 6, 1)

2. Discutir el sistema según los valores de a


 x − 2y + 2z = 2
2x + y + 3z = 2
5x + ay + 8z = 6


Solución: Calculamos el determinante de la matriz A del sistema e igualamos a

cero:
1 −2 2
|A| = 2 1 3 =a=0
5 a 8
Tenemos las siguientes alternativas

a) a 6= 0, entonces rg(A) = 3 y entonces también rg(A∗ ) = 3. Por el Teorema de

Rouché-Frobenius, el sistema es compatible determinado. Podemos calcular la
60 CAPÍTULO 2. ÁLGEBRA

solución aplicando la Regla de Cramer

2 −2 2
2 1 3
6 a 8 16 + 4a − 36 − (12 + 6a − 32)
x = = = −2
|A| a
1 2 2
2 2 3
5 6 8 16 + 30 + 24 − (20 + 18 + 32)
y = = =0
|A| a
1 −2 2
2 1 2
5 a 6 6 − 20 + 4a − (10 + 2a − 24)
z = = =2
|A| a

y la solución del sistema para a 6= 0 es (−2, 0, 2)

b) a = 0 Estudiamos la matriz ampliada del sistema, teniendo en cuenta que para
este valor de a se tiene rg(A) = 2
 
1 −2 2 2
A∗ =  2 1 3 2 
5 0 8 6

y vemos que el determinante de la submatriz formada por las tres últimas

columnas es cero. Resulta entonces que rg(A∗ ) = 2 y el sistema es compatible
indeterminado. Podemos prescindir de una de las ecuaciones del sistema, y el
sistema que tenemos que resolver es:
 
x − 2y + 2z = 2 x − 2y = 2 − 2z
⇔
2x + y + 3z = 2 2x + y = 2 − 3z

que podemos resolver por la regla de Cramer

1 2 − 2z
2 2 − 3z −8z + 6
x = =
5 5
1 2 − 2z
2 2 − 3z z−2
y = =
5 5
y si tomamos z como parámetro podemos expresar el conjunto de soluciones
 
−8λ + 6 λ − 2
, ,λ
5 5

3. Discutir el sistema 
 ax + y + bz = 1
x + ay + z = 0
x + y + bz = 1

2.10. EJERCICIOS 61

Calculamos el determinante de la matriz A del sistema e igualamos a cero:

a 1 b
|A| = 1 a 1 = a2 b + b + 1 − (ab + a + b) = 0
1 1 b
= ab(a − 1) − (a − 1)
= (ab − 1)(a − 1) = 0

Tenemos las siguientes alternativas

a) a 6= 1, ab 6= 1 ⇒ rg(A) = rg(A∗ ) = 3, se trata de u sistema compatible

determinado y cuya solución es:

1 1 b
0 a 1
1 1 b ab + 1 − (ab + 1)
x = = =0
|A| (ab − 1)(a − 1)
1 2 2
2 2 3
5 6 8 1 + b − (a + b) 1−a
y = = =
|A| (ab − 1)(a − 1) (ab − 1)(a − 1)
1 −2 2
2 1 2
5 a 6 a2 − a
z = =
|A| (ab − 1)(a − 1)

b) a = 1, b 6= 1 ⇒ rg(A) = 2, estudiamos la matriz ampliada

 
1 1 b 1
A∗ =  1 1 1 0 
1 1 b 1

donde se ve claramente que hay que estudiar la submatriz formada por las tres
últimas columnas. Esta submatriz tiene determinante cero, por lo que el vector
columna añadido a la matriz A es combinación lineal de los otros vectores
columna, y entonces tenemos rg(A∗ ) = 2 y se trata de un sistema compatible
indeterminado.
Para calcular las soluciones basta tener en cuenta que la tercera ecuación es
exactamente igual que la primera, y entonces tenemos el sistema

x + y + bz = 1
x+y+z = 0
de donde tenemos, restando la segunda ecuación a la primera

1
(b − 1)z = 1 ⇒ z =
b−1
62 CAPÍTULO 2. ÁLGEBRA

y sustituyendo en la segunda
1
x= −y
1−b
el conjunto de soluciones es, tomando como parámetro libre λ = y:
 
1 1
− λ, λ,
1−b b−1

c) a = 1, b = 1 ⇒ rg(A) = 1, rg(A∗ ) = 2, en este caso el sistema es incompatible.

d ) a 6= 1, ab = 1 ⇒ rg(A) = 2, estudiamos la matriz ampliada, más concretamente
la submatriz que resulta al suprimir la tercera columna
 
a 1 b 1 a 1 1
∗
A =  1 a 1 0  −→ 1 a 0 = a2 − a 6= 0
1 1 b 1 1 1 1

y por lo tanto la matriz ampliada tiene rango 3 y el sistema es incompatible.

2.11. DIAGONALIZACIÓN DE APLICACIONES LINEALES SIMÉTRICAS 63

2.11. Diagonalización de aplicaciones lineales

simétricas
A lo largo de esta sección utilizaremos notación matricial (letras mayúsculas) para los
vectores. Esto lo haremos porque en algunas demostraciones nos interesa distinguir al
vector escrito en forma de fila o columna.

Vectores y valores propios

Sea V un K-espacio vectorial y
A:V →V
una aplicación lineal. Un elemento V ∈ V es un vector propio de A si existe un número
λ ∈ K tal que
AV = λV
en este caso, dado que λ está determinado de forma única (si V 6= 0, λ1 V = λ2 V implica
que λ1 = λ2 ) decimos que λ es un valor propio.

Si A es una matriz real n × n, según la definición de vector propio, éste es un vector

columna X de Rn de manera que existe λ ∈ R tal que

AX = λX

Ejemplo 2.8. Sea A : V → V una aplicación lineal, y V un vector propio de A, entonces

para cualquier escalar µ 6= 0, µV también es un vector propio de A, asociado al mismo
valor propio.

Teorema 2.17. Sea V un espacio vectorial y sea V → V una aplicación lineal. Sea λ ∈ K
y sea Vλ el subespacio de V generado por todos los vectores propios de A cuyo valor propio
es λ. Entonces todo elemento no nulo de Vλ es un vector propio de A que tiene a λ como
valor propio.

Demostración. Se deja como ejercicio.

Teorema 2.18. Sea V un espacio vectorial y sea A : V → V una aplicación lineal.

Sean v1 , . . . , vm vectores propios de A. Sean λ1 , . . . , λm los valores propios respectivos. Si
λi 6= λj para i 6= j, entonces los vectores v1 , . . . , vm son linealmente independientes.

Demostración. Los vi son distintos de cero por ser autovectores. Supongamos que la afir-
mación no es cierta. En ese caso los vectores v1 , . . . , vm serı́an linealmente dependientes.
64 CAPÍTULO 2. ÁLGEBRA

Podemos suponer que el primero de ellos es combinación lineal de los otros; esto es, existen
a2 , · · · , am tal que

0 6= v1 = a2 v2 + · · · + am vm (2.11)

y además podemos suponer que los vectores v2 , . . . , vm son linealmente independientes, de

lo contrario nos podrı́amos quedar con mayor el subconjunto de ellos que fuese linealmente
independiente. Entonces

λ1 v1 = λ1 a2 v2 + · · · + λ1 am vm

Por otra parte, se tiene

A(v1 ) = λ1 v1 = a2 λ2 v2 + · · · am λm vm

Restando las dos expresiones que tenemos del vector λ1 v1 resulta

0 = a2 (λ1 − λ2 )v2 + · · · + am (λ1 − λm )vm

y como los vectores v2 , . . . , vm son linealmente independientes tendrı́an que ser ai (λ1 −
λi ) = 0 para i = 2, . . . , m, y como los autovalores son distintos, serı́a a2 = · · · = am = 0,
y de la igualdad (2.11) serı́a v1 = 0.

Observación 2.3. Si V es un espacio vectorial de dimensión n y A : V → V una aplica-

ción lineal con n vectores propios v1 , . . . , vn con valores propios λ1 , . . . , λn son distintos
entre sı́. Entonces {v1 , . . . , vn } es una base de V .

El polinomio caracterı́stico
En este apartado estudiaremos la forma de obtener los valores propios de una matriz
Teorema 2.19. Sea V un espacio vectorial de dimensión finita y sea λ un escalar. Sea
A : V → V una aplicación lineal. Entonces λ es un valor propio de A si, y sólo si, A − λI
es singular.

Definición 2.24. Sea A = (aij ) una matriz n×n. Se define el polinomio caracterı́stico
de A y lo denotamos PA al determinante

PA (t) = det (tI − A)

2.11. DIAGONALIZACIÓN DE APLICACIONES LINEALES SIMÉTRICAS 65

Teorema 2.20. Sea A una matriz n × n. Un escalar λ es un valor propio de A si, y sólo
si, λ es una raı́z del polinomio caracterı́stico de A.

Teorema 2.21. Sean A y B dos matrices n × n y supongamos que B es invertible.

Entonces el polinomio caracterı́stico de A es el mismo que el de B −1 AB.

Demostración. Tenemos

det(tI − A) = det(B −1 (tI − A)B) = det(tB −1 B − B −1 AB)

= det(tI − B −1 AB).

En el contexto que nos estamos moviendo vamos a enunciar aunque sin demostración el
Teorema Fundamental del Álgebra, que citaremos para asegurar la existencia de autovalo-
res y, más concretamente en el estudio de las matrices simétricas, donde éstos son siempre
reales.

Teorema 2.22 (Teorema Fundamental del Álgebra). Todo polinomio no constante con
coeficientes en C tiene una raı́z en C.

Valores propios y vectores propios de matrices reales simétricas

Matrices Ortogonales

Definición 2.25. Una matriz cuadrada M con entradas reales se dice que es ortogonal
si
MMT = I
donde I representa la matriz identidad.
Observación 2.4. Las matrices ortogonales representan aplicaciones lineales isométri-
cas dentro del espacio vectorial Rn . Esto es, si consideramos la métrica inducida por el
producto escalar, de manera que la distancia entre dos vectores es

d(u, v) = ku − vk

entonces una matriz ortogonal determina una aplicación lineal que preserva esa métrica
y, por ser lineal, fija el origen.
66 CAPÍTULO 2. ÁLGEBRA

Proposición 2.6. Sea A una matriz cuadrada con entradas reales, entonces son equiva-
lentes:

1. A es ortogonal.

2. det A = ±1.

3. Tanto los vectores fila como los vectores columna forman una base ortonormal.

Demostración. La equivalencia entre 1) y 2) es evidente

det (AAT ) = det A · det AT = det I = 1

Para ver la equivalencia entre 1) y 3) escribimos la matriz A

 
a11 a12 · · · a1n
 a21 a22 · · · a2n 

A =  .. ..  = A1 A2 · · · An
 
.. ..
 . . . . 
an1 an2 · · · ann

entonces la matriz identidad se tiene haciendo los distintos productos escalares entre estos
vectores, concretamente
 
A1

 A2 

I = AAT = A1 A2 · · · An  .. 
 . 
An

de donde se tiene que Ai · Ai = 1 = a21i + a22i + · · · + a2ni y Ai · Aj = 0, i 6= j.

Teorema 2.23. Sea A una matriz real simétrica de orden n × n. Entonces los autovalores
de A son todos reales.

Demostración. Sea λ un autovalor de A, que existe por el Teorema Fundamental del

Álgebra, entonces para algún vector 0 6= V ∈ Cn se tiene

AV = λV (2.12)

conjugando los dos miembros de 2.12 tenemos (1 )

AV = λV (2.13)
1
Teniendo en cuenta que zw = z w para z, w ∈ C.
2.11. DIAGONALIZACIÓN DE APLICACIONES LINEALES SIMÉTRICAS 67

y como A tiene entradas reales se tiene A = A; es decir:

AV = λ V . (2.14)

Por otra parte, si trasponemos 2.14 tenemos

T T
V AT = V λ (2.15)

y como A es simétrica se tiene

T T
V A=V λ (2.16)
multiplicando los dos miembros de esta igualdad a la derecha por v y los dos miembros
de 2.12 por vT tenemos
T T
V AV = V λV
T T
V AV = V λV

y entonces
T T T
0 = V λV − V λV = (λ − λ)V V = (λ − λ)(|v1 |2 + · · · + |vn |2 )

y como V 6= 0 tiene que ser λ = λ; es decir, λ ∈ R.

Observación 2.5. De las ecuaciones 2.12 y 2.13 se deduce que los valores propios de
matrices reales vienen en pares conjugados.

Teorema 2.24. Los autovectores de una matriz simétrica A correspondientes a autova-

lores distintos son ortogonales entres ellos.

Demostración. Sean λ1 y λ2 dos autovalores distintos de la matriz A, y sean los autovec-

tores asociados V1 y V2 , entonces de las igualdades AVi = λi Vi , i = 1, 2 se tiene

V2T AV1 = V2T λ1 V1

V1T AV2 = V1T λ2 V2

y restando miembro a miembro, resulta de la simetrı́a de A que

(λ1 − λ2 )V1T V2 = 0.

de donde resulta que V1 y V2 son ortogonales.

68 CAPÍTULO 2. ÁLGEBRA

Proposición 2.7. Si todos los autovalores de una matriz A son distintos, entonces la
matriz X cuyas columnas son los autovectores asociados a esos autovalores, es una matriz
ortogonal.

Demostración. Una matriz es ortogonal si sus filas o columnas constituyen una base orto-
normal. Equivalentemente X es ortogonal si X T = X −1 . El resultado se tiene del teorema
anterior al considerar la base formada por autovectores, de manera que éstos sean de
módulo unidad, ya que entonces los vectores columna de X forman una base ortonor-
mal.

Teorema 2.25. Sea A una matriz simétrica y λ una raı́z del polinomio caracterı́stico
PA (t) con multiplicidad m ≥ 2. Entonces existen m autovectores ortonormales asociados
a λ.

Demostración. Por ser λ raı́z de PA (t), existe un autovector X1 asociado a este valor. Po-
demos formar con otros vectores Y2 , . . . , Yn una base ortonormal de Rn (ver demostración
del teorema 2.13). Consideremos la matriz

M = [X1 Y2 . . . Yn ] = [X1 Y ]

cuyas columnas son las coordenadas de estos vectores.

Se tiene entonces
     
T λX1T X1 X1T AY λ λX1T Y λ 0
M AM = = =
λY T X1 Y T AY λY T X1 Y T AY 0 Y T AY

teniendo en cuenta que X1 es ortogonal a los vectores Yi .

Las matrices A y M −1 AM tienen el mismo polinomio caracterı́stco (teorema 2.21), y

como M es ortogonal se tiene que M −1 = M T , entonces

det (M T AM − tIn ) = (λ − t)det (Y T AY − tIn−1 ).

Si la raı́z λ tiene multiplicidad m ≥ 2 tendremos

det (Y T AY − λIn−1 ) = 0

lo que significa que para el valor λ la matriz rg (M −1 AM ) ≤ 2, esto es, el espacio Sλ

de soluciones del sistema homogéneo correspondiente tiene dimensión mayor o igual a 2.
A este subespacio Sλ pertenece el vector X1 . Entonces podemos encontrar otro vector
X2 ∈ Sλ ortogonal a X1 . Igual que antes, podemos encontrar n − 2 vectores Yi tales que:

{X1 , X2 , Y3 , . . . , Yn }
2.11. DIAGONALIZACIÓN DE APLICACIONES LINEALES SIMÉTRICAS 69

forma una base ortonormal de Rn . Se repite el argumento hasta llegar una base

{X1 , . . . , Xm , Ym+1 , . . . , Yn }

donde los vectores X1 , . . . , Xm forman una base ortonormal del espacio Sλ .

Diagonalización de una aplicación lineal simétrica

Si A es una matriz (no necesariamente simétrica) cuyos autovalores son λ1 , . . . , λn y sean
vectores propios X1 , . . . , Xn , tales que kX1 k = 1, entonces:

Ai Xi = λi Xi , i = 1, . . . , n

consideremos la matriz
X = [X1 · · · Xn ]

cuyas columnas representan los vectores propios normalizados de A. Se tienen entonces

AX = [AX1 · · · AXn ]
= [λn X1 · · · λn Xn ]
= XD

donde  
λ1 0 · · · 0
 0 λ2 · · · 0 
D=
 
.. .. . . .
. ..

 . . 
0 0 0 λn

es la matriz diagonal formada por los autovalores (siguiendo el mismo orden en que se
han puesto los autovectores).

Si la matriz X tiene inversa, entonces

D = X −1 AX. (2.17)

Si la matriz A es simétrica, los autovectores Xi normalizados elegidos antes son ortonor-

males y entonces X es una matriz ortogonal, entonces X −1 = X T y 2.17 queda en este
caso
D = X T AX.
70 CAPÍTULO 2. ÁLGEBRA

Formas Cuadráticas

Definición 2.26. Sea A una matriz simétrica y X ∈ Rn un vector columna. La función

f : Rn → R
f (X) = X t AX
se llama forma cuadrática asociada a la matriz A. X ∈ Rn

Ejemplo 2.9. Sea  

a b
b c
una matriz 2×2 con entradas reales. Entonces la expresión de la forma cuadrática asociada
es:   

a b x
x y = ax2 + 2bxy + dy 2
b c y

Teorema 2.26. Sea A una matriz simétrica real y sea f (X) = X t AX su forma cuadrática
asociada. Sea P un punto sobre la esfera unitaria tal que f (P ) es un máximo para f sobre
la esfera. Entonces P es un vector propio para A. Dicho de otro modo, existe un número
λ tal que AP = λA.

Corolario 2.5. El valor máximo de f sobre la esfera de radio uno es igual al valor propio
más grande de A.
Capı́tulo 3

Geometrı́a

3.1. Espacio afı́n

La estructura de espacio afı́n está asociada a la de espacio vectorial de manera natural. En
el espacio vectorial, hemos visto como definir un producto escalar, lo que permite hablar
de medidas de águlos y distancias. Los subespacios de un espacio vectorial son variedades
lineales: rectas, planos, etc que pasan por el origen por ser espacios vectoriales.

El espacio afı́n permite ver un espacio vectorial como espacio de puntos además de su
estructura vectorial, y de este modo podremos considerar copias de subespacios lineales
que no necesariamente pasan por el origen.

El desarrollo formal que haremos a continuación es para R3 pero su generalización a otras

dimensiones es casi inmediata.

Consideremos el conjunto R3 . Cada elemento de este conjunto es una terna ordenada de

números reales.
R3 = {(x1 , x2 , x3 ) : xi ∈ R}
los elementos de R3 son puntos y se denotan con letra mayúscula P = (p1 , p2 , p3 ). Cada
punto se identifica con un vector de R3 .

Más generalmente, dos puntos P y Q determinan un vector del espacio vectorial R3 :

P Q = v = (q1 − p1 , q2 − p2 , q3 − p3 )

Definición 3.1. Sean u y v dos vectores de R3 . Diremos que u y v son equipolentes si

tienen el mismo módulo, dirección y sentido.

Definición 3.2. La relación de equipolencia entre vectores es una relación de equivalencia.

Un clase de equivalencia de esta relación se llama vector libre.

El conjunto de vectores equipolentes construidos a partir del espacio R3 es un espacio

vectorial isomorfo a R3 y lo notaremos V (R3 ).

71
72 CAPÍTULO 3. GEOMETRÍA

Espacio afı́n A3
Lo dicho antes permite considerar en el espacio R3 simultaneamente la estructura de
espacio vectorial con la de espacio de puntos de forma que ambas interactúan, como
veremos a continuación.
Definición 3.3. El espacio afı́n A3 asociado al espacio V (R3 ) es el espacio R3 junto con
la aplicación
ϕ : R3 × R3 → V (R3 )
tal que ϕ(A, B) es el vector libre representado por el vector AB. Además la aplicación ϕ
verifica:

∀A ∈ R3 , ∀v ∈ V (R3 ) exist un único B ∈ R3 tal que AB = v.

∀A, B, C ∈ R3 se tiene AC = AB + BC.

Sistema de referencia afı́n en R3 .

Un sistema de referencia afı́n de R3 lo constituyen una cuaterna de puntos {O, A, B, C}
de manera que los vectores
OA, OB, OC
formen una base de V (R3 ).

Un punto cualquiera P ∈ R3 determina de manera unı́voca el vector

p = OP = p1 u1 + p2 u2 + p3 u3
donde
OA = u1 , OB = u2 , OC = u3

Equivalentemente podemos escribir el sistema de referencia anterior como

{O; u1 , u2 , u3 }

Ecuaciones de la recta
Sea r una recta en el espacio. Un vector v es paralelo a la recta r si existen dos puntos A
y B de r tal que AB es equipolente a v. Diremos entonces que v es un vector director de
la recta r.

Consideremos un punto A del espacio afı́n R3 y v un vector libre no nulo. Hay una única
recta que pasa por A y tiene como vector director v.

A continuación vamos a obtener una condición necesaria y suficiente para que un punto
X pertenezca a la recta que pasa por A y tiene vector de dirección v. A partir de la
definición de vector director de la recta podemos establecer la siguiente
Afirmación 3.1. X pertenece a r definida por el punto A y con vector director v si y
sólo si existe un número λ ∈ R tal que AX = λv.
3.1. ESPACIO AFÍN 73

Tomemos un sistema de referencia {O; u1 , u2 , u3 } y sean a = OA y x = OX los vectores

de posición de los puntos A y X respecto del origen O del sistema de referencia. Entonces
se tiene
OA + AX = OX ⇐⇒ AX = OX − OA = x − a = λv
y finalmente la condición de pertenencia se puede escribir

x = a + λv (3.1)

Ahora, si (x, y, z) son las coordenadas del vector x en la base {u1 , u2 , u3 }, o lo que
es lo mismo, las coordenadas del punto X en el sistema de referencia {O; u1 , u2 , u3 }.
Análogamente, sean (a1 , a2 , a3 ) las coordenadas del vector a o del punto A y (v1 , v2 , v3 )
las coordenadas del vector v, la ecuaión 3.1 puede escribirse

(x, y, z) = (a1 , a2 , a3 ) + λ(v1 , v2 , v3 ) (3.2)

que a su vez podemos desglosar como

x = a1 + λv1
y = a2 + λv2
z = a3 + λv3

Estas ecuaciones se llaman ecuaciones paramétricas de la recta r, que está definida por el
punto A y el vector v.

Si al sistema de ecuaciones paramétricas anterior, lo consideramos un sistema de ecuacio-

nes lineales en la incógnita λ; es decir

λv1 = x − a1
λv2 = y − a2
λv3 = z − a3

entonces el sistema tiene solución cuando

   
v1 v1 x − a1
rango  v2  = rango  v2 y − a2  = 1
v3 v3 z − a3

ya que se considera el vector v no nulo. Esto quiere decir que las dos columnas de la
segunda matriz son proporcionales, lo que podemos expresar como
x − a1 y − a2 z − a3
= =
v1 v2 v3
Observación 3.1. Si una coordenada del vector de dirección v de la recta es igual a cero,
la expresión anterior tiene sentido entendiendo que expresa que los vectores (x − a1 , y −
a2 , z − a3 ) y (v1 , v2 , v3 ) son proporcionales, pero no es correcto escribir en una ecuación
un denominador igual a cero. Si, por ejemplo tenemos v = (v1 , 0, v3 ), la recta r está
74 CAPÍTULO 3. GEOMETRÍA

contenida en el plano y − a2 = 0, ya que el vector v es paralelo a este plano y la recta

pasa por A. Por otra parte la ecuación restringida a las variables x y z es la ecuación
(continua) de una recta en el plano XZ, por lo que considerada como ecuación en R3
es la ecuación de un plano perpendicular al plano XZ. Resulta entonces que la ecuación
continua de la recta r se transforma en
x − a1 z − a3
= ; y − a2 = 0
v1 v3
El mismo tipo de razonamiento podemos hacer si tenemos dos coordenadas del vector de
dirección iguales a cero.

Ecuaciones del plano

Un plano en el espacio vectorial R3 es un subespacio vectorial de dimensión 2. Si en espacio
afı́n R3 tomamos como sistema de referencia {O; e1 , e2 , e3 } de manera que el punto O sea
(0, 0, 0) y {e1 , e2 , e3 } la base canónica, los puntos P = (p1 , p2 , p3 ) del espacio afı́n R3 se
identifican con los vectores del espacio vectorial R3 que tienen las mismas coordenadas.
Un subespacio vectorial de dimensión 2 se identifica con un plano del espacio afı́n que
pasa por el origen.

La definición de espacio afı́n permite considerar en el espacio R3 planos en situación

más general, es decir, que no necesariamente pasan por el origen. El plano queda ası́
determinado si consideramos un subespacio vectorial de dimensión 2 y un punto del espacio
afı́n por el que pasa.

Sean v y w dos vectores de R3 linealmente independientes y sea M un punto de R3 . Estos

datos determinan un único plano π. Veamos que condiciones tiene que verifcar un punto
X para pertenecer al plano π.

Si X ∈ π se tiene de lo dicho antes

M X = λv + µw (3.3)

y recı́procamente, es decir, si se tiene esta condición entonces X ∈ π. Si llamamos x y m

a los vectores de posición de los puntos X y M se tiene que M X = x − m y la ecuación
3.3 puede escribirse

x = m + λv + µw (3.4)

que se llama ecuación vectorial del plano. Al igual que en el caso de la recta podemos
desglosar esta ecuación por componentes

x = m1 + λv1 + µw1
y = m2 + λv2 + µw2
z = m3 + λv3 + µw3
3.1. ESPACIO AFÍN 75

y estas son las ecuaciones paramétricas del plano π. Estas ecuaciones pueden verse como
un sistema en las incógnitas λ y µ. Escribimos el sistema de tres ecuaciones

λv1 + µw1 = x − m1
λv2 + µw2 = y − m2
λv3 + µw3 = x − m3

entonces el sistema tiene solución cuando

   
v1 w 1 v1 w1 x − m1
rango  v2 w2  = rango  v2 w2 y − m2  (3.5)
v3 w3 v3 w3 z − m3

La primera matriz tiene rango 2 ya que los vectores v y w son linealmente independientes.
Entonces el sistema tiene solución si la otra matriz tiene rango 2, para lo cuál se tiene
que tener
v1 w1 x − m1
v2 w2 y − m2 = 0
v3 w3 z − m3
y tenemos ası́ una ecuación del plano. Podemos transformar más esta ecuación. Si desa-
rrollamos el determinante por la última columna tenemos

v2 w2 v1 w1 v1 w1
(x − m1 ) − (y − m2 ) + (z − m3 ) =0 (3.6)
v3 w3 v3 w3 v2 w2

y si llamamos

v2 w2 v w1 v1 w1
a= , b=− 1 , c= , d = −m1 a − m2 b − m3 c
v3 w3 v3 w 3 v2 w2

podemos escribir 3.6 como

ax + by + cz + d = 0
que es la ecuación implı́cita del plano.

Observación 3.2. En la ecuación anterior no puede ser (a, b, c) = (0, 0, 0) porque si fuese
ası́ serı́a el determinante de la primera matriz de 3.5 igual a cero, ya que a, b y c son todos
los menores de orden dos de esta matriz.

Ejemplo 3.1. Hallar la ecuación del plano que pasa por el punto M = (1, 0, −1) y es
paralelo a los vectores v = (−1, 3, 2) y w = (1, 0, 5).

Solución: Las coordenadas se suponen con respecto a un sistema de referencia afı́n

{O; u1 , u2 , u3 } y en este sistema se tiene

v = −u1 + 3u2 + 2u, w = u1 + 5u3

76 CAPÍTULO 3. GEOMETRÍA

y si x es el vector de posición de un punto X del plano la ecuación de éste es

x = m + λv + µw

que se puede desarrollar, obteniendo las ecuaciones paramétricas

x = 1−λ+µ
y = 3λ
z = −1 + 2λ + 5µ

y la ecuación implı́cita la obtenemos de

−1 1 x − 1
3 0 y =0
2 5 z+1

que desarrollando se tiene

15x + 7y − 3z − 18 = 0
3.2. PROBLEMAS DE INCIDENCIA 77

3.2. Problemas de incidencia

Recta incidente con dos puntos
Consideremos dos puntos distintos A = (a1 , a2 , a3 ) y B = (b1 , b2 , b3 ). El vector AB es un
vector de dirección de la recta que pasa por A y B. Las coordenadas de este vector son
AB = (b1 − a1 , b2 − a2 , b3 − a3 )

y una ecuación continua de la recta será

x − a1 y − a2 z − a3
= =
b 1 − a1 b 2 − a2 b 3 − a3
o también las ecuaciones paramétricas
x = a1 + λ(b1 − a1 )
y = a2 + λ(b2 − a2 )
z = a3 + λ(b3 − a3 )

Rectas paralelas
Consideremos las rectas r y r0 de ecuaciones
r : x = a + λv r0 : x = b + µw
donde a es el vector de posición del punto A y b es el vector de posición del punto B.

Si las r y r0 son paralelas entonces los vectores v y w son proporcionales. En este caso
cualquier recta que pase por un punto H = (h1 , h2 , h3 ) y paralela a r (y por lo tanto a r0 )
tiene como ecuación
x − h1 x − h2 x − h3
= =
v1 v2 v3

Plano incidente con tres puntos

Dados los puntos M = (m1 , m2 , m3 ), N = (n1 , n2 , n3 ) y Q = (q1 , q2 , q3 ) no alineados,
hallaremos la ecuación del plano que determinan. Este problema podemos reducirlo al
anterior considerando los vectores linealmente independientes M N y M Q, que tienen
coordenadas
M N = (n1 − m1 , n2 − m2 , n3 − m3 )
M Q = (q1 − m1 , q2 − m2 , q3 − m3 )

y tenemos la ecuación
n1 − m1 q1 − m1 x − m1
n2 − m2 q2 − m2 y − m2 =0
n3 − m3 q3 − m3 z − m3
78 CAPÍTULO 3. GEOMETRÍA

Planos paralelos
Supongamos que tenemos dos planos π y π 0 paralelos. Podemos suponer dos vectores
v = (v1 , v2 , v3 ) y w = (w1 , w2 , w3 ) linealmente independientes que son paralelos a los dos
planos. Para tener definidos a π y π 0 necesitamos además dos puntos M = (m1 , m2 , m3 ) ∈
π y N = (n1 , n2 , n3 ) ∈ π 0 .

Operando como antes tenemos las ecuaciones implı́citas de los dos planos. Si llamamos
v2 w 2 v w1 v1 w1
a= , b=− 1 , c=
v3 w 3 v3 w3 v2 w2
y
d = −m1 a − m2 b − m3 c, d0 = −n1 a − n2 b − n3 c
entonces las ecuaciones implı́citas de los planos son
π : ax + by + cz + d = 0, π 0 : ax + by + cz + d0 = 0

Con el razonamiento anterior se puede hallar la ecuación del plano π 0 paralelo a uno dado
π : ax + by + cz + d = 0, que pasa por el punto N = (n1 , n2 , n3 ).

La ecuación de π 0 es ax + by + cz + d0 = 0, donde d0 = −n1 a − n2 b − n3 c, sustituyendo en

la ecuación de π tenemos
ax + by + cz + d0 = 0
ax + by + cz − n1 a − n2 b − n3 c = 0
a(x − n1 ) + b(y − n2 ) + c(z − n3 ) = 0

y esta última expresión refleja una situación geométrica: las coordenadas de un punto
X = (x, y, z) arbitrario son las de un punto del plano si y aólo si el vector N X es
ortogonal al vector (a, b, c).

Recta y plano paralelos

Consideremos el plano π de ecuación ax + by + cz + d = 0 y la recta r de ecuación
x = p + tv. Entonces la recta r es paralela a π si existe una recta r0 contenida en π que
es paralela a r.

Consideremos dos puntos distintos de la recta r0 , M = (m1 , m2 , m3 ) y N = (n1 , n2 , n3 ).

Los dos puntos pertenecen a π y además el vector M N 6= 0 es paralelo al vector v. De la
pertenencia al plano se tiene
am1 + bm2 + cm3 + d = 0, an1 + bn2 + cn3 + d = 0
ahora restando la primera ecuación de la segunda tenemos
a(n1 − m1 ) + b(n2 − m2 ) + c(n3 − m3 ) = 0 (3.7)
3.2. PROBLEMAS DE INCIDENCIA 79

y de la segunda condición: el vector M N = (n1 − m1 , n2 − m2 , n3 − m3 ) es paralelo al

vector v = (v1 , v2 , v3 ), se tiene entonces:

n1 − m1 n2 − m2 n3 − m3
= =
v1 v2 v3

si a esta razón de proporcionalidad la llamamos k 6= 0 tenemos

n1 − m1 = kv1 , n2 − m2 = kv2 , n3 − m3 = kv3

y sustituyendo en 3.7 se tiene

akv1 + bkv2 + ckv3 = k(av1 + bv2 + cv3 ) = 0

de donde se tiene que

av1 + bv2 + cv3 = 0
y esta es la condición de paralelismo entre recta y plano. Esta ecuación es susceptible de
interpretar geométricamente en términos del producto escalar, o sea, de la perpendicula-
ridad.

Intersección de planos
Consideremos los planos π y π 0 de ecuaciones ax + by + cz + d = 0 y a0 x + b0 y + c0 z + d0 = 0
respectivamente. Supongamos que los dos planos no son paralelos, esto implica que
 
a b c
rango =2 (3.8)
a0 b 0 c 0

que es lo mismo que afirmar que el sistema formado por las ecuaciones de los dos planos
tiene infinitas soluciones. Estas soluciones son las coordenadas de los puntos de la recta
intersección de los dos planos.

Observación 3.3. Las dos ecuaciones anteriores definen la recta de intersección de los
planos π y π 0 . Esto es, el par de ecuaciones equivale a cualquier ecuación de la recta π ∩π 0 ,
como veremos a continuación.

De la condición de 3.8 deducimos que existe un menor de orden 2 distinto de cero, supon-
gamos
a b
6= 0
a0 b 0
Podemos escribir el sistema

ax + by = −cz − d
a0 x + b0 y = −c0 s − d0
80 CAPÍTULO 3. GEOMETRÍA

y aplicando la regla de Cramer encontramos la solución

−cz − d b c b d b
− z− 0 0
−c0 z − d0 b0 c 0 b0 d b
x= = = Az + B
a b a b
a0 b 0 a0 b 0

a −cz − d a c d b
− 0 z −
a0 −c0 z − d0 0
a c d0 b0
y= = = Cz + D
a b a b
a0 b 0 a0 b 0

Entonces resultan las ecuaciones paramétricas de la recta


 x = B + At
y = D + Ct
z = t


y es la recta que pasa por el punto de coordenadas (B, D, 0) y tiene la dirección del vector
(A, C, 1).

Haz de planos
Dada una recta queremos hallar la ecuación paramétrica de todos los planos que contienen
a esa recta. Este conjunto de planos se llama haz de planos y la recta dada se llama arista.
Esta ecuación puede ser útil para la resolución de algunos problemas.

Consideremos, por ejemplo la recta

ax + by + cz + d = 0
a x + b0 y + c 0 z + d 0 = 0
0

ahora hacemos la combinación lineal de las ecuaciones de esos planos con coeficientes s y
t:
s(ax + by + cz + d) + t(a0 x + b0 y + c0 z + d0 ) = 0
que podemos escribir

(sa + ta0 )x + (sb + tb0 )y + (sc + tc0 )z + (sd + td0 ) = 0

y para (s, t) 6= (0, 0) tenemos la ecuación de un plano. Veamos que este plano pasa por la
recta dada.

Si P = (p1 , p2 , p3 ) es un punto de la recta, tenemos que probar que P pertenece al plano.

En primer lugar, por pertenecer P a la recta, se tiene que verificar

ap1 + bp2 + cp3 + d = 0, a0 p 1 + b 0 p 2 + c 0 p 3 + d 0 = 0

3.3. EL ESPACIO EUCLIDIANO TRIDIMENSIONAL 81

y ahora para ver que P está en el plano (s, t) arbitrario sustituimos

(sa + ta0 )p1 + (sb + tb0 )p2 + (sc + tc0 )p3 + (sd + td0 ) =
= s(ap1 + bp2 + cp3 + d) + t(a0 p1 + b0 p2 + c0 p3 + d0 ) = 0

y entonces P está en la recta, como es un punto arbitrario, se tiene que toda la recta
está contenida en el plano. También el par (s, t) es arbitrario, esto significa que todos los
planos del haz contienen a la recta.

3.3. El espacio euclidiano tridimensional

En el capı́tulo anterior se ha estudiado el espacio afı́n tridimensional y problemas que se re-
ferı́an a incidencia, intersección y paralelismo, pero no los problemas métricos: longitudes,
áeas, águlos. Si dotamos al espacio afı́n de un producto escalar, tenemos la herramienta
que necesitamos para hacer esos cálculos. El modelo que se construye en el espacio afı́n
R3 con el producto escalar habitual es el del espacio euclidiano.

3.4. Producto escalar

Si consideramos dos vectores libres a y b y consideramos dos representantes de ellos con
el mismo origen O, OA y OB. Entonces los dos vectores definen en el plano OAB dos
ángulos cuya suma es la de dos rectos. Llamaremos ángulo de los vectores a y b al menor
de ellos.

Definición 3.4. Llamamos producto escalar de dos vectores a y b al producto de los

módulos de los vectores por el coseno del ángulo que forman.

a · b = |a| |b| cos (a, b)

Proposición 3.1. El producto escalar cumple las siguientes propiedades:

1. a · b = b · a.

2. a · (b + c) = a · b + a · c.

3. (λa) · b = a · (λb) = λ(a · b).

4. a · a = |a|2 .

5. a · 0 = 0, donde 0 es el vector nulo.

6. Si a y b son distintos de cero, la condición necesaria y suficiente para que sean

perpendiculares (el ángulo que forman es de 90◦ ) es que a · b = 0.
82 CAPÍTULO 3. GEOMETRÍA

3.5. Espacio métrico

Una vez que hemos introducido el producto en el espacio vectorial R3 , y consideremos
una base {u, v, w} de manera que cada uno de los vectores tiene módulo 1 y entre ellos
son perpendiculares.

Sean ahora dos vectores a y b que con respecto a la base anterior tienen coordenadas
(a1 , a2 , a3 ) y (b1 , b2 , b3 ):

a = a1 u + a2 v + a3 w, b = b1 u + b2 v + b3 w

ahora calculamos el producto escalar teniendo en cuenta las propiedades de la proposición

anterior:

a · b = (a1 u + a2 v + a3 w) · (b1 u + b2 v + b3 w)
= a1 b1 (u · u) + a2 b2 (v · v) + a3 b3 (w · w) + · · ·

donde el resto de los sumandos se anulan porque cada uno de ellos lleva un producto
escalar de vectores de la base distintos, y por ser perpendiculares éste es cero.

Resulta ası́ la expresión del producto escalar en coordenadas

a · b = a1 b 1 + a2 b 2 + a3 b 3

(comparar con la sección 2.7).

Lema 3.1. El coseno del ángulo de dos vectores a y b de coordenadas (a1 , a2 , a3 ) y
(b1 , b2 , b3 ) respecto de una base ortonormal es:
a1 b 1 + a2 b 2 + a3 b 3
cos (a, b) = p 2 p
a1 + a22 + a23 b21 + b22 + b23

Demostración. Basta tener en cuenta que

a·b
cos (a, b) =
|a| · |b|

3.6. El espacio euclidiano

El espacio afı́n tridimensional R3 vamos a considerarlo construido sobre el espacio vecto-
rial dotado de producto escalar. Entonces podemos considerar un sistema de referencia
métrico u ortonormal que será {O; u1 , u2 , u3 } siendo la base {u1 , u2 , u3 } una base orto-
normal.
3.6. EL ESPACIO EUCLIDIANO 83

Definición 3.5. Dado el espacio euclidiano R3 con un sistema de referencia métrico, la

distancia d(A, B) entre dos puntos A y B es igual al módulo del vector AB con coordenadas
en este sistema.
p
d(A, B) = AB = (b1 − a1 )2 + (b2 − a2 )2 + (b3 − a3 )2

Supongamos que un plano tiene de ecuación con respecto del sistema de referencia métrico

ax + by + cz + d = 0

fijemos un punto P de ese plano cuyas coordenadas (p1 , p2 , p3 ) verifican la ecuación. En-
tonces
ap1 + bp2 + cp3 + d = 0
y restando la ecuación del plano

a(x − p1 ) + b(y − p2 ) + c(z − p3 ) = 0

y esto lo podemos interpretar del siguiente modo: si llamamos n = (a, b, c) este vector es
distinto de cero. El vector P X = (x − p1 , y − p2 , z − p3 ) es distinto de cero si P 6= X y
X es un punto arbitrario que satisface las coordenadas del plano. Entonces la ecuación
anterior nos dice que estos vectores son perpendiculares cualquiera que sea X.

Ángulo de dos rectas

Dos rectas no paralelas en el espacio pueden no tener ningún punto en común, estas son
rectas que se cruzan. Sin embargo, podemos hablar del ángulo de dos rectas que se cruzan.

La manera más natural de definir este ángulo será la del ángulo que forman dos rectas
paralelas a ambas que pasan por un punto dado.

Hemos trasladado de esta forma el problema de determinar el ángulo a determinar el

ángulo en caso de que las rectas se corten en el plano. Cuando tenemos dos rectas que
se cortan en el plano, se entiende por ángulo que forman al menor de los dos ángulos
suplementarios que forman. Sin embargo el ángulo que forman los vectores de dirección
puede ser el menor de las rectas o el mayor.

Entonces el ángulo que buscamos será el menor de entre los dos siguientes: que forman
sus vectores directores o su suplementario.

Si r y r0 son las rectas que tenemos y v y w sus vectores de dirección, para tomar el
ángulo menor de entre los dos anteriores, si nos referimos al coseno tenemos

cos (r, r0 ) = | cos (v, w)|

84 CAPÍTULO 3. GEOMETRÍA

Ángulo de dos planos

Dos planos π y π 0 que se cortan forman cuatro diedros que son iguales dos a dos. Como
en el caso de dos rectas que se cortan en el plano los diedros distintos son suplementarios
(su suma es de 180◦ ) y tomaremos el menor de éstos para medir el ángulo de los planos.

Si las ecuaciones de los planos son

π : ax + by + cz + d = 0
π 0 : a0 x + b0 y + c0 z + d0 = 0

Los vectores n y n0 de coordenadas (a, b, c) y (a0 , b0 , c0 ) son perpendiculares a los planos

π y π 0 . La recta de intersección de los dos planos es perpendicular a los vectores n y n0
y como el ángulo que se mide entre los dos planos se hace en un plano perpendicular a
ambos, es decir un plano π 00 perpendicular a la recta de intersección. El plano π 00 define
dos rectas π ∩ π 00 y π 0 ∩ π 00 y medir el ángulo de estas dos rectas es lo mismo que medir el
ángulo entre π y π 0 . Los vectores n y n0 son ortogonales a estas rectas y nos sirven para
hallar el ángulo buscado.

|aa0 + bb0 + cc0 |

cos (π, π 0 ) = | cos (n, n0 )| = √ √
a2 + b2 + c2 a02 + b02 + c02

Ángulo de recta y plano

Supongamos que la recta r y el plano π se cortan en un punto P . Sea r0 la proyección
ortogonal de la recta r sobre el plano π. Esta proyección la podemos obtener intersecando
el plano π con el plano π 0 que es perpendicular a π y que contiene a la recta r. Entonces
el ángulo buscado α es el formado por las rectas r y r0 .

Los vectores n = (a, b, c) normal al plano π y v de dirección de la recta r son paralelos al

plano π 0 . El ángulo formado por los vectores n y v es el complementario del formado por
las rectas r y r0 .

Como el coseno de un ángulo es igual al coseno de su complementario tenemos que

sen (α) = | cos (n, v)|

y en coordenadas
|av1 + cv2 + dv3 |
sen (r, π) = √ p
a + b2 + c2 v12 + v22 + v32
2
3.6. EL ESPACIO EUCLIDIANO 85

Distancia de un punto a un plano

Sea π un plano de ecuación ax + by + cz + d = 0 y P un punto de coordenadas (p1 , p2 , p3 )
con respecto a un sistema de referencia métrico. La perpendicular de P a π corta a este
plano en un punto M y lo que queremos hallar es la longitud del segmento M P . Si las
coordenadas del punto M , que no conocemos, son (m1 , m2 , m3 ), lo que queremos hallar
es p
d(P, π) = |M P | = (p1 − m1 )2 + (p2 − m2 )2 + (p3 − m3 )2

Del punto M sabemos que pertenece al plano π y que el vector M P es perpendicular a

este plano, y por tanto paralelo al vector (a, b, c), entonces

p1 − m1 p2 − m2 p3 − m3
= = =k
a b c
donde k es la razón de proporcionalidad. De aquı́

p1 − m1 = ka, p2 − m2 = kb, p3 − m3 = kc, (3.9)

y sumando los cuadrados de los miembros de la izquierda de 3.9 tenemos

(p1 − m1 )2 + (p2 − m2 )2 + (p3 − m3 )2 = k 2 (a2 + b2 + c2 ) (3.10)

El punto M pertenece al plano π y entonces

am1 + bm2 + cm3 = 0

despejando los valores de m1 , m2 , m3 en 3.9 y sustituyendo en la ecuación anterior se tiene

a(p1 − ka) + b(p2 − kb) + c(p3 − kc) + d = 0

ap1 + bp2 + cp3 + d − k(a2 + b2 + c2 ) = 0

y de aquı́
ap1 + bp2 + cp3 + d
k=
a2 + b 2 + c 2
por otra parte de 3.10 se tiene
p √
(p1 − m1 )2 + (p2 − m2 )2 + (p3 − m3 )2 = ±k a2 + b2 + c2

donde la anbigüedad del signo del segundo miembro depende de si k es mayor o menor
que 0. Vemos que el primer miembro es la distancia que estábamos buscando. Entonces
si sustituimos en el segundo el valor de k calculado antes, tenemos

|ap1 + bp2 + cp3 + d|

d(P, π) = √
a2 + b 2 + c 2
86 CAPÍTULO 3. GEOMETRÍA

Área de un triángulo
Sea el triángulo determinado por los puntos A = (a1 , a2 , a3 ), B = (b1 , b2 , b3 ) y C =
(c1 , c2 , c3 ). De la definición de producto vectorial de dos vectores se tiene

1
Area(ABC) = AB × AC
2

Distancia de un punto a una recta

Supongamos la recta dada por las ecuaciones
x − a1 y − a2 z − a3
= =
v1 v2 v3
y consideremos un punto P = (p1 , p2 , p3 ) que no pertenece a la recta.

De la ecuación sabemos que la recta pasa por el punto A = (a1 , a2 , a3 ). Consideremos el

representante AV del vector v de dirección de la recta, el punto V estará en la recta. Y
sea P H el segmento perpendicular desde P hasta la recta. Estamos buscando la longitud
de P H.

Consideremos el triángulo P AV . El área de este triángulo será

1
Area(AP V ) = AP × v
2
que en coordenadas resulta
s
2 2 2
1 p2 − a2 v2 p − a1 v1 p − a1 v2
Area(AP V ) = + 1 + 1
2 p3 − a3 v3 p3 − a3 v3 p2 − a2 v2

Por otra parte el área del triángulo AP V es igual a

1
Area(AP V ) = AV · P H
2
y q
AV = |v| = v12 + v22 + v32
igualando las áreas se tiene
1 1
AP × v = |v| · P H
2 2
de donde
AP × v
d(P, r) = P H =
|v|
que en coordenadas es
s
2 2 2
p 2 − a2 v 2 p1 − a1 v1 p1 − a1 v2
+ +
p 3 − a3 v 3 p3 − a3 v3 p2 − a2 v2
d(P, r) =
v12 + v22 + v32
3.6. EL ESPACIO EUCLIDIANO 87

Distancia entre dos rectas que se cruzan

Consideremos las rectas dadas por las ecuaciones

r : x = a + tv, r : x = b + sw

La recta r pasa por el punto A cuyo vector de posición es a = (a1 , a2 , a3 ) y tiene dirección
del vector v = (v1 , v2 , v3 ). La recta r0 pasa por el punto B cuyo vector de posición es
b = (b1 , b2 , b3 ) y tiene dirección del vector w = (w1 , w2 , w3 )

Sea π el plano que contiene a la recta r y es paralelo a r0 , este plano pasa por el punto A
y es paralelo a los vectores v y w. Entonces la ecuación de π será

v1 w1 x − a1
v2 w2 y − a2 =0
v3 w3 z − a3

y desarrollando se tienen que los coeficientes de las variables son los menores

v2 w2 v w1 v1 w1
, − 1 ,
v3 w3 v3 w3 v2 w2

Todos los puntos de la recta r0 distan lo mismo del plano π y esa es la distancia mı́nima
entre las rectas r y r0 . Entonces esta distancia es la distancia de un punto cualquiera de
la recta r0 a π. Si elegimos como punto B tenemos

v1 w1 b1 − a1
v2 w2 b2 − a2
v3 w3 b3 − a3
d(r, r0 ) = s
2 2 2
v2 w2 v w1 v w1
+ 1 + 1
v3 w3 v3 w3 v2 w2

Volumen de un tetraedro
Ver ejercicio en la sección siguiente.
88 CAPÍTULO 3. GEOMETRÍA

3.7. Ejercicios

1. Sean u y v dos vectores del espacio R3 que forman entre sı́ un ángulo de π/4. Probar
que u · v = ku × vk.

2. Sean u, v y w tres vectores del espacio R3 tales que kuk = 1, kvk = 2, kwk = 3 y
el ángulo entre dos cualesquiera de ellos es π/3. Calcular ku + v + wk.
Solución: Se tiene que cos (π/3) = 1/2. Entonces

ku + v + wk2 = (u + v + w) · (u + v + w)
= (u, u) + (v, v) + (w, w) + 2(u, v) + 2(u, w) + 2(v, w)
π π π
= kuk2 + kvk2 + kwk2 + 2kukkvk cos + 2kukkvk cos 2kukkvk cos
3 3 3
= 14 + 4 + 3 + 4
= 25
√
y entonces ku + v + wk = 25 = 5.

3. Encontrar la ecuación del plano que:

a) Es perpendicular a v = (1, 1, 1) y pasa por (1, 0, 0).

b) Es perpendicular a la recta l(t) = (5, 0, 2)t + (3, −1, 1) y pasa por el punto
(5, −1, 0).

Solución:

a) El vector normal al plano es n = i + j + k y la ecuación de éste será

(x − 1) + (y − 0) + (z − 0) = 0

es decir
x+y+z =1

b) El vector de dirección de la recta tiene que ser perpendicular al plano, es decir

podemos tomar como vector normal n = (5, 0, 2), y la ecuación del plano se
obtiene finalmente al imponer que el punto (5, −1, 0) pertenece a éste:

5 · (x − 5) + 0 · (y + 1) + 2 · (z − 0) = 0

o lo que es lo mismo
5x + 2z = 25

4. Hallar la intersección de los planos x + 2y + z = 0 y x − 3y − z = 0.

Solución: Se trata de encontrar los puntos que satisfacen el siguiente sistema

x + 2y + z = 0
x − 3y − z = 0
3.7. EJERCICIOS 89

que es equivalente al sistema


x + 2y = −z
x − 3y = z

De la segunda ecuación tenemos

x = z + 3y (3.11)

y sustituyendo en la primera
2
z + 5y = −z ⇒ y = − z
5
y sustituyendo este valor en (1) tenemos
6 1
x=z− z=− z
5 5
esto nos da la siguiente ecuación paramétrica de la intersección
 
1 2
(x, y, z) = − λ, − λ, λ
5 5
que es la de una recta que pasa por el origen, ya que los dos planos pasan por el
origen.
El sistema del inicio lo podrı́amos haber resuelto del siguiente modo: sumando las
dos ecuaciones y restando la segunda a la primera, tenemos respectivamente

2x − y = 0 y 5y = −2z

e igualando se tiene
5 x y z
2x = y = − z ≡ = =
2 1/2 1 −5/2
que es la ecuación
x y z
= =
−5 −2/5 1
que equivale la paramétrica calculada antes.
Tambien podrı́amos haber resuelto el sistema utilizando la regla de Cramer

−z 2 1 −z
z −3 1 z
x= , y=
1 2 1 2
1 −3 1 −3
.
5. Sean u, v ∈ R3 dos vectores con normas 4 y 2 respectivamente. Si ku × vk = 5,
calcular u · v.
90 CAPÍTULO 3. GEOMETRÍA

6. Calcular la distancia desde el origen de R3 al plano de ecuación x + 2y + 3z = 4.

7. Encontrar la ecuación del plano que:

a) Es perpendicular a v = (1, 1, 1) y pasa por (1, 0, 0).

b) Es perpendicular a la recta l(t) = (5, 0, 2)t + (3, −1, 1).

Solución:

a) El vector normal al plano es n = i + j + k y la ecuación de éste será

(x − 1) + (y − 0) + (z − 0) = 0

es decir
x+y+z =1

b) El vector de dirección de la recta tiene que ser perpendicular al plano, es decir

podemos tomar como vector normal n = (5, 0, 2), y la ecuación del plano se
obtiene finalmente al imponer que el punto (5, −1, 0) pertenece a éste:

5 · (x − 5) + 0 · (y + 1) + 2 · (z − 0) = 0

o lo que es lo mismo
5x + 2z = 25

8. Hallar la intersección de los planos x + 2y + z = 0 y x − 3y − z = 0.

Solución: Se trata de encontrar los puntos que satisfacen el siguiente sistema

x + 2y + z = 0
x − 3y − z = 0

que es equivalente al sistema


x + 2y = −z
x − 3y = z

De la segunda ecuación tenemos

x = z + 3y (3.12)

y sustituyendo en la primera
2
z + 5y = −z ⇒ y = − z
5
y sustituyendo este valor en (1) tenemos

6 1
x=z− z=− z
5 5
3.7. EJERCICIOS 91

esto nos da la siguiente ecuación paramétrica de la intersección

 
1 2
(x, y, z) = − λ, − λ, λ
5 5

que es la de una recta que pasa por el origen, ya que los dos planos pasan por el
origen.
El sistema del inicio lo podrı́amos haber resuelto del siguiente modo: sumando las
dos ecuaciones y restando la segunda a la primera, tenemos respectivamente

2x − y = 0 y 5y = −2z

e igualando se tiene
5 x y z
2x = y = − z ≡ = =
2 1/2 1 −5/2

que es la ecuación
x y z
= =
−5 −2/5 1
que equivale la paramétrica calculada antes.
Tambien podrı́amos haber resuelto el sistema utilizando la regla de Cramer

−z 2 1 −z
z −3 1 z
x= , y=
1 2 1 2
1 −3 1 −3

9. Obtener un vector unitario con componente k positiva que sea perpendicular a los
vectores u = 2i − j + 2k y v = 2i − 3j − k.
Solución: Un vector perpendicular a los vectores u y v es el producto vectorial
u × v:

i j k
u×v = 2 −1 2 = 7i + 6j − 4k
2 −3 −1
√
Tenemos que |u × v| = 101, pero el vector u × v tiene la componente k negativa,
por lo que el vector pedido es
1 1
−√ |u × v| = √ (−7i − 6j + 4k)
101 101

10. Verificar las siguientes igualdades

a) u × u = 0.
92 CAPÍTULO 3. GEOMETRÍA

b) u × v = −v × u.
c) (u + v) × w = u × w + v × w.
d) (tu) × v = u × (tv) = t(u × v)
e) u · (u × v) = v · (u × v) = 0.

Solución:

a) El vector u forma un ángulo 0 consigo mismo. Por ello tenemos |u × u| =

|u| · |u| · sen 0 = 0 y entonces u × u = 0. Tambien se puede ver que

i j k
u×u= u1 u2 u3 =0
u1 u2 u3
ya que en el determinante tenemos dos filas iguales.
b) Tenemos

i j k i j k
u×v = u1 u2 u3 = − v1 v2 v3 = −v × u
v1 v2 v3 u1 u2 u3

teniendo en cuenta que permutamos dos filas.

c) Teniendo en cuenta de nuevo las propiedades de los determinantes

i j k
(u + v) × w = u1 + v1 u2 + v2 u3 + v3 =
w1 w2 w3
i j k i j k
= u1 u2 u3 + v1 v2 v3 =u×w+v×w
w1 w2 w3 w1 w2 w3

d)

i j k i j k i j k
(tu) × v = tu1 tu2 tu3 = t u1 u2 u3 = u1 u2 u3 =
v1 v2 v3 v1 v2 v3 tv1 tv2 tv3
= t(u × v) = u × (tv)

e)

u · (u × v) =
u2 u3 u1 u3 u1 u2
= u1 − u2 + u3 =
v2 v3 v1 v3 v1 v2
= u1 u2 u3 − u1 v2 u3 − u2 u1 v3 + u2 v1 u3 + u3 u1 v2 − u3 v1 u2 = 0.

v · (u × v) = −v · (v × u) = 0.
3.7. EJERCICIOS 93

11. El volumen de un tetraedro es 31 A · h, donde A es el área de la base y h la altura

medida perpendicularmente desde el vértice opuesto a la base. Sean u, v y w tres
vectores que coinciden con las aristas del tetraedro, que concurren en un vértice.
Demostrar que el volumen del tetraedro puede expresarse como
u1 u2 u3
1 1
V = |u · (v × w)| = | v1 v2 v3 |
6 6
w1 w2 w3

Solución:
La base del tetraedro es el triángulo que definen los dos vectores v y w, cuya área
es
1
A = |v × w|
2

v×w

h w

La altura h del tetraedro es igual a la longitud de la proyección de u sobre el vector

v × w (este vector es perpendicular a la base). Tenemos
|u · (v × w)|
h=
v×w
Y tenemos que el volumen del tetraedro es
u1 u2 u3
1 1 1
V = Ah = |u · (v × w)| = | v1 v2 v3 |.
3 6 6
w1 w2 w3

12. Demostrar que el plano que pasa por los tres puntos A = (a1 , a2 , a3 ), B = (b1 , b2 , b3 )
y C = (c1 , c2 , c3 ) está formado por los puntos P = (x, y, z) tales que
a1 − x a 2 − y a 3 − z
b1 − x b2 − y b3 − z = 0.
c1 − x c2 − y c3 − z
94 CAPÍTULO 3. GEOMETRÍA

−→ −→
Solución: El vector AB × AC es ortogonal al plano determinado por A, B y C.
−→
Este vector define el plano como el conjunto de puntos P = (x, y, z) tal que P C
−→ −→ −→ −→ −→
es un vector combinación lineal de AB y AC y por tanto (AB × AC) · P C = 0.
Expresando esta última ecuación tenemos

b1 − a1 b2 − a2 b3 − a3
c 1 − a1 c 2 − a2 c 3 − a3 =0
c1 − x c2 − y c3 − z

y restando a la segunda fila, la tercera y cambiándola de signo se tiene

b1 − a1 b2 − a2 b3 − a3
a1 − x a2 − y a3 − z =0
c1 − x c2 − y c3 − z

ahora procediendo de manera análoga con la primera fila y reordenando después,

tenemos la condición del enunciado.
Capı́tulo 4

Análisis

4.1. Funciones elementales

Funciones afines
Son las de la forma f (x) = mx + n, cuando n = 0 se llaman lineales y cuando m = 0 se
llaman constantes

Funciones potenciales
Son las de la forma f (x) = ±xn . Su representación gráfica es sencilla.

Funciones polinómicas
Son las que tienen la forma f (x) = an xn + an−1 xn−1 + · · · + a1 x + a0 y las anteriores son
un caso particular de éstas.

La función de Dirichlet
Es la función


0 si; x ∈ Q
f (x) =
1 si; x ∈ R − Q

La función valor absoluto

Es la función


x si; x > 0
f (x) =
−x si; x < 0

95
96 CAPÍTULO 4. ANÁLISIS

Funciones trigonométricas
Si consideramos la circunferencia C : x2 + y 2 = 1 en el plano, con centro el origen de
coordenadas y radio 1, un punto (a, b) de C determina un ángulo rectángulo con el eje X
que podemos medir en radianes x. Entonces tenemos

b
sen x = b cos x = a tg x =
a

Relaciones de las funciones trigonométricas

La relación fundamental es
sen2 x + cos2 x = 1
Además también se tiene

1
1 + tg2 x =
cos2 x
Las de adición

sen (x ± y) = sen x cos y ± cos xsen y

cos (x ± y) = cos x cos y ∓ sen xsen y
tg x ± tg y
tg (x ± y) =
1 ∓ tg xtg y
4.1. FUNCIONES ELEMENTALES 97

Funciones exponenciales
La función exponencial general es
f (x) = ax

donde a es un número real mayor que cero.

En la siguiente figura aparecen las gráficas de ex , 2x , 4x , e−x , 2−x y 4−x . La primera de

ellas está dibujada en color rojo. A qué colores corresponden las otras cinco?

Figura 4.1: Funciones exponenciales. Las crecientes (que se aproximan al eje X negativo)
tienen base positiva. Las decrecientes tienen base negativa. En rojo la función f (x) = ex .
Obsérvese que la escala está distorsionada.

La función exponencial que tiene un interés especial es la que tiene como base el número
e; esta es f (x) = ex .

Definición El número e es:

 n
1
e = lı́m 1 +
n n

Nota En la siguiente sección se trata de manera breve el concepto de sucesión.

98 CAPÍTULO 4. ANÁLISIS

Figura 4.2: El dibujo muestra las gráficas de las funciones anteriores con su escala real, lo
que da una idea de la rapidez de su crecimiento.

Funciones logarı́tmicas
La función logaritmo en base a (a > 0), digamos loga x es la inversa de ax . Nosotros
utilizaremos, salvo que se indique lo contrario, la función log x, que es la función logaritmo
en base e y que a veces se escribe ln x.

Propiedades de las funciones logarı́tmicas

Las propiedades, que se cumplen para cualquier base, son:

x
log (xy) = log x + log y log = log x − log y log xn = n log x
y

Nota la función logaritmo es la inversa de la eponencial, ası́ sus gráficas son simétricas
respecto a la recta y = x.
4.2. SUCESIONES 99

4.2. Sucesiones
Definición Se llama sucesión de elementos de un cierto conjunto X a una aplicación

s : N −→ X

Los elementos de una sucesión se llaman términos y si el conjunto X es el conjunto R de

los números reales, tendremos una sucesón de números reales.

Los términos de una sucesión se representan con subı́ndices a1 , a2 , a2 , · · · y el término

que ocupa un lugar cualquiera se denota an y se llama término general. Las sucesiones se
representan por su término general: (an ) o simplemente an .

Las sucesiones de más fácil manejo son aquellas que su término general puede expresarse
por medio de una fórmula, por ejemplo:

1. an = 3n − 7 → a1 = −4, a2 = −1, a3 = 2, a4 = 5, · · ·
n2 − 2n + 1 4 9
2. → b 1 = 0, b2 = 1, b3 = , b4 = ,···
n2 − 3 6 13
 n
1 9 64 625
3. xn = 1 + → x1 = 2, x2 = , x3 = , x4 = ,···
n 4 27 256

Lı́mite de una sucesión

Lo que interesa de una sucesión normalmente es su tendencia, es decir, su comporta-
miento en el infinito. Esto nos proporciona una herramienta muy interesante ya que la
propia sucesión determina un valor cuya obtención puede requerir un proceso infinito. El
ejemplo clásico es el de la construcción de los números reales como lı́mites de sucesiones
de números racionales. Más concretamente, en el punto 3 anterior tenemos la sucesión xn
que representa al número e, que juega un papel muy destacado en Matemáticas.

La relación entre la sucsión xn anterior y el número e es que a medida que los términos de
la sucesión crecen, ésta se va acercando cada vez más a este número, que por otra parte
se puede demostrar que no es un número racional.

Necesitamos precisar el concepto de lı́mite de una sucesión.

El lı́mite de la sucesión an es un número l. Podemos obtener un valor de an

tan próximo a l como queramos, eligiendo un término an0 suficientemente alto de la
sucesión, y además todos los términos de ı́ndice mayor que n0 estarán por lo menos
tan cerca como an0 de l.
lı́m an = l
n

Dicho con precisión lo anterior serı́a: ∀ε > 0, existe un n0 tal que

si n > n0 , entonces l − ε < an < l + ε, o bien |an − l| < ε.

100 CAPÍTULO 4. ANÁLISIS

El lı́mite de la sucesión an es ∞. Podemos conseguir que an sea tan grande como

queramos tomando n suficientmente grande
lı́m an = ∞
n

Dicho en lenguaje matemático serı́a: ∀K, existe n0 tal que

si n > n0 , entonces an > K

El lı́mite de la sucesión an es −∞. Este caso es análogo al anterior, basta

considerar K negativo.

El número e.
Consideremos la sucesión:

 n
1
xn = 1 +
n

Los valores de esta sucesión para n grande parecen aproximarse a un valor, aunque como se
puede ver muy lentamente. A continuación escribimos estos valores que tienen un número
finito de cifras decimales, pero no las escribimos todas.

 10
1
x10 = 1 + = 10 110 = 20 59374...
10
 100
1
x100 = 1 + = 10 1100 = 20 70481...
100
 1000
1
x1000 = 1 + = 10 11000 = 20 71692...
1000
 1000000
1
x1000000 = 1 + = 10 110 = 20 71828047...
1000000
 1000000000
1
x1000000000 = 1 + = 10 11000000000 = 20 728181827...
1000000000

Su lı́mite es un número con infinitas cifras decimales en las que no existe ningún perı́odo
por lo tanto es un número irracional :
 n
1
e = lı́m 1 + = 20 7182818284...
n n
También se puede expresar e como una suma infinita (la correspondiente a evaluar el
desarrollo de Taylor de la función f (x) = ex en un entorno de 0.)
1 1 1 1 1
e=1+ + + + + ··· + ···
1! 2! 3! 4! n!
4.3. SENO Y COSENO DE UN ÁNGULO 101

4.3. Seno y coseno de un ángulo

El seno y el coseno son medidas angulares aunque se puede pasar a una representación
en x ∈ R si se expresan los ángulos en radianes.

Si tomamos en el plano R2 con coordenadas (x, y), la circunferencia de radio 1 centrada

en el origen, cada punto P ≡ (a, b) de esta circunferencia determina un ángulo: el que
forma el eje X con el radio que define P , medido por este orden. Entonces si este ángulo
es θ tenemos
sen θ = b y cos θ = a

Dicho de manera geométrica, si el punto P está en el cuadrante (+, +) el radio OP y los

segmentos P A y OA forman un triángulo rectángulo, donde A ≡ (a, 0) y O es el origen
de coordenadas.

Coordenadas polares
Vista de la definición anterior de seno y coseno, podemos determinar la posición de un
punto P ≡ (x, y) cualquiera del plano, salvo (0, , 0) por sus coordenadas polares ρ y θ,
donde:

x = ρ cos θ
y = ρ sen θ
102 CAPÍTULO 4. ANÁLISIS

4.4. El concepto de lı́mite de una función en un

punto
Este concepto es fundamental en esta área de las Matemáticas. Empezaremos por el lı́mite
de una función en un punto, pero antes de dar su definición formal explicaremos la idea
intuitiva, idea que incluso nos servirá para el cálculo efectivo.

Si una función f (x) se acerca al valor l, a medida que x se acerca al valor x0 decimos que
el lı́mite de f cuando x tiende a x0 es l y se escribe :

lı́m f (x) = l
x→x0

Observaciones

1. El lı́mite de f cuando x → x0 no tiene nada que ver con el valor f (x0 ), estos pueden
coincidir o no, incluso pueden existir independientemente.

2. Para determinar el lı́mite de f cuando x → x0 tenemos que analizar la función f a

ambos lados de x0 .

3. Puede darse el caso de que mientras que x se acerca a x0 , los valores f (x) no se
aproximan a ningún valor fijo, con lo que el lı́mite de f cuando x → x0 no existe.

Propiedades básicas de los lı́mites

Supongamos que existen lı́m f (x) y lı́m g(x), entonces:

x→x0 x→x0

1. lı́m [f (x) + g(x)] = lı́m f (x) + lı́m g(x)

x→x0 x→x0 x→x0
2. lı́m [f (x) · g(x)] = lı́m f (x) · lı́m g(x)
x→x0   x→x0 x→x0
1 1
3. lı́m = si lı́m f (x) 6= 0
x→x0 f (x) lı́m f (x) x→x0
x→x0

Observaciones:

i) El lı́mite de la función constante es el propio valor de la función, o sea lı́mx→x0 c = c.

ii) Si dos funciones f y g toman los mismos valores en un intervalo abierto (a, b) que
contiene a x0 , aunque las dos funciones tomen distintos valores en x0 , se tiene:
4.4. EL CONCEPTO DE LÍMITE DE UNA FUNCIÓN EN UN PUNTO 103

lı́m f (x) = lı́m g(x)

x→x0 x→x0

Con la noción de lı́mite que tenemos, aun sin haber dado la definición precisa, y con las
propiedades anteriores, ya podemos hacer algunos cálculos.

|x|
Ejemplo 16 ¿Existe lı́m ?
x→0 x
Solución: La función puede escribirse


x/x si x > 0
f (x) =
−x/x si x < 0

La función f , cuya gráfica aparece a continuación, tiende a 1 cuando x tiende a 0 pero

es mayor que 0 y tiende a −1 cuando x tiende a 0 pero es menor que 0. Ası́ resulta que
lı́mx→0 f (x) no existe.
104 CAPÍTULO 4. ANÁLISIS

Ejemplo 17 Calcular el valor de:

x3 − 2x2 − 5x + 6
i) lı́m
x→0 x2 − 3x + 2
x3 − 2x2 − 5x + 6
ii) lı́m
x→1 x2 − 3x + 2

Solución: En el primer caso calculamos los lı́mites del numerador y denominador apli-
cando las propiedades expuestas antes y tenemos:

lı́m (x3 − 2x2 − 5x + 6) = lı́m x3 − lı́m 2x2 − lı́m 5x + lı́m 6

x→0 x→0 x→0 x→0 x→0
 3  2
= lı́m x − 2 lı́m x − 5 lı́m x + lı́m 6
x→0 x→0 x→0 x→0
= 6

ya que lı́mx→0 x = 0. Análogamente se tiene para el denominador

lı́m (x2 − 3x + 2) = 2
x→0

y como consecuencia tenmos que

x3 − 2x2 − 5x + 6
lı́m =3
x→0 x2 − 3x + 2

En el segundo caso, calculando los lı́mites del numerador y denominador, se comprueba

que ambos son igual a 0. Como se trata de polinomios, deducimos que x = 1 es raı́z, tanto
del numerador como del denominador, entonces podemos poner:

x3 − 2x2 − 5x + 6 (x2 − x − 6)(x − 1) x2 − x − 6

= =
x2 − 3x + 2 (x − 2)(x − 1) x−2
donde en esta última expresión el denominador no toma el valor 0 para x = 1. Lo que
hemos obtenido es la expresión de un cociente que toma los mismos valores que el pri-
mero en un intervalo (a, b) que contiene a x = 1, por ejemplo (00 8, 10 2), salvo para el
punto x = 1, donde precisamente el primer denominador no está definido (a/0, en este
caso concretamente (0/0)). Entonces con el denominador simplificado podemos calcular
el lı́mite ya que no hay que dividir por cero y tenemos
4.5. CONTINUIDAD DE UNA FUNCIÓN 105

x3 − 2x2 − 5x + 6 x2 − x − 6 −6
lı́m 2
= lı́m = =6
x→1 x − 3x + 2 x→1 x−2 −1

4.5. Continuidad de una función

Este concepto va asociado al de lı́mite y se define puntualmente, aunque luego se puede
extender a un conjunto.

Definición Se dice que una función f es continua en x = x0 si

lı́m = f (x0 )
x→x0

Observación De la definición se deduce que f (x0 ) tiene que estar definido, de lo contrario,
f no serı́a continua en ese punto.

4.6. Lı́mites infinitos y lı́mites en el infinito

Consideremos la función f (x) = 1/x, cuya gráfica aparece representada en la figura 4.6.
Observamos que cuando x se aproxima a 0 y es mayor que éste, el valor f (x) = 1/x se
hace cada vez más grande y cuando se aproxima a este valor pero siendo menor que él, el
valor de la imagen se hace cada vez más grande pero en sentido negativo. Esta situación
cuya definición se precisa como la de lı́mite, se conoce como lı́mites infinitos, y en este
caso concreto escribirı́amos:

lı́m = +∞ y lı́m = −∞
x→0+ x→0−

Por otra parte, si nos fijamos en la gráfica de la figura 4.6, vemos que a medida que x
crece en sentido positivo, su imagen por f se acerca a 0. Lo mismo si x crece en sentido
negativo. Esta situación se puede precisar al igual que la anterior y se conoce como lı́mites
en el infinito. En este ejemplo tendrı́amos:

lı́m = 0 y lı́m = 0
x→+∞ x→−∞
106 CAPÍTULO 4. ANÁLISIS

Observación: Las propiedades de los lı́mites cuando x → ±∞ son las mismas que en el
caso x → x0 .

x2 − 3x + 3
Ejemplo 18 Calcular lı́m
x→+∞ 3x2 + 4x − 5

Solución: Si dividimos el numerador y el denominador por x2 tenemos:

x2 − 3x + 3 1 − 3/x + 3/x2 1
lı́m = lı́m =
x→+∞ 3x2 + 4x − 5 x→+∞ 3 + 4/x − 5/x2 3

1 1
ya que lı́mx→+∞ x
= 0 y también, como consecuencia de lo anterior, lı́mx→+∞ x2
= 0.
4.7. LA DERIVADA 107

4.7. La derivada
La derivada de una función f en un punto a es la inclinación de la recta tangente a la
gráfica de f en el punto (a, f (a)).

Ejemplo Consideremos la función f (x) = x2 − 1, y supongamos que queremos calcular

la derivada en el punto a = 1. Si consideramos una recta que pase por el punto (a, f (a))
y que corta a la gráfica de f en otro punto arbitrario, digamos (x, f (x)), el cociente
f (x) − f (a)
x−a
nos da la inclinación de esta recta secante. Para llegar a la situación que nos interesa,
pasamos al lı́mite cuando x → a y tenemos la derivada en es punto:

f (x) − f (a)
f 0 (a) = lı́m
x→a x−a

En nuestro ejemplo tendremos:

x2 − 1 − 0 (x + 1)(x − 1)
f 0 (a) = lı́m = lı́m = lı́m (x + 1) = 2
x→1 x−1 x→1 x−1 x→1

La figura anterior muestra el proceso de paso de la secante a la tangente en el ejemplo

tratado.
108 CAPÍTULO 4. ANÁLISIS

Derivada de una función en un punto y función derivada

Acabamos de ver que la derivada es un concepto puntual. En general suele escribirse la
derivada de la función f en el punto x = a como:

f (a + h) − f (a)
f 0 (a) = lı́m
h→0 h

donde hemos tomado a + h = x.

Si una función f está definida en un conjunto A ⊂ R y existe la derivada de esta función

en los puntos de B ⊂ A, a la función que a cada b ∈ B asigna el valor f 0 (b) le llamamos
función derivada, que en este caso estará definida en B.

Derivadas laterales
Como ya hemos visto, la derivada es un lı́mite, y la existencia de este implica la de los
lı́mites laterales, veamos un ejemplo sencillo:

Ejemplo 19 Calcular la derivada de la función f (x) = |x|.

Solución: La función f puede escribirse de la siguiente forma:


x si x > 0
f (x) =
−x si x < 0

entonces podemos calcular f 0 (x) en los intervalos (−∞, 0) y (0, +∞) y tenemos:


0 1 si x > 0
f (x) =
−1 si x < 0

Ahora, para calcular la derivada en x = 0, lo hacemos aplicando la definición, y tenemos

que calcular el lı́mite: lı́mx→0 f (x)−f
x−0
(0)
. La expresión de la función, nos obliga a calcular
los lı́mites laterales:

f (x) − f (0) f (x)

lı́m+ = lı́m+ =1
x→0 x−0 x→0 x
4.7. LA DERIVADA 109

f (x) − f (0) f (−x)

lı́m− = lı́m− = −1
x→0 x−0 x→0 x

Y resulta entonces, que el lı́mite buscado no existe, por lo tanto f no es derivable en

x = 0.

Derivada de la función compuesta

La derivada de f se puede entender como una medida del crecimiento de la función f en
cada punto. Este crecimiento está referido a la variable x. Si esta variable dependiese a
su vez de otra, x(u), se nos plantea la situación de expresar la derivada de f en función
de la variable u.

Lo que tenemos es entonces una composición de funciones, y la regla de la cadena nos

dice cual es la derivada de esta composición con respecto a la variable independiente.

Regla de la cadena
Si f es una función derivable en g(x) y g lo es en x, la derivada de f ◦ g en x es:

(f ◦ g)0 (x) = f 0 (g(x)) · g 0 (x)

√
Ejemplo 4.1. Hallar la derivada de h(x) = x3 + 2.

√
Solución: Tenemos la composición h(x) = (f ◦ g)(x), donde f (x) = x y g(x) = x3 + 2.
Entonces, como f 0 (x) = 2√1 x y g 0 (x) = 3x2 , se tiene f 0 (g(x)) = 2√x13 +2 y la derivada
buscada es:

1
f 0 (x) = √ · 3x2
3
2 x +2

Máximos y mı́nimos locales

Sea f una función definida en un conjunto A ⊂ R. La función f tiene un extremo local,
máximos o mı́nimo, en un punto x = a si f (a) > f (x) para x ∈ (x − ε, x + ε) ∩ A, en cuyo
caso serı́a un máximo, o bien si f (a) 6 f (x) para x ∈ (x − ε, x + ε) ∩ A, en cuyo caso
serı́a un mı́nimo.
110 CAPÍTULO 4. ANÁLISIS

En el caso de que f tenga un extremo en x = a y f sea derivable en ese punto, entonces

f 0 (a) = 0. El recı́proco no es cierto.

4.8. Teoremas de Rolle, Cauchy y del Valor Medio

El teorema de Rolle refleja analı́ticamente un hecho geométrico evidente sobre la gráfica
de una función derivable, puede enunciarse de la siguiente manera.
Teorema 4.1 (Rolle). Sea f una función continua en [a, b] y derivable en (a, b). Entonces,
si f (a) = f (b), se puede asegurar que existe un punto ξ ∈ (a, b) de forma que f 0 (ξ) = 0.

Este teorema puede considerarse una versión más general del teorema de Rolle.
Teorema 4.2 (Valor Medio). Sea f una función continua en [a, b] y derivable en (a, b).
Entonces, existe ξ ∈ (a, b) de manera que
f (b) − f (a)
f 0 (ξ) =
b−a
4.8. TEOREMAS DE ROLLE, CAUCHY Y DEL VALOR MEDIO 111

2
En la figura se muestra la gráfica de la función f (x) = x x+3 , sobre el intervalo [1, 5], el
Teorema del Valor Medio, asegura que existe un punto ξ ∈ (1, 5) donde la tangente a la
gráfica tiene la misma inclinación que la recta que pasa por los puntos (1, f (1)) y (5, f (5)).

Ejercicio 4.1. Probar que si una función tiene derivada nula, entonces es constante.

on

Teorema 4.3 (de Cauchy). Sean f y g son funciones continuas en [a, b] y derivables en
(a, b). Entonces existe un punto ξ ∈ (a, b) tal que:

f 0 (c)[g(b) − g(a)] = g 0 (c)[f (b) − f (a)]

Teorema 4.4 (Regla de L’Hôpital). Sean f y g son funciones derivables y tenemos la

indeterminación
lı́m f (x) = 0 y lı́m f (x) = 0
x→a x→a

donde a representa un número real o incluso ±∞, y

f 0 (x)
lı́m =l
x→a g 0 (x)

entonces
f (x)
lı́m =l
x→a g(x)
112 CAPÍTULO 4. ANÁLISIS

Ejemplo

1
ln x h ∞ i L0 H 1 1
lı́m = = lı́m x lı́m = = 0.
x→+∞ x ∞ x→+∞ 1 x→+∞ x +∞

Ejemplo
√
x+2 h∞i 0
LH 1 x2 − 3
lı́m √ = = lı́m = lı́m =
x→+∞ x2 − 3 ∞ x→+∞ √2x x→+∞ x
2 x2 −3
r r
x2 − 3 1 − 3/x2
lı́m = lı́m =1
x→+∞ x2 x→+∞ 1

Ejemplo
 
0 lg x L0 H 1/x
lı́m (tan x · lg x) = = lı́m 1 = lı́m
x→0 0 x→0
tan x
x→0 −1/ sin2 x
sin2 x L0 H 2 sin x cos x
= lı́m = − lı́m
x→0 x x→0 1
= 0.

Ejemplo

(2 − x)ex − x − 2 −ex + (2 − x)ex − 1

 
0 L0 H
lı́m 3
= = lı́m
x→0 x 0 x→0 3x2
0 L0 H −ex − ex + (2 − x)ex
 
= lı́m =
x→0 0 6x
0 L0 H −e − e − ex + (2 − x)ex
x x
 
= lı́m =
x→0 0 6
1
= −
6

Ejemplo

3x h ∞ i L0 H 3x log 3 h ∞ i L0 H 3x (log 3)2

lı́m = = lı́m = = lı́m
x→∞ x3 ∞ x→∞ 3x2 ∞ x→∞ 6x
h∞i 0 x 3
LH 3 (log 3) h ∞ i
= = lı́m = =∞
∞ x→∞ 6 6
4.9. LA DERIVADA SEGUNDA. CONCAVIDAD Y CONVEXIDAD. 113

Consecuencias del Teorema del Valor Medio

Teorema 4.5. Sea f continua en [a, b] y derivable en (a, b).

Si f 0 (x) > 0 para todo x ∈ (a, b), f es estrictamente creciente en [a, b].

Si f 0 (x) < 0 para todo x ∈ (a, b), f es estrictamente decreciente en [a, b].

4.9. La derivada segunda. Concavidad y convexidad.

Ejemplo Sean f (x) = x2 y g(x) = x3 . Las gráficas de estas dos funciones son bien
conocidas, y sabemos que son funciones derivables en todo R. Se cumple que f 0 (0) =
g 0 (0) = 0 y, sin embargo, f tiene un mı́nimo en ese punto pero g no.

En el caso de que una función sea dos veces derivable en un punto x = a podemos utilizar
la derivada segunda en ese punto, siempre que esta no valga 0 para determinar el tipo de
extremo:

Si f 00 (a) > 0 entonces f tiene un mı́nimo en x = a.

Si f 00 (a) < 0 entonces f tiene un máximo en x = a.

En el caso de que f 00 (a) = 0, tendrı́amos que recurrir a la tercera derivada y si esta fuese
cero, a la siguiente hasta llegar a la primera que fuese distinta de 0. El criterio a seguir
en este caso se explicará en el apartado correspondiente al Teorema de Taylor.

4.10. Teorema de Taylor

En esta sección vamos a tratar de aproximar una función en el entorno de un punto x = a
por otra más sencilla: un polinomio.

Esta aproximación se puede estimar en función de la variable independiente x, en el

siguiente sentido: Si la función que queremos aproximar es f y el polinomio aproximador
es P , la diferencia

|f (a + h) − P (a + h)|
es el error cometido al sustituir la función por el polinomio. Lo primero que podemos pedir
a nuestro polinomio es que por lo menos coincida con el valor de la función en el punto
a, entonces si el error anterior lo vemos como una función E de h, que es la distancia de
un punto x = a + h al punto a tendrı́amos que
114 CAPÍTULO 4. ANÁLISIS

lı́m E(h) = 0
h→0

Pero además podemos exigir a nuestro polinomio no sólo esto, sino que tienda a 0 con
cierta ’velocidad’. La palabra velocidad va entre comillas porque este es un concepto
fı́sico que quiere reflejar la diferencia que existe en cuanto a la forma de tender a cero,
por ejemplo, de dos funciones como f (x) = x y g(x) = x2 . Ambas tienden a 0 cuando x
tiende a 0, pero g lo hace más ’rápido’. Esto quiere decir que

g(x)
lı́m =0
x→0 f (x)

Infinitésimo
Se llama infinitésimo a una función ϕ(x) en x = a si lı́mx→a ϕ(x) = 0.

Orden de un infinitésimo
Decimos que un infinitésimo ϕ(x) tiene orden r si

ϕ(x)
lı́m
x→0 xr
existe y es distinto de cero. Del mismo modo podemos comparar dos infinitésimos f (x) y
g(x) entre si, calculando

f (x)
lı́m
x→0 g(x)

Si además este lı́mite vale 1 se dice que los infinitésimos son equivalentes. Si el lı́mite vale
0 será f (x) de orden superior a g(x), o sea, que tiende a 0 con mayor rapidez. Si lı́mite es
±∞ entonces será g(x) de orden superior a f (x), o sea, que tiende a 0 con mayor rapidez.

Para medir la rapidez de nuestro error, definimos:

Orden de aproximación
Para dos funciones f y g definidas en un entorno de a se dice que tienen orden de apro-
ximación r si

f (a + h) − g(a + h)
lı́m
h→0 xr
es un número distinto de 0.
4.10. TEOREMA DE TAYLOR 115

Teorema de Taylor
Sea f con derivadas hasta el orden n + 1 en un intervalo (a − k, a + k) de a. Entonces
tenemos:
f 0 (a) f 00 (a) f (n) (a) f (n+1) (ξ)
f (x) = f (a) + (x − a) + (x − a)2 + · · · + (x − a)n + (x − a)n+1
1! 2! n! (n + 1)!

donde ξ ∈ (x, a) ó ξ ∈ (a, x) dependiendo de si x < a ó x > a.

Si en la expresión del teorema anterior se suprime el último sumando nos queda el poli-
nomio de Taylor Pn (x; a) de grado n para la función f en el punto x = a.

El polinomio deTaylor de grado n aproxima a la función f , como mı́nimo hasta el orden

n + 1 en x = a: si se tiene f (n+1) (a) 6= 0, entonces:

f (x) − Pn (x; a)
lı́m = L 6= 0
x→a (x − a)n+1

x3 −4x+1
En las siguientes figuras se muestra en rojo la gráfica de la función f (x) = x−2
y en
azul los polinomios de Taylor de grados 1, 2, 4 y 8 en el punto x = 0.

Figura 4.3: P1 (x; 0) = − 21 + 74 x

116 CAPÍTULO 4. ANÁLISIS

Figura 4.4: P2 (x; 0) = − 12 + 74 x + 87 x2

Figura 4.5: P4 (x; 0) = − 12 + 47 x + 78 x2 − 1 3

16
x − 1 4
32
x

Figura 4.6: P8 (x; 0) = − 12 + 74 x + 87 x2 − 1 3

16
x − 1 4
32
x − 1 5
64
x − 1
128
x6 − 1
256
x7 − 1
312
x8

Ejemplo Comprobemos que los infinitésimos x, sen x, ex − 1, tg x y log (1 + x) son equi-

valentes en x = 0.
4.10. TEOREMA DE TAYLOR 117

Para ver que los dos primeros lo son, calculamos:

sen x cos x
lı́m = lı́m =1
x→0 x x→0 1
para lo que hemos aplicado la regla de L’Hôpital. Del mismo modo tenemos
ex −1 ex
lı́mx→0 x
= lı́mx→0 1
=1
1
tgx cos2 x
lı́mx→0 x = lı́mx→0 1
=1
1
lı́mx→0 log(1+x)
x
= lı́mx→0 1+x
1
=1

Ejercicio Estudiar si los siguientes infinitésimos en x = 0 son equivalentes: x2 , sen2 x y

1 − cos x.

Solución: Comparamos los dos primeros, y tenemos:

sen2 x L0 H 2 sen x cos x L0 H 2 cos2 x − 2 sen2 x

lı́m = lı́m = lı́m =1
x→0 x2 x→0 2x x→0 2
donde hemos aplicado dos veces la regla de L’Hôpital y entonces resulta que estos dos
infinitésimos son equivalentes.

Por otra parte

1 − cos x L0 H sen x L0 H cos x 1

lı́m 2
= lı́m = t lı́m =
x→0 x x→0 2x x→0 2 2
y entonces estos dos infinitésimos son del mismo orden pero no son equivalentes.

Ejemplo Si f (x) y g(x) son infinitésimos equivalentes a (x−a)r y (x−a)s respectivamente

cuando x → a, entonces:

f (x) f (x) (x − a)s (x − a)r (x − a)r

lı́m = lı́m = lı́m
x→a g(x) x→a (x − a)r g(x) (x − a)s x→a (x − a)s
y lo que resulta es que sustituimos las funciones por sus infinitésimos correspondientes en
un entorno del punto a. Este lı́mite es, en general, más fácil de calcular.

Importante Hay que tener en cuenta que la equivalencia de infinitésimos es local, es

decir, que tiene lugar en el entorno de un punto a, y entonces la sustitución que se hace
para calcular un lı́mite cuando x → a, como en el ejemplo anterior, sólo se puede hacer
cuando la equivalencia de los infinitésimos se tiene en el punto a.
118 CAPÍTULO 4. ANÁLISIS

Ejercicio Calcular el lı́mite

2x2 + 1 tg2 x + 1
 
lı́m −
x→0 sen2 x x2

Solución: Tenemos

2x2 + 1 tg2 x + 1 x2 (2x2 + 1) − sen2 x(tg2 x + 1)

 
lı́m − = lı́m
x→0 sen2 x x2 x→0 x2 sen2 x
y como x → 0, sustituimos en primer lugar sen2 x por el infinitésimo equivalente en 0,
que es x2 , posteriormente se sustituye el de tg2 , resulta

x2 (2x2 + 1) − sen2 x(tg2 x + 1) x2 (2x2 + 1) − sen2 x(tg2 x + 1)

lı́m = lı́m
x→0 x2 sen2 x x→0 x4
x (2x + 1) − sen2 x tg2 x + sen2 x
2 2
= lı́m
x→0
 4 x4
2x + x − sen2 x sen2 x tg2 x
2

= lı́m −
x→0 x4 x4
 4
2x + x2 − sen2 x sen2 x tg2 x

= lı́m − lı́m
x→0 x4 x→0 x4
 4
2x + x2 − sen2 x x2 x2

= lı́m − lı́m
x→0 x4 x→0 x4
 4 2 2

2x + x − sen x
= lı́m −1
x→0 x4
L0 H 8x3 + 2x − 2sen x cos x
= lı́m −1
x→0 4x3
L0 H 24x2 − 2(cos2 x − sen2 x)
= lı́m −1
x→0 12x2
L0 H 48x + 8 sen x cos x
= lı́m −1
x→0 24x
L0 H 48 + 8(cos2 x − sen2 x)
= lı́m −1
x→0 24
56 4
= −1=
24 3

4.11. Aplicaciones de la fórmula de Taylor

Estudio local de f
Supongamos que la función f es derivable n + 1 veces en el punto x = a, y tenemos
4.11. APLICACIONES DE LA FÓRMULA DE TAYLOR 119

f 0 (a) = f 00 (a) = · · · = f (n−1) (a) = 0 y f (n) (a) 6= 0

entonces:

1. n es impar y f (n) (a) > 0, entonces f es creciente en a.

2. n es impar y f (n) (a) < 0, entonces f es decreciente en a.

3. n es par y f (n) (a) > 0, entonces f tiene un mı́nimo en a.

4. n es par y f (n) (a) < 0, entonces f tiene un máximo en a.

además en los casos 1 y 2, f tiene un punto de inflexión en a, en el caso 3, f es convexa

en a y en el caso 4, f es cóncava en a.

Ejemplo: Estudiar el comportamiento de la función f (x) = ex (x − 1)3 en el punto x = 1.

Las derivadas primera y segunda de f son

f 0 (x) = ex (x − 1)3 + 3ex (x − 1)2

f 00 (x) = ex (x − 1)3 + 6ex (x − 1)2 + 6ex (x − 1)

y ambas se anulan en x = 1. Hacemos la tercera derivada y tenemos

f 000 (x) = ex (x − 1)3 + 9ex (x − 1)2 + 18ex (x − 1) + 6ex

y se tiene f 000 (1) = 6e > 0, entonces f tiene un punto de inflexión en x = 1, siendo además
creciente en ese punto.
120 CAPÍTULO 4. ANÁLISIS

4.12. Cálculo de primitivas

Consideremos una función f : [a, b] → R continua en [a, b]. Decimos que F es una primitiva
de f en [a, b] cuando F es derivable y se cumple F 0 (x) = f (x) para todo x ∈ [a, b].

Observemos que si F es una primitiva de f en [a, b], entonces cualquier primitiva de f es

de la forma F (x) + C, donde C es una constante.

Nosotros nos ocuparemos del cálculo de primitivas utilizando las técnicas básicas. Este
cálculo, en general, puede ser muy complicado e incluso irresoluble, ya que lo que busca-
mos es expresar la primitiva como suma, producto y composiciń de un número finito de
funciones racionales, logarı́tmicas, exponenciales o trigonométricas y sus inversas. Esto no
2
es posible, por ejemplo con funciones como e−x y senx x .

El primer paso que daremos surge de manera natural si aplicamos las tablas de derivación
en sentido inverso, por ejemplo:

Z Z Z Z Z
5 x 5
ex dx
 
3 + x + cos x + e dx = 3 dx + x dx + cos x dx +

A continuación se muestra una pequeña tabla con las integrales inmediatas, más fáciles

n+1
xn dx = xn+1 (n 6= −1)
R R
C dx = Cx
x−1 dx = log |x|
R R
ex dx = ex
ax
R R
ax dx = log a
sen x dx = − cos x
R R
cos x dx = sen x tg x dx = − log | cos x|
dx
dx = a1 arctg √
x dx 1
log x−a
R R
x2 +a2 a x2 −a2
dx = 2a x+a
(a 6= 0)
√ dx dx x
R R
x2 +a2
= log |x + x 2 + a2 | √
a2 −x2
= arcsen a

Integración por cambio de variable

Supongamos una función ϕ diferenciable, tal que

x = ϕ(t)

entonces

Z Z
f (x)dx = f [ϕ(t)]ϕ0 (t)dt
4.12. CÁLCULO DE PRIMITIVAS 121

Ejemplo: Para calcular la integral

√
Z
x3x−1

√
hacemos t = 3
x − 1, de donde resulta x = t3 + 1 y dx = 3t2 dt y tenemos

√
Z Z Z
3
x x−1 = (t + 1) · t · 3t · dt = 3 (t6 + t3 )dt
3 2

 7
t4

t 3 4 3 7
= 3 + = (x − 1) 3 + (x − 1) 3 .
7 4 4 7

Ejemplo: Calculemos Z
dx
ex +1
Si hacemos el cambio x = − log t, resulta t = e−x y dx = − 1t dt. Entonces

− 1t dt
Z Z Z
dx dt
x
= 1 =−
e +1 t
+1 t+1
= − log |t + 1| = − log |e−x + 1|.

En el ejemplo anterior, parece que lo natural hubiese sido hacer el cambio x = log t,
la razón de elegir el signo menos, es para que una vez hecho el cambio quede en el
denominador 1t + 1 que simplifica más la fracción.

Integración por partes

Si u y v son funciones con derivadas continuas en un intervalo, entonces
Z Z
u dv = u v − v du

Ejemplo: Calculemos la integral Z

log xdx

hacemos

u = log x ⇒ du = dx

x
dv = dx ⇒ v = x
y entonces
122 CAPÍTULO 4. ANÁLISIS

Z Z
log xdx = x log x − dx = x log x − x.

Ejemplo: La integral

Z
x log xdx
R
puede resolverse teniendo en cuenta, como acabamos de ver, que log xdx = x log x − x
haciendo


u=x ⇒ du =R dx
dv = log xdx ⇒ v = log xdx = x log x − x

Integración de funciones racionales

En este apartado estudiamos el método para resolver integrales

Z
P (x)
dx
Q(x)
donde P y Q son polinomios con coeficientes reales. Lo primero que tenemos que obsevar
es el grado de ambos. Si el grado del numerador es mayor que el del denominador se
efectúa la división

P (x) R(x)
= C(x) +
Q(x) Q(x)
y el problema se reduce a calcular la integral
Z
R(x)
dx
Q(x)

donde el grado del numerador es estrictamente menor que el del denominador.

Definición: Llamamos fracción simple a un cociente del tipo

A Mx + N
α
o
(x − a) ((x − p)2 + q 2 )α
donde A, M , N , a, p y q son números reales y α es un número natural.
4.12. CÁLCULO DE PRIMITIVAS 123

Teorema: Todo cociente de polinomios

R(x)
Q(x)
con gr(R(x)) < gr(Q(x)), se puede expresar de la siguiente forma

R(x) A1 Aα
= + ··· + + ···
Q(x) x−a (x − a)α
B1 Bβ
+ ··· + + ···
x−b (x − b)β
M1 x + N1 Mγ x + Nγ
+ 2 2
+ ··· + + ···
(x − m) + n ((x − m)2 + n2 )γ
S1 x + T1 Sδ x + Tδ
+ 2 2
+ ··· +
(x − s) + t ((x − s)2 + t2 )δ

donde

Q(x) = (x − a)α · · · (x − b)β · · · ((x − m)2 + n2 )γ · · · ((x − s)2 + t2 )δ

y los coeficientes A1 , . . . , Aα , B1 , . . . , Bβ , M1 , N1 , . . . , Mγ , Nγ , . . . , S1 , T1 , . . . , Sδ Tδ se cal-
culan por el método de los coeficientes indeterminados.

Esta descomposición facilita como vermos en algunos ejemplos sencillos, el cálculo de la

primitiva.

Ejemplo: Calculemos la integral

Z
dx
x(x − 1)2

Por el metodo de descomposición en fracciones simples tenemos que

1 A B1 B2
2
= + +
x(x − 1) x x − 1 (x − 1)2
y para determinar los coeficientes multiplicamos ambos miembros de la igualdad por
x(x − 1)2 , ası́ resulta
124 CAPÍTULO 4. ANÁLISIS

1 = A(x − 1)2 + B1 x(x − 1) + B2 x

y dando valores convenientes a x, e este caso, por ejemplo x = 0, 1, 2, se tiene el sistema

1 = A
1 = B2
1 = A + 2B1 + 2B2
e inmediatamente los valores A = 1, B1 = −1 y B2 = 1. Ası́ la integral queda de la
siguiente forma

−1
Z Z Z Z
dx 1 1
= dx + dx + dx
x(x − 1)2 x x−1 (x − 1)2
1
= log |x| − log |x − 1| − .
x−1

Ejemplo: Calculemos la integral

x4 − 6x3 + 12x2 + 6
Z
dx
x3 − 6x2 + 12x − 8

Tenemos

x4 − 6x3 + 12x2 + 6
Z Z  
8x + 6
dx = x+ 3 dx
x3 − 6x2 + 12x − 8 x − 6x2 + 12x − 8
Z Z
1 2 8 22
= x + 2
+
2 (x − 2) (x − 2)3
1 2 8 11
= x − − .
2 x − 2 (x − 2)2
4.13. LA INTEGRAL DE RIEMANN 125

4.13. La integral de Riemann

Particiones de un intervalo
Dado un intervalo cerrado [a, b], se dice que el conjunto finito

P = {x0 , x1 , · · · , xn }
donde
a = x0 < x1 < · · · < xn = b
es una partición del intervalo [a, b]. Denotaremos P([a, b]) al conjunto de todas las parti-
ciones del intervalo [a, b]. Dadas dos particiones P1 , P2 ∈ P([a, b]), decimos que P1 es más
fina que P2 si P1 contiene a P2 , esto lo escribiremos P2 < P1 .

Se llama diámetro de la partición P , al número

d(P ) = máx {|xi − xi−1 | : i = 1, 2, · · · , n}

Ejemplo: Consideramos el intervalo [a, b]. Sean

b−a
xi = a + i, i = 0, 1, 2, · · · , n
n
Los números xi forman una partición de [a, b] de diámetro 1/n cuyos subintervalos tienen
la misma amplitud.

Sumas superior e inferior

Sea f una función real definida en el intervalo [a, b] y acotada. Sea P una partición de
[a, b]. Se llama suma superior de Riemann correspondiene a f y P , al número
n
X
S(f, P ) = Mi (xi − xi−1 )
i=1

donde Mi = sup {f (x), x ∈ [xi−1 , xi ]}. Análogamente se llama suma inferior de Riemann
correspondiene a f y P , al número
n
X
S(f, P ) = mi (xi − xi−1 )
i=1

donde mi = ı́nf {f (x), x ∈ [xi−1 , xi ]}

126 CAPÍTULO 4. ANÁLISIS

En las figuras anteriores se pueden ver las sumas superior e inferiort para la función sen x
en el intervalo [0, π] usando una partición de 10 subintervalos de la misma amplitud.

Propiedades de las sumas superior e inferior

Sean P, P1 , P2 ∈ P([a, b]), entonces:

1. S(f, P ) > S(f, P ).

2. Si P2 < P1 , entonces
a) S(f, P1 ) 6 S(f, P2 )
b) S(f, P1 ) > S(f, P2 )
3. S(f, P1 ) > S(f, P2 )


De la propiedad 3 se deduce que los conjuntos S(f, P ) P ∈P([a,b]) y {S(f, P )}P ∈P([a,b])
tienen cota inferior y superior, y por lo tanto, existen los números

Z b 
Z b
f = ı́nf S(f, P ) P ∈P([a,b])
y f = ı́nf {S(f, P )}P ∈P([a,b])
a a

Estos dos números se llaman integral superior e inferior de Riemann de f en [a, b].

Definición. Sea f : [a, b] → R una función acotada, entonces se dice que f es integrable
(en el sentido Riemann), si se tiene
4.13. LA INTEGRAL DE RIEMANN 127

Z b Z b
f= f
a a

este valor se escribe Z b

f.
a

Rb
En el caso particular que f (x) > 0, ∀x ∈ [a, b], la integral a f mide el área de la región
limitada por el eje OX, la gráfica de f y las rectas x = a y x = b.

Criterio de integrabilidad de Cauchy Sea f : [a, b] → R, una función acotada.

Entonces f es integrable en [a, b] si y sólo si dado ε > 0, ∃P ∈ P([a, b]) tal que

S(f, P ) − S(f, P ) < ε

Ejemplo: Veamos que la función de Dirichlet


0 si x ∈ Q
f (x) =
1 si x ∈
/Q

no es integrabe en [0, 1], pues siempre se tiene Mi = 1 y mi = 0 para cualquier partición

P . Entonces se tiene S(f, P ) = 1 y S(f, P ) = 0, ası́ se ve como no se cumple el criterio
Cauchy.

Ejemplo: Consideremos la función parte entera, E(x), que es una función monótona.
Tenemos

Z 4 Z 4
E(x) = E(x) = S(E(x), P ) = 0 + 1 + 2 + 3 = 6
0 0

donde P es una partición que contiene a 0, 1, 2, 3 y 4.

Teorema 4.6. Si f es una función continua en [a, b], entonces f es integrable en [a, b]

La condición del teorema anterior es suficiente, pero no necesaria, como se vio en el ejemplo
de la función parte entera.

Teorema 4.7. Si f es una función monótona en [a, b], entonces f es integrable en [a, b].
128 CAPÍTULO 4. ANÁLISIS

Teorema fundamental del Cálculo

En esta sección se encuentra la relación entre dos temas: el cálculo diferencial y el cálculo
integral, que se han desarrollado independientemente.

Consideremos una función f : [a, b] → R, integrable en ese intervalo. entoncesRtiene sentido

x
considerar la función F : [a, b] → R, definida del siguiente modo: F (x) = a f (t)dt. La
función F es continua en [a, b] y se llama función integral de f en el intervalo [a, b].

Rx
Teorema 4.8 (Fundamental del Cálculo)). Sea f integrable en [a, b] y F (x) = a f (t)dt.
Entonces, si f es continua en [a, b], se tiene que F es dervable en (a, b) y su derivada es
F 0 (x) = f (x).

Teorema 4.9 (Regla de Barrow). Sea f integrable en [a, b] y sea G una primitiva de f ,
en [a, b]. Entonces
Z b
f (x)dx = G(b) − G(a)
a

Rx 2 +1)
Ejemplo: Dada la función F (x) = 0
ecos (t dt, la derivada en el intervalo (0, 2π) de F
2
es F 0 (x) = ecos (x +1) .

Ejemplo: El área A de la región limitada por la gráfica de la función sen x, ej eje 0X y

las rectas x = 0 y x = π es
Z π
A= sen x dx = [− cos x]π0 = 1 + 1 = 2
0

Ejemplo: Calcular la integral definida

Z 3
(|x| + xex ) dx.
−3

Solución: Podemos escribir la integral del siguiente modo:

4.13. LA INTEGRAL DE RIEMANN 129

Z 3 Z 0 Z 3
x x
(|x| + xe ) dx = (−x + xe ) dx + (x + xex ) dx
−3 −3 0
Z 0 Z 0 Z 3 Z 3
x
= −xdx + xe dx + xdx + xex dx
−3 −3 0 0

R
Calculamos la integral indefinida xex dx, por partes:


u=x ⇒ du =R dx
dv = e dx ⇒ v = ex dx = ex
x

entonces, tenemos

Z Z
x x
xe dx = xe − ex dx = xex − ex .

Aplicamos la regla de Barrow y tenemos

3 0  2 3
−x2
Z 
x 0 x
(|x| + xe ) dx =x x
+ [xe − e ]−3 + + [xex − ex ]30
−3 2 −3 2 0
9 9
= 0 + − 1 + 4e−3 + − 0 + 2e3 + 1 = 9 + 4e−3 + 2e3 .
2 2
Índice alfabético

adjunto, 27 menor, 27
aplicación lineal, 45 máximo, 109
imagen de, 45 mı́nimo, 109
núcleo de, 45
número
cofactor, 27 e, 100
número entero, 5
derivada
número natural, 5
de una función en un punto, 107
número racional, 5
ecuación
lineal, 17 operaciones
extremo, 109 elementales, 18

función polinomio caracterı́stico, 64

continua, 105
de Dirichlet, 95 Regla
exponencial, 97 de L’Hôpital, 111
logarı́tmica, 98 Regla de la cadena, 109
trigonométrica, 96
sistema
lı́mite compatible, 18
de una función en un punto, 102 de ecuaciones lineales, 18
de una sucesión, 99 equivalente, 18
incompatible, 18
matrices sistema de ecuaciones lineales, 22
equivalentes, 22 en forma matricial, 22
producto de, 21 sucesión, 99
suma de, 21
matriz, 21 Teorema
traspuesta, 21 de Cauchy, 111
cuadrada, 21 de Rolle, 110
invertible, 21 de Rouché-Frobenius, 31
ortogonal, 65 del valor medio, 110

131