Code for Hull-White LU Decomposition

Hello! I'm reading this paper: A direct LU solver for pricing American bond options under Hull–White model . I was wondering if anyone has the MATLAB code for the LU decomposition? I find trying to work this out alone really challenging as an undergraduate, so any kind of help would be really appreciated!