Use Matlab

10.19 Develop a user-friendly program for LU decomposition, in cluding the capability to evaluate the matrix inverse. Base the pro- gram on Figs. 10.2 and 10.5 LU DECOMPOSITION AND MATRIX INVERSION SUB Ludecomp (a, b, n, tol. x. er) END IF DIM On. S er = 0 CALL Decompose(a, n. tol, o, s, er) IF er -1 THEN END Decompose SUB Pivot (a, o. s. n. k) big = ABS(ao(k),k/so(k)) DOFOR 11k+1. rn CALL Substitute(a, o. n, b. x) END IF END Ludecomp dumy ABS(aoror) IF dumy > big THEN SUB Decompose (a, n, tol, o, s. er) big - dummy END IF ABS(a,1) END DO IF ABS(a, J)s, THEN S,ABS(a j)O END DO END DO DOFOR k = 1, n-1 END Pivot SUB Substitute (a, o. n. b. x) CALL Pivot(a, o. s. n. k) IF ABS(ao(k).ds(k))