Análisis de Programación Lineal Paramétrica
Análisis de Programación Lineal Paramétrica
Análisis Paramétrico
Departamento de Matemáticas
Universidad Nacional de Trujillo
28 de diciembre de 2021
Introducción
Introducción
Introducción
BxB + NxN = b
BxB + NxN = b
0 0
Actualizando la tabla y denotando cB yj por zj , tenemos:
BxB + NxN = b
0 0
Actualizando la tabla y denotando cB yj por zj , tenemos:
X 0 0 0
z+ [(zj − cj ) + λ(zj − cj )]xj = cB b + λcB b
j∈R
P
xB + j∈R yj xj = b,
X 0 0 0
z+ [(zj − cj ) + λ(zj − cj )]xj = cB b + λcB b
j∈R
P
xB + j∈R yj xj = b,
X 0 0 0
z+ [(zj − cj ) + λ(zj − cj )]xj = cB b + λcB b
j∈R
P
xB + j∈R yj xj = b,
X 0 0 0
z+ [(zj − cj ) + λ(zj − cj )]xj = cB b + λcB b
j∈R
P
xB + j∈R yj xj = b,
X 0 0 0
z+ [(zj − cj ) + λ(zj − cj )]xj = cB b + λcB b
j∈R
P
xB + j∈R yj xj = b,
X 0 0 0
z+ [(zj − cj ) + λ(zj − cj )]xj = cB b + λcB b
j∈R
P
xB + j∈R yj xj = b,
( )
−(zj − cj ) −(zk − ck )
λ̂ = mı́n 0 0 = 0 0
j∈S zj − cj zk − ck
( )
−(zj − cj ) −(zk − ck )
λ̂ = mı́n 0 0 = 0 0
j∈S zj − cj zk − ck
( )
−(zj − cj ) −(zk − ck )
λ̂ = mı́n 0 0 = 0 0
j∈S zj − cj zk − ck
( )
−(zj − cj ) −(zk − ck )
λ̂ = mı́n 0 0 = 0 0
j∈S zj − cj zk − ck
( )
−(zj − cj ) −(zk − ck )
λ̂ = mı́n 0 0 = 0 0
j∈S zj − cj zk − ck
( )
−(zj − cj ) −(zk − ck )
λ̂ = mı́n 0 0 = 0 0
j∈S zj − cj zk − ck
Ejemplo 1
Considere el siguiente problema:
Minimizar −x1 − 3x2
Sujeto a x1 + x2 ≤ 6
−x1 + 2x2 ≤ 6
x1 ≥ 0, x2 ≥ 0
Ejemplo 1
Considere el siguiente problema:
Minimizar −x1 − 3x2
Sujeto a x1 + x2 ≤ 6
−x1 + 2x2 ≤ 6
x1 ≥ 0, x2 ≥ 0
Se desea encontrar las soluciones óptimas y los valores objetivos
óptimos de la clase de problemas cuya función objetivo es
(−1 + 2λ, −3 + λ) para λ ≥ 0; es decir, perturbando el vector de
costos a lo largo del vector (2, 1).
Ejemplo 1
Considere el siguiente problema:
Minimizar −x1 − 3x2
Sujeto a x1 + x2 ≤ 6
−x1 + 2x2 ≤ 6
x1 ≥ 0, x2 ≥ 0
Se desea encontrar las soluciones óptimas y los valores objetivos
óptimos de la clase de problemas cuya función objetivo es
(−1 + 2λ, −3 + λ) para λ ≥ 0; es decir, perturbando el vector de
costos a lo largo del vector (2, 1).
Solución.
Ejemplo 1
Considere el siguiente problema:
Minimizar −x1 − 3x2
Sujeto a x1 + x2 ≤ 6
−x1 + 2x2 ≤ 6
x1 ≥ 0, x2 ≥ 0
Se desea encontrar las soluciones óptimas y los valores objetivos
óptimos de la clase de problemas cuya función objetivo es
(−1 + 2λ, −3 + λ) para λ ≥ 0; es decir, perturbando el vector de
costos a lo largo del vector (2, 1).
Solución.
Primero resolvemos el problema con λ = 0, donde x3 y x4 son
variables de holgura. La tabla óptima para λ = 0 es la siguiente:
Ejemplo 1
Considere el siguiente problema:
Minimizar −x1 − 3x2
Sujeto a x1 + x2 ≤ 6
−x1 + 2x2 ≤ 6
x1 ≥ 0, x2 ≥ 0
Se desea encontrar las soluciones óptimas y los valores objetivos
óptimos de la clase de problemas cuya función objetivo es
(−1 + 2λ, −3 + λ) para λ ≥ 0; es decir, perturbando el vector de
costos a lo largo del vector (2, 1).
Solución.
Primero resolvemos el problema con λ = 0, donde x3 y x4 son
variables de holgura. La tabla óptima para λ = 0 es la siguiente:
z x1 x2 x3 x4 RHS
z 1 0 0 − 35 − 23 −14
2
x1 0 1 0 3 − 13 2
1 1
x2 0 0 1 3 3 4
z x1 x2 x3 x4 RHS
z 1 0 0 − 35 − 23 −14
2
x1 0 1 0 3 − 13 2
1 1
x2 0 0 1 3 3 4
z x1 x2 x3 x4 RHS
z 1 0 0 − 35 − 23 −14
2
x1 0 1 0 3 − 13 2
1 1
x2 0 0 1 3 3 4
z x1 x2 x3 x4 RHS
z 1 0 0 − 35 − 23 −14
2
x1 0 1 0 3 − 13 2
1 1
x2 0 0 1 3 3 4
z x1 x2 x3 x4 RHS
z 1 0 0 − 35 − 23 −14
2
x1 0 1 0 3 − 13 2
1 1
x2 0 0 1 3 3 4
z x1 x2 x3 x4 RHS
z 1 0 0 − 35 − 23 −14
2
x1 0 1 0 3 − 13 2
1 1
x2 0 0 1 3 3 4
−(z3 − c3 ) − (−5/3)
λ̂ = 0 0 = =1
z3 − c3 5/3
Ası́, λ1 = 1 y para λ ∈ [0, 1] la base (a1 , a2 ) permanece óptima.
−(z3 − c3 ) − (−5/3)
λ̂ = 0 0 = =1
z3 − c3 5/3
Ası́, λ1 = 1 y para λ ∈ [0, 1] la base (a1 , a2 ) permanece óptima.
−(z3 − c3 ) − (−5/3)
λ̂ = 0 0 = =1
z3 − c3 5/3
Ası́, λ1 = 1 y para λ ∈ [0, 1] la base (a1 , a2 ) permanece óptima.
−(z3 − c3 ) − (−5/3)
λ̂ = 0 0 = =1
z3 − c3 5/3
Ası́, λ1 = 1 y para λ ∈ [0, 1] la base (a1 , a2 ) permanece óptima.
−(z3 − c3 ) − (−5/3)
λ̂ = 0 0 = =1
z3 − c3 5/3
Ası́, λ1 = 1 y para λ ∈ [0, 1] la base (a1 , a2 ) permanece óptima.
−(z3 − c3 ) − (−5/3)
λ̂ = 0 0 = =1
z3 − c3 5/3
Ası́, λ1 = 1 y para λ ∈ [0, 1] la base (a1 , a2 ) permanece óptima.
−(z3 − c3 ) − (−5/3)
λ̂ = 0 0 = =1
z3 − c3 5/3
Ası́, λ1 = 1 y para λ ∈ [0, 1] la base (a1 , a2 ) permanece óptima.
0 0 2 1
(z4 − c4 ) + λ(z4 − c4 ) = − − λ
3 3
Uceda, R.A. Programación Lineal Paramétrica
Introducción Vector de Costos
Análisis Paramétrico Vector de Recursos
−(z3 − c3 ) − (−5/3)
λ̂ = 0 0 = =1
z3 − c3 5/3
Ası́, λ1 = 1 y para λ ∈ [0, 1] la base (a1 , a2 ) permanece óptima.
0 0 2 1
(z4 − c4 ) + λ(z4 − c4 ) = − − λ
3 3
Uceda, R.A. Programación Lineal Paramétrica
Introducción Vector de Costos
Análisis Paramétrico Vector de Recursos
z x1 x2 x3 x4 RHS
z 1 0 0 − 53 + 35 λ − 23 − 13 λ −14 + 8λ
2
x1 0 1 0 3 − 13 2
1 1
x2 0 0 1 3 3 4
z x1 x2 x3 x4 RHS
z 1 0 0 − 53 + 35 λ − 23 − 13 λ −14 + 8λ
2
x1 0 1 0 3 − 13 2
1 1
x2 0 0 1 3 3 4
z x1 x2 x3 x4 RHS
z 1 0 0 − 53 + 35 λ − 23− 13 λ −14 + 8λ
2
x1 0 1 0 3 − 13 2
1 1
x2 0 0 1 3 3 4
z x1 x2 x3 x4 RHS
z 1 0 0 0 −1 −6
3
x3 0 2 0 1 − 21 3
x2 0 − 12 1 0 1
2 3
z x1 x2 x3 x4 RHS
z 1 0 0 − 53 + 35 λ − 23− 13 λ −14 + 8λ
2
x1 0 1 0 3 − 13 2
1 1
x2 0 0 1 3 3 4
z x1 x2 x3 x4 RHS
z 1 0 0 0 −1 −6
3
x3 0 2 0 1 − 21 3
x2 0 − 12 1 0 1
2 3
z x1 x2 x3 x4 RHS
z 1 0 0 − 53 + 35 λ − 23− 13 λ −14 + 8λ
2
x1 0 1 0 3 − 13 2
1 1
x2 0 0 1 3 3 4
z x1 x2 x3 x4 RHS
z 1 0 0 0 −1 −6
3
x3 0 2 0 1 − 21 3
x2 0 − 12 1 0 1
2 3
!
3 5
z1 − c1 = cB y1 − c1 = (0, 3) 2 +1=
− 12 2
!
− 21 3
z4 − c4 = cB y4 − c4 = (0, 3) 1 −0=−
2 2
!
3 5
z1 − c1 = cB y1 − c1 = (0, 3) 2 +1=
− 12 2
!
− 21 3
z4 − c4 = cB y4 − c4 = (0, 3) 1 −0=−
2 2
!
0 0 0 0
3 5
z1 − c1 = cB y1 − c1 = (0, 1) 2 −2=−
1
−2 2
!
3 5
z1 − c1 = cB y1 − c1 = (0, 3) 2 +1=
− 12 2
!
− 21 3
z4 − c4 = cB y4 − c4 = (0, 3) 1 −0=−
2 2
!
0 0 0 0
3 5
z1 − c1 = cB y1 − c1 = (0, 1) 2 −2=−
1
−2 2
!
0 0 0 0 − 21 1
z4 − c4 = cB y4 − c4 = (0, 1) 1 −0=
2 2
!
3 5
z1 − c1 = cB y1 − c1 = (0, 3) 2 +1=
− 12 2
!
− 21 3
z4 − c4 = cB y4 − c4 = (0, 3) 1 −0=−
2 2
!
0 0 0 0
3 5
z1 − c1 = cB y1 − c1 = (0, 1) 2 −2=−
1
−2 2
!
0 0 0 0 − 21 1
z4 − c4 = cB y4 − c4 = (0, 1) 1 −0=
2 2
Por tanto, los precios sombra para las variables no básicas x1 y x4
serán:
!
3 5
z1 − c1 = cB y1 − c1 = (0, 3) 2 +1=
− 12 2
!
− 21 3
z4 − c4 = cB y4 − c4 = (0, 3) 1 −0=−
2 2
!
0 0 0 0
3 5
z1 − c1 = cB y1 − c1 = (0, 1) 2 −2=−
1
−2 2
!
0 0 0 0 − 21 1
z4 − c4 = cB y4 − c4 = (0, 1) 1 −0=
2 2
Por tanto, los precios sombra para las variables no básicas x1 y x4
serán:
0 0 5 5
(z1 − c1 ) + λ(z1 − c1 ) = − λ
2 2
0 0 3 1
(z4 − c4 ) + λ(z4 − c4 ) = − + λ
2 2
0 0 5 5
(z1 − c1 ) + λ(z1 − c1 ) = − λ
2 2
0 0 3 1
(z4 − c4 ) + λ(z4 − c4 ) = − + λ
2 2
De este modo, para λ ∈ [1, 3] los precios sombra son no positivos y
la base consistente de a3 y a2 es óptima.
0 0 5 5
(z1 − c1 ) + λ(z1 − c1 ) = − λ
2 2
0 0 3 1
(z4 − c4 ) + λ(z4 − c4 ) = − + λ
2 2
De este modo, para λ ∈ [1, 3] los precios sombra son no positivos y
la base consistente de a3 y a2 es óptima.
0 0 5 5
(z1 − c1 ) + λ(z1 − c1 ) = − λ
2 2
0 0 3 1
(z4 − c4 ) + λ(z4 − c4 ) = − + λ
2 2
De este modo, para λ ∈ [1, 3] los precios sombra son no positivos y
la base consistente de a3 y a2 es óptima.
0 0 5 5
(z1 − c1 ) + λ(z1 − c1 ) = − λ
2 2
0 0 3 1
(z4 − c4 ) + λ(z4 − c4 ) = − + λ
2 2
De este modo, para λ ∈ [1, 3] los precios sombra son no positivos y
la base consistente de a3 y a2 es óptima.
0 0 5 5
(z1 − c1 ) + λ(z1 − c1 ) = − λ
2 2
0 0 3 1
(z4 − c4 ) + λ(z4 − c4 ) = − + λ
2 2
De este modo, para λ ∈ [1, 3] los precios sombra son no positivos y
la base consistente de a3 y a2 es óptima.
0 0 5 5
(z1 − c1 ) + λ(z1 − c1 ) = − λ
2 2
0 0 3 1
(z4 − c4 ) + λ(z4 − c4 ) = − + λ
2 2
De este modo, para λ ∈ [1, 3] los precios sombra son no positivos y
la base consistente de a3 y a2 es óptima.
0 0 5 5
(z1 − c1 ) + λ(z1 − c1 ) = − λ
2 2
0 0 3 1
(z4 − c4 ) + λ(z4 − c4 ) = − + λ
2 2
De este modo, para λ ∈ [1, 3] los precios sombra son no positivos y
la base consistente de a3 y a2 es óptima.
0 0 5 5
(z1 − c1 ) + λ(z1 − c1 ) = − λ
2 2
0 0 3 1
(z4 − c4 ) + λ(z4 − c4 ) = − + λ
2 2
De este modo, para λ ∈ [1, 3] los precios sombra son no positivos y
la base consistente de a3 y a2 es óptima.
z x1 x2 x3 x4 RHS
5
z 1 2 − 52 λ 0 0 − 23 + 12 λ −9 + 3λ
3
x3 0 2 0 1 − 12 3
x2 0 − 12 1 0 1
2 3
z x1 x2 x3 x4 RHS
5
z 1 2 − 52 λ 0 0 − 23 + 12 λ −9 + 3λ
3
x3 0 2 0 1 − 12 3
x2 0 − 12 1 0 1
2 3
z x1 x2 x3 x4 RHS
5
z 1 2 − 52 λ 0 0 − 23 + 12 λ −9 + 3λ
3
x3 0 2 0 1 − 12 3
x2 0 − 12 1 0 1
2 3
z x1 x2 x3 x4 RHS
z 1 −5 0 0 0 0
x3 0 1 1 1 0 6
x4 0 −1 2 0 1 6
z x1 x2 x3 x4 RHS
5
z 1 2 − 52 λ 0 0 − 23 + 12 λ −9 + 3λ
3
x3 0 2 0 1 − 12 3
x2 0 − 12 1 0 1
2 3
z x1 x2 x3 x4 RHS
z 1 −5 0 0 0 0
x3 0 1 1 1 0 6
x4 0 −1 2 0 1 6
z x1 x2 x3 x4 RHS
5
z 1 2 − 52 λ 0 0 − 23 + 12 λ −9 + 3λ
3
x3 0 2 0 1 − 12 3
x2 0 − 12 1 0 1
2 3
z x1 x2 x3 x4 RHS
z 1 −5 0 0 0 0
x3 0 1 1 1 0 6
x4 0 −1 2 0 1 6
0 0 0 0
z1 − c1 = cB y1 − c1 = −2
0 0 0 0
z2 − c2 = cB y2 − c2 = −1
0 0 0 0
z1 − c1 = cB y1 − c1 = −2
0 0 0 0
z2 − c2 = cB y2 − c2 = −1
0 0 0 0
z1 − c1 = cB y1 − c1 = −2
0 0 0 0
z2 − c2 = cB y2 − c2 = −1
0
Supongamos que el vector de recursos b es sustituido por b + λb ,
con λ ≥ 0.
0
Supongamos que el vector de recursos b es sustituido por b + λb ,
con λ ≥ 0.
0
Supongamos que el vector de recursos b es sustituido por b + λb ,
con λ ≥ 0.
0
Supongamos que el vector de recursos b es sustituido por b + λb ,
con λ ≥ 0.
z + (cB B −1 N − cN )xN = cB B −1 b
xB + B −1 NxN = B −1 b
z + (cB B −1 N − cN )xN = cB B −1 b
xB + B −1 NxN = B −1 b
donde cB B −1 N − cN ≤ 0.
z + (cB B −1 N − cN )xN = cB B −1 b
xB + B −1 NxN = B −1 b
donde cB B −1 N − cN ≤ 0.
0
Si se sustituye b por b + λb , el vector cB B −1 N − cN no se
afectará; es decir, no se afecta la factibilidad dual.
z + (cB B −1 N − cN )xN = cB B −1 b
xB + B −1 NxN = B −1 b
donde cB B −1 N − cN ≤ 0.
0
Si se sustituye b por b + λb , el vector cB B −1 N − cN no se
afectará; es decir, no se afecta la factibilidad dual.
0
El único cambio es sustituir B −1 b por B −1 (b + λb ) y, en
consecuencia, el valor objetivo se convierte en
z + (cB B −1 N − cN )xN = cB B −1 b
xB + B −1 NxN = B −1 b
donde cB B −1 N − cN ≤ 0.
0
Si se sustituye b por b + λb , el vector cB B −1 N − cN no se
afectará; es decir, no se afecta la factibilidad dual.
0
El único cambio es sustituir B −1 b por B −1 (b + λb ) y, en
consecuencia, el valor objetivo se convierte en
0
cB B −1 (b + λb )
z + (cB B −1 N − cN )xN = cB B −1 b
xB + B −1 NxN = B −1 b
donde cB B −1 N − cN ≤ 0.
0
Si se sustituye b por b + λb , el vector cB B −1 N − cN no se
afectará; es decir, no se afecta la factibilidad dual.
0
El único cambio es sustituir B −1 b por B −1 (b + λb ) y, en
consecuencia, el valor objetivo se convierte en
0
cB B −1 (b + λb )
0
Mientras que B −1 (b + λb ) sea no negativo, la actual base sigue
siendo óptima.
z + (cB B −1 N − cN )xN = cB B −1 b
xB + B −1 NxN = B −1 b
donde cB B −1 N − cN ≤ 0.
0
Si se sustituye b por b + λb , el vector cB B −1 N − cN no se
afectará; es decir, no se afecta la factibilidad dual.
0
El único cambio es sustituir B −1 b por B −1 (b + λb ) y, en
consecuencia, el valor objetivo se convierte en
0
cB B −1 (b + λb )
0
Mientras que B −1 (b + λb ) sea no negativo, la actual base sigue
siendo óptima.
Ejemplo 2
Considere el siguiente problema:
Minimizar −x1 − 3x2
Sujeto a x1 + x2 ≤ 6
−x1 + 2x2 ≤ 6
x1 ≥ 0, x2 ≥ 0
Ejemplo 2
Considere el siguiente problema:
Minimizar −x1 − 3x2
Sujeto a x1 + x2 ≤ 6
−x1 + 2x2 ≤ 6
x1 ≥ 0, x2 ≥ 0
Se desea encontrar la solución óptima y las bases óptimas cuando
se perturba el vector de recursos a lo largo de la dirección
(−1, 1)T , es decir, si b = (6, 6)T se sustituye por
0
b + λb = (6, 6)T + λ(−1, 1)T .
Ejemplo 2
Considere el siguiente problema:
Minimizar −x1 − 3x2
Sujeto a x1 + x2 ≤ 6
−x1 + 2x2 ≤ 6
x1 ≥ 0, x2 ≥ 0
Se desea encontrar la solución óptima y las bases óptimas cuando
se perturba el vector de recursos a lo largo de la dirección
(−1, 1)T , es decir, si b = (6, 6)T se sustituye por
0
b + λb = (6, 6)T + λ(−1, 1)T .
Solución.
Ejemplo 2
Considere el siguiente problema:
Minimizar −x1 − 3x2
Sujeto a x1 + x2 ≤ 6
−x1 + 2x2 ≤ 6
x1 ≥ 0, x2 ≥ 0
Se desea encontrar la solución óptima y las bases óptimas cuando
se perturba el vector de recursos a lo largo de la dirección
(−1, 1)T , es decir, si b = (6, 6)T se sustituye por
0
b + λb = (6, 6)T + λ(−1, 1)T .
Solución.
La solución óptima para λ = 0 se muestra a continuación, siendo
x3 y x4 variables de holgura.
Ejemplo 2
Considere el siguiente problema:
Minimizar −x1 − 3x2
Sujeto a x1 + x2 ≤ 6
−x1 + 2x2 ≤ 6
x1 ≥ 0, x2 ≥ 0
Se desea encontrar la solución óptima y las bases óptimas cuando
se perturba el vector de recursos a lo largo de la dirección
(−1, 1)T , es decir, si b = (6, 6)T se sustituye por
0
b + λb = (6, 6)T + λ(−1, 1)T .
Solución.
La solución óptima para λ = 0 se muestra a continuación, siendo
x3 y x4 variables de holgura.
z x1 x2 x3 x4 RHS
z 1 0 0 − 35 − 23 −14
2
x1 0 1 0 3 − 13 2
1 1
x2 0 0 1 3 3 4
z x1 x2 x3 x4 RHS
z 1 0 0 − 35 − 23 −14
2
x1 0 1 0 3 − 13 2
1 1
x2 0 0 1 3 3 4
z x1 x2 x3 x4 RHS
z 1 0 0 − 35 − 23 −14
2
x1 0 1 0 3 − 13 2
1 1
x2 0 0 1 3 3 4
z x1 x2 x3 x4 RHS
z 1 0 0 − 35 − 23 −14
2
x1 0 1 0 3 − 13 2
1 1
x2 0 0 1 3 3 4
b1 2
λ1 = 0 = =2
−b 1 −(−1)
z x1 x2 x3 x4 RHS
z 1 0 0 − 35 − 23 −14
2
x1 0 1 0 3 − 13 2
1 1
x2 0 0 1 3 3 4
b1 2
λ1 = 0 = =2
−b 1 −(−1)
z x1 x2 x3 x4 RHS
z 1 0 0 − 53 − 23 −14 + λ
2
x1 0 1 0 3 − 13 2−λ
1 1
x2 0 0 1 3 3 4
z x1 x2 x3 x4 RHS
z 1 0 0 − 53 − 23 −14 + λ
2
x1 0 1 0 3 − 13 2−λ
1 1
x2 0 0 1 3 3 4
z x1 x2 x3 x4 RHS
z 1 −2 0 −3 0 −12
x4 0 -3 0 −2 1 0
x2 0 1 1 1 0 4
z x1 x2 x3 x4 RHS
z 1 0 0 − 53 − 23 −14 + λ
2
x1 0 1 0 3 − 13 2−λ
1 1
x2 0 0 1 3 3 4
z x1 x2 x3 x4 RHS
z 1 −2 0 −3 0 −12
x4 0 -3 0 −2 1 0
x2 0 1 1 1 0 4
z x1 x2 x3 x4 RHS
z 1 0 0 − 53 − 23 −14 + λ
2
x1 0 1 0 3 − 13 2−λ
1 1
x2 0 0 1 3 3 4
z x1 x2 x3 x4 RHS
z 1 −2 0 −3 0 −12
x4 0 -3 0 −2 1 0
x2 0 1 1 1 0 4
! ! !
−2 1 6 −6
b= B −1 b = =
1 0 6 6
! ! !
0 0 −2 1 −1 3
b = B −1 b = =
1 0 1 −1
! ! !
−2 1 6 −6
b= B −1 b = =
1 0 6 6
! ! !
0 0 −2 1 −1 3
b = B −1 b = =
1 0 1 −1
! ! !
−2 1 6 −6
b= B −1 b = =
1 0 6 6
! ! !
0 0 −2 1 −1 3
b = B −1 b = =
1 0 1 −1
b2 6
λ2 = = =6
−b 2 −(−1)
! ! !
−2 1 6 −6
b= B −1 b = =
1 0 6 6
! ! !
0 0 −2 1 −1 3
b = B −1 b = =
1 0 1 −1
b2 6
λ2 = = =6
−b 2 −(−1)
! ! !
−2 1 6 −6
b= B −1 b = =
1 0 6 6
! ! !
0 0 −2 1 −1 3
b = B −1 b = =
1 0 1 −1
b2 6
λ2 = = =6
−b 2 −(−1)
! ! !
−2 1 6 −6
b= B −1 b = =
1 0 6 6
! ! !
0 0 −2 1 −1 3
b = B −1 b = =
1 0 1 −1
b2 6
λ2 = = =6
−b 2 −(−1)
! ! !
−2 1 6 −6
b= B −1 b = =
1 0 6 6
! ! !
0 0 −2 1 −1 3
b = B −1 b = =
1 0 1 −1
b2 6
λ2 = = =6
−b 2 −(−1)
! ! !
0 −6 3 −6 + 3λ
b + λb = +λ =
6 −1 6−λ
! ! !
0 −6 3 −6 + 3λ
b + λb = +λ =
6 −1 6−λ
z x1 x2 x3 x4 RHS
z 1 −2 0 −3 0 −18 + 3λ
x4 0 -3 0 −2 1 −6 + 3λ
x2 0 1 1 1 0 6−λ
! ! !
0 −6 3 −6 + 3λ
b + λb = +λ =
6 −1 6−λ
z x1 x2 x3 x4 RHS
z 1 −2 0 −3 0 −18 + 3λ
x4 0 -3 0 −2 1 −6 + 3λ
x2 0 1 1 1 0 6−λ
! ! !
0 −6 3 −6 + 3λ
b + λb = +λ =
6 −1 6−λ
z x1 x2 x3 x4 RHS
z 1 −2 0 −3 0 −18 + 3λ
x4 0 -3 0 −2 1 −6 + 3λ
x2 0 1 1 1 0 6−λ