Solution Manual Jaan Kiusalaas Numerical — Methods In Engineering With Matlab 2nd 58

b_perm = P*b; y = forwardSub(L, b_perm); x = backSub(U, y); disp(x); ( x \approx [0.7234, -0.6809, -1.1064]^T ) Step 3: Compute inverse using LU decomposition For ( A^-1 ), solve ( A X = I ), column by column, reusing ( L, U, P ):

[ P = \beginbmatrix 0 & 0 & 1 \ 1 & 0 & 0 \ 0 & 1 & 0 \endbmatrix, \quad L = \beginbmatrix 1 & 0 & 0 \ 0.6 & 1 & 0 \ -0.4 & 0.5455 & 1 \endbmatrix, \quad U = \beginbmatrix 5 & 2 & -3 \ 0 & -2.2 & 3.8 \ 0 & 0 & 4.2727 \endbmatrix ] b_perm = P*b; y = forwardSub(L, b_perm); x

(Values are approximate, matching typical pivot choices.) First permute ( b ): ( b' = P b ). Then forward substitution: ( L y = b' ). Then back substitution: ( U x = y ). b_perm = P*b