Answer: D
To solve the system using LU \text{LU} LU decomposition, we first write the system in matrix form A x = b Ax = b A x = b :
( 1 1 − 2 1 3 − 1 2 1 − 5 ) ( x 1 x 2 x 3 ) = ( 4 7 7 ) \begin{pmatrix}
1 & 1 & -2 \\
1 & 3 & -1 \\
2 & 1 & -5
\end{pmatrix}
\begin{pmatrix}
x_1 \\ x_2 \\ x_3
\end{pmatrix}
=
\begin{pmatrix}
4 \\ 7 \\ 7
\end{pmatrix} 1 1 2 1 3 1 − 2 − 1 − 5 x 1 x 2 x 3 = 4 7 7
Using Doolittle's algorithm, we set the diagonal elements of L L L to 1 1 1 (L 11 = L 22 = L 33 = 1 L_{11} = L_{22} = L_{33} = 1 L 11 = L 22 = L 33 = 1 ). Thus, we have A = L U A = LU A = LU :
( 1 1 − 2 1 3 − 1 2 1 − 5 ) = ( 1 0 0 L 21 1 0 L 31 L 32 1 ) ( U 11 U 12 U 13 0 U 22 U 23 0 0 U 33 ) \begin{pmatrix}
1 & 1 & -2 \\
1 & 3 & -1 \\
2 & 1 & -5
\end{pmatrix}
=
\begin{pmatrix}
1 & 0 & 0 \\
L_{21} & 1 & 0 \\
L_{31} & L_{32} & 1
\end{pmatrix}
\begin{pmatrix}
U_{11} & U_{12} & U_{13} \\
0 & U_{22} & U_{23} \\
0 & 0 & U_{33}
\end{pmatrix} 1 1 2 1 3 1 − 2 − 1 − 5 = 1 L 21 L 31 0 1 L 32 0 0 1 U 11 0 0 U 12 U 22 0 U 13 U 23 U 33
By multiplying the matrices and equating the elements:
First row: U 11 = 1 U_{11} = 1 U 11 = 1 , U 12 = 1 U_{12} = 1 U 12 = 1 , U 13 = − 2 U_{13} = -2 U 13 = − 2
Second row:
L 21 U 11 = 1 ⟹ L 21 = 1 L_{21}U_{11} = 1 \implies L_{21} = 1 L 21 U 11 = 1 ⟹ L 21 = 1
L 21 U 12 + U 22 = 3 ⟹ 1 ( 1 ) + U 22 = 3 ⟹ U 22 = 2 L_{21}U_{12} + U_{22} = 3 \implies 1(1) + U_{22} = 3 \implies U_{22} = 2 L 21 U 12 + U 22 = 3 ⟹ 1 ( 1 ) + U 22 = 3 ⟹ U 22 = 2
L 21 U 13 + U 23 = − 1 ⟹ 1 ( − 2 ) + U 23 = − 1 ⟹ U 23 = 1 L_{21}U_{13} + U_{23} = -1 \implies 1(-2) + U_{23} = -1 \implies U_{23} = 1 L 21 U 13 + U 23 = − 1 ⟹ 1 ( − 2 ) + U 23 = − 1 ⟹ U 23 = 1
Third row:
L 31 U 11 = 2 ⟹ L 31 = 2 L_{31}U_{11} = 2 \implies L_{31} = 2 L 31 U 11 = 2 ⟹ L 31 = 2
L 31 U 12 + L 32 U 22 = 1 ⟹ 2 ( 1 ) + L 32 ( 2 ) = 1 ⟹ L 32 = − 1 2 L_{31}U_{12} + L_{32}U_{22} = 1 \implies 2(1) + L_{32}(2) = 1 \implies L_{32} = -\frac{1}{2} L 31 U 12 + L 32 U 22 = 1 ⟹ 2 ( 1 ) + L 32 ( 2 ) = 1 ⟹ L 32 = − 2 1
L 31 U 13 + L 32 U 23 + U 33 = − 5 ⟹ 2 ( − 2 ) + ( − 1 2 ) ( 1 ) + U 33 = − 5 ⟹ U 33 = − 1 2 L_{31}U_{13} + L_{32}U_{23} + U_{33} = -5 \implies 2(-2) + \left(-\frac{1}{2}\right)(1) + U_{33} = -5 \implies U_{33} = -\frac{1}{2} L 31 U 13 + L 32 U 23 + U 33 = − 5 ⟹ 2 ( − 2 ) + ( − 2 1 ) ( 1 ) + U 33 = − 5 ⟹ U 33 = − 2 1
Now, we solve the system of equations to find x 1 x_1 x 1 :
x 1 + x 2 − 2 x 3 = 4 x_1 + x_2 - 2x_3 = 4 x 1 + x 2 − 2 x 3 = 4
x 1 + 3 x 2 − x 3 = 7 x_1 + 3x_2 - x_3 = 7 x 1 + 3 x 2 − x 3 = 7
2 x 1 + x 2 − 5 x 3 = 7 2x_1 + x_2 - 5x_3 = 7 2 x 1 + x 2 − 5 x 3 = 7
Subtracting equation (1) from (2) gives:
2 x 2 + x 3 = 3 — (4) 2x_2 + x_3 = 3 \quad \text{--- (4)} 2 x 2 + x 3 = 3 — (4)
Multiplying equation (1) by 2 2 2 and subtracting from (3) gives:
− x 2 − x 3 = − 1 ⟹ x 2 + x 3 = 1 — (5) -x_2 - x_3 = -1 \implies x_2 + x_3 = 1 \quad \text{--- (5)} − x 2 − x 3 = − 1 ⟹ x 2 + x 3 = 1 — (5)
Solving equations (4) and (5) yields x 2 = 2 x_2 = 2 x 2 = 2 and x 3 = − 1 x_3 = -1 x 3 = − 1 . Substituting these back into equation (1):
x 1 + 2 − 2 ( − 1 ) = 4 ⟹ x 1 = 0 x_1 + 2 - 2(-1) = 4 \implies x_1 = 0 x 1 + 2 − 2 ( − 1 ) = 4 ⟹ x 1 = 0
Thus, L 32 = − 1 2 L_{32} = -\frac{1}{2} L 32 = − 2 1 , U 33 = − 1 2 U_{33} = -\frac{1}{2} U 33 = − 2 1 , and x 1 = 0 x_1 = 0 x 1 = 0 .