I have a question on if there is a way to convert a LU decomposition into a crout decomposition. I have a Ti inspire cx CAS that can very easily return the upper and lower triangular matrices (with the lower containing the diagonal of 1s) using the LU['''matrix'''],lo,up,perm command. I am looking for a method to find the crout decomposition (which has the diagonal of 1s on the upper triangle matrix), either by converting from the Doolittle decomposition matrices or something else entirely. Most problems are 3x3 or 4x4. Thanks.