A matrix with sums of Catalan numbers -- LU-decomposition and determinant

Following Benjamin et al., a matrix with entries being sums of two neighbouring Catalan numbers is considered. Its LU-decomposition is given, by guessing the results and later prove it by computer algebra, with lots of human help. Specializing a parameter, the determinant turns out to be a Fibonacci number with odd index, confirming earlier results, obtained back then by combinatorial methods.