% Inputs: % Coefficient matrix A % % Outputs: % an lower and upper triangular matrix; % Usage: % [l, u] = LUdecompCrout(a); % Notes: % % ***********************************************************************…