Develop a subprogram to achieve the LU decomposition of a banded matrix with l sub diagonals and super diagonals with either the Crout or the Doolittle factorization. Simplify the subprogram for the cases of and a symmetric matrix.