Proving the inverse (if any) of a lower triangular matrix is lower triangular

Let's write $$L^{-1}=[y_1\:\cdots\:y_n],$$ where each $y_k$ is an $n\times 1$ matrix.

Now, by definition, $$LL^{-1}=I=[e_1\:\cdots\:e_n],$$ where $e_k$ is the $n\times 1$ matrix with a $1$ in the $k$th row and $0$s everywhere else. Observe, though, that $$LL^{-1}=L[y_1\:\cdots\:y_n]=[Ly_1\:\cdots\: Ly_n],$$ so $$Ly_k=e_k\qquad(1\leq k\leq n)$$

By the proposition, since $e_k$ has only $0$s above the $k$th row and $L$ is lower triangular and $Ly_k=e_k$, then $y_k$ has only $0$s above the $k$th row. This is true for all $1\leq k\leq n$, so since $$L^{-1}=[y_1\:\cdots\:y_n],$$ then $L^{-1}$ is lower triangular, too.

$$********$$

Here's an alternative (but related) approach.

Observe that a lower triangular matrix is nonsingular if and only if it has all nonzero entries on the diagonal. Let's proceed by induction on $n$. The base case ($n=1$) is simple, as all scalars are trivially "lower triangular". Now, let's suppose that all nonsingular $n\times n$ lower triangular matrices have lower triangular inverses, and let $A$ be any nonsingular $(n+1)\times(n+1)$ lower triangular matrix. In block form, then, we have $$A=\left[\begin{array}{c|c}L & 0_n\\\hline x^T & \alpha\end{array}\right],$$ where $L$ is a nonsingular $n\times n$ lower triangular matrix, $0_n$ is the $n\times 1$ matrix of $0$s, $x$ is some $n\times 1$ matrix, and $\alpha$ is some nonzero scalar. (Can you see why this is true?) Now, in compatible block form, we have $$A^{-1}=\left[\begin{array}{c|c}M & b\\\hline y^T & \beta\end{array}\right],$$ where $M$ is an $n\times n$ matrix, $b,y$ are $n\times 1$ matrices, and $\beta$ some scalar. Letting $I_n$ and $I_{n+1}$ denote the $n\times n$ and $(n+1)\times(n+1)$ identity matrices, respectively, we have $$I_{n+1}=\left[\begin{array}{c|c}I_n & 0_n\\\hline 0_n^T & 1\end{array}\right].$$ Hence, $$\left[\begin{array}{c|c}I_n & 0_n\\\hline 0_n^T & 1\end{array}\right]=I_{n+1}=A^{-1}A=\left[\begin{array}{c|c}ML+by^T & M0_n+b\alpha\\\hline x^TM+\alpha y^T & y^T0_n+\beta\alpha\end{array}\right]=\left[\begin{array}{c|c}ML+by^T & \alpha b\\\hline x^TM+\alpha y^T & \beta\alpha\end{array}\right].$$ Since $\alpha$ is a nonzero scalar and $\alpha b=0_n$, then we must have $b=0_n$. Thus, $$A^{-1}=\left[\begin{array}{c|c}M & 0_n\\\hline y^T & \beta\end{array}\right],$$ and $$\left[\begin{array}{c|c}I_n & 0_n\\\hline 0_n^T & 1\end{array}\right]=\left[\begin{array}{c|c}ML & 0_n\\\hline x^TM+\alpha y^T & \beta\alpha\end{array}\right].$$ Since $ML=I_n$, then $M=L^{-1}$, and by inductive hypothesis, we have that $M$ is then lower triangular. Therefore, $$A^{-1}=\left[\begin{array}{c|c}M & 0_n\\\hline y^T & \beta\end{array}\right]$$ is lower triangular, too, as desired.

Suppose you have an invertible lower-triangular matrix $L$. To find its inverse, you must solve the matrix equation $LX = I$, where $I$ denotes the $n$-by-$n$ identity matrix.

Based on how matrix multiplication works, the $i^{\text{th}}$ column of $LX$ is equal to $L$ times the $i^{\text{th}}$ column of $X$. In order for $LX = I$, it must be that the first $i-1$ entries in the $i^{\text{th}}$ column of $LX$ are all zero. The hint is that you can prove that this implies that the first $i-1$ entries in the $i^{\text{th}}$ column of $X$ must all be zero. To do this, you can explicitly write out your calculation, using your assumption that $L$ is lower-triangular. You'll get a fairly easy system of linear equations to analyze.

In simple form, we can write A = D*(I+L); where A is lower triangular matrix, D is diagonal matrix, I is identity matrix and L is lower triangular with all zeros in diagonal. Since $A^{-1} = (I+L)^{-1}*D^{-1}$ and inverse of D is simply inverse of diagonal element. And for very large n $L^{-n} = 0$ since it is having only lower triangular elements. And we can write $(I+L)^{-1} = I - L + L^2 - L^3 + .... (-1)^n*L^n$ which itself is lower triangular matrix.