Crout method Crout type (u_ii == 1) Input a: matrix, n: matrix size Output a: matrix for(k=0;k