Expansión de Laplace para calcular Determinantes & derivada de Lanczos

1. Método de Laplace (1772):

Consiste en obtener el valor de un Determinante reduciéndolo a Determinantes de menor orden. Este

| |
3 7 4
D 1= 1 9 6 , y en todos los cálculos mantener la imagen de signos: ¿, (1)
2 5 8

entonces determinamos D1 vía la 1ra. fila:

D1=¿ (2)

donde se observa que en D1 eliminamos la fila y columna de 3 para obtener el 1er. determinante 2 x 2; el
2do. determinante 2 x 2 se obtuvo al cancelar en D1 la fila y columna de 7, finalmente, al suprimir en D1
la fila y columna de 4 resulta el 3er. determinante 2 x 2. Por lo tanto:

D1=3 ( 9 ∙ 8−6∙ 5 ) −7 ( 1∙ 8−2∙ 6 ) + 4 ( 1 ∙5−2∙ 9 ) =102. (3)

Este valor también se logra si ahora realizamos el cálculo vía, por ejemplo, la 2da. columna:

| |
3 7 4
D1= 1 9 6 =¿
2 5 8

¿−7 ( 1∙ 8−2∙ 6 ) + 9 ( 3 ∙ 8−2∙ 4 )−5 ( 3 ∙ 6−1 ∙ 4 )=102.

Nótese que en todo momento tomamos en cuenta los signos indicados en (1). Otro ejemplo, que
desarrollaremos por la 3ra. fila:

| |
10 0 21
D2= 5 −4 13 =¿
2 −7 16

¿ 2 ( 0 ∙13−(−4 ) ∙21 ) +7 (10 ∙ 13−5 ∙ 21 ) +16 ( 10 ∙ (−4 )−5 ∙0 )=−297. (4)

Ahora consideremos la situación 4 x 4, que calcularemos vía la 4ta. fila [porque ahí hay muchos ceros]:

| |
1 2 0 1
0 −1 1 2
D 3= , con la imagen de signos: ¿ (5)
1 1 0 1
0 0 1 0
así:
 D3=¿

| |
1 2 1
¿−1 0 −1 2 éste es el único determinante 3 x 3 que sobrevive porque los demás están
1 1 1

multiplicados por cero. Ahora como ejercicio dejo que calculen este último determinante 3 x 3 como se
mostró para D1 y D2 . El resultado debe ser D3=−2.

Ahora propongo calcular los siguientes Determinantes mediante el Método de Laplace, y luego usar una
Calculadora o algún Paquete para comprobar cada uno de los valores obtenidos:

| | | |
1 −1 2 −2

| || |
5 0 9 1
−1 0 3 8 −11 0
6 0 3 −6 11 −4 7 8
6 −2 1 , −9 3 5, , . (6)
8 −3 −8 7 0 6 −1 3
0 4 5 1 10 0
10 0 5 0 2 0 −3 1

2. Derivada de Lanczos.

A este concepto también se le llama “Derivada Ortogonal” o “Derivada Generalizada”, el cual fue
introducido por Cioranescu, Haslam-Jones y Lanczos. Recuerden que el matemático húngaro Lanczos fue
colaborador de A. Einstein, y ya lo mencionamos en el tema de “Economización de Potencias” mediante
Polinomios de Chebyshev.

En esencia, este enfoque muestra cómo calcular una derivada empleando una integral! La técnica es sencilla,
y de importancia en diversas aplicaciones en Análisis Numérico e Ingeniería, espero les agrade. Para este
tópico deben tener una calculadora o un Paquete que integre, en lo personal me acostumbré a la máquina
Casio fx-991ES.

En la fórmula de Lanczos hay una cantidad que debe ser muy pequeña, y después de muchísimos ejemplos
observamos que se logran buenos resultados si utilizamos 10−4, sin embargo, si alguien desea puede emplear
−15
10 pero no es necesario un valor demasiado pequeño. Entonces:

0.0001
'
f ( x 0 )=1.5 x 10
12
∫ t ∙ f ( t + x 0 ) dt , (7)
−0.0001

lo cual quedará claro si hacemos varios ejemplos:

a). f ( x )=tan x .

De la Voca [o las Tablas] sabemos que f ' ( x )=Sec 2 x , entonces la derivada en x 0=1 es fácil de obtener con
la Calculadora en Modo de Radianes:

' 2
f ( 1 )=Sec 1=3.42551882. (8)

Lo anterior es la manera tradicional de obtener una derivada en determinado punto, ahora lo novedoso es
ejecutar el cálculo mediante la integral (7) propuesta por Lanczos. Antes, note lo siguiente:
 f ( x )=tan x ∴ f ( t + x 0 )=f ( t+1 ) =tan ( t+ 1 ) ,

que se encuentra en el integrando de (7), en consecuencia vuestra Calculadora o Paquete debe computar la
siguiente integral de Lanczos:

0.0001
'
f ( 1 )=1.5 x 10
12
∫ t ∙ tan ( t+ 1 ) dt=3.42551887 , (9)
−0.0001

bastante cercano al valor (8).

b). f ( x )=x 2 .

' '
Es sencilla la derivada de esta función: f ( x )=2 x ∴ f ( 1 )=2. (10)

2
Ahora también se tiene x 0=1, entonces f ( t+ x 0 )=f ( t+1 )=(t +1) que debe aplicarse en (7):

0.0001
'
f ( 1 )=1.5 x 10
12
∫ 2
t (t+1) dt=2, (11)
−0.0001

en acuerdo con el valor exacto (10).

1
c). f ( x )= . Deseamos su derivada en x 0=0.5 .
1+ x

El cálculo directo de la Prepa da:

' −1 −1
f ( x )= 2
∴ f ' ( 0.5 )= 2 =−0.4444 4444 44 … (12)
(1+ x ) 1.5

1 1
y en este caso: f ( t+ x 0 )=f ( t+ 0.5 )= = que debe emplearse en (7):
1+t+ 0.5 t +1.5

0.0001
t
'
f ( 0.5 )=1.5 x 10
12
∫ dt=−0.4444 4444 56 … (13)
−0.0001 t+1.5

muy cercano a (12).

d). f ( x )=ln x . Nos interesa su derivada en x 0=0.25 .

' 1 ' 1
Es conocida la fórmula: f ( x )= ∴ f ( 0.25 )= =4 , (14)
x 0.25

además: f ( t+ x 0 )=f ( t+ 0.25 )=ln (t+ 0.25) y así (7) implica el valor:
 0.0001
'
f ( 0.25 )=1.5 x 10
12
∫ t ∙ ln ( t +0.25 ) dt =4.0000 0012 8 (15)
−0.0001

que puede compararse con el cálculo exacto (14).

Entonces recomiendo realizar los siguientes ejercicios sobre la Derivada de Lanczos, y recuerden que la
Calculadora debe estar en Modo de Radianes:

I. f ( x )=x 4 . Queremos su derivada en x 0=1.23

' 3 ' 3
Por cálculo directo obtenemos: f ( x )=4 x ∴ f ( 1.23 )=4(1.23) =7.4434 68. (16)

y con el enfoque de Lanczos es necesario obtener la integral en (7):

0.0001
'
f ( 1.23 )=1.5 x 10
12
∫ 4
t ∙(t +1.23) dt , (17)
−0.0001

y comparar con (16). Calcule (17).

II. f ( x )=x 2 e x . Necesitamos su derivada en x 0=2.16

Es fácil derivar este producto para deducir que:

f ' ( x )=( x 2 +2 x ) e x ∴ f ' ( 2.16 )=[ ( 2.16 )2 +2 ( 2.16 ) ] e 2.16=77.9153 7454 , (18)

entonces como ejercicio obtenga el valor de la integral de Lanczos:

0.0001
'
f ( 2.16 )=1.5 x 10
12
∫ t (t+ 2.16) e
2 t +2.16
dt , (19)
−0.0001

y compare con (18).

III. f ( x )=x sen x . Nos interesa su derivada en x 0=3.62

En este caso:

' '
f ( x )=sen x+ x cos x ∴ f ( 3.62 )=sen 3.62+3.62 cos 3.62=−3.6739 4580 6 , (20)

entonces con la Calculadora o con algún Paquete obtenga el valor de (7):

0.0001
'
f ( 3.62 )=1.5 x 10
12
∫ t ( t +3.62 ) sen ( t +3.62 ) dt , (21)
−0.0001

y compare con (20).

Al método de Lanczos también se le llama: Diferenciación vía Integración.

En mis Cursos también enseño el proceso inverso: Integración vía Diferenciación,

Es útil saber muchos trucos para evaluar una integral, pero en la vida real aparecen integrales tan difíciles
que es casi imposible realizar un cálculo exacto, por eso deben practicar la Tecla de Integración en su
Calculadora y también el manejo de un Paquete Computacional, lo cual les será útil en diversos Cursos.

Se cuidan.

Prof. José Luis