31     puts(
"############ Test the LU decomposition algorithm ###############");
 
   32     matrix_t A[5][5] = { { 0.8147, 0.1576, 0.6557, 0.7060, 0.4387 },
 
   33                          { 0.9058, 0.9706, 0.0357, 0.0318, 0.3816 },
 
   34                          { 0.1270, 0.9572, 0.8491, 0.2769, 0.7655 },
 
   35                          { 0.9134, 0.4854, 0.9340, 0.0462, 0.7952 },
 
   36                          { 0.6324, 0.8003, 0.6787, 0.0971, 0.1869 }, };
 
   51         { 0.4387,  0.6797,  0.5060,  0.2435,  0.9172,  0.1299,  0.2630,  0.9961,  0.2599,  0.1450,
 
   53         { 0.3816,  0.6551,  0.6991,  0.9293,  0.2858,  0.5688,  0.6541,  0.0782,  0.8001,  0.8530,
 
   55         { 0.7655,  0.1626,  0.8909,  0.3500,  0.7572,  0.4694,  0.6892,  0.4427,  0.4314,  0.6221,
 
   57         { 0.7952,  0.1190,  0.9593,  0.1966,  0.7537,  0.0119,  0.7482,  0.1067,  0.9106,  0.3510,
 
   59         { 0.1869,  0.4984,  0.5472,  0.2511,  0.3804,  0.3371,  0.4505,  0.9619,  0.1818,  0.5132,
 
   61         { 0.4898,  0.9597,  0.1386,  0.6160,  0.5678,  0.1622,  0.0838,  0.0046,  0.2638,  0.4018,
 
   63         { 0.4456,  0.3404,  0.1493,  0.4733,  0.0759,  0.7943,  0.2290,  0.7749,  0.1455,  0.0760,
 
   65         { 0.6463,  0.5853,  0.2575,  0.3517,  0.0540,  0.3112,  0.9133,  0.8173,  0.1361,  0.2399,
 
   67         { 0.7094,  0.2238,  0.8407,  0.8308,  0.5308,  0.5285,  0.1524,  0.8687,  0.8693,  0.1233,
 
   69         { 0.7547,  0.7513,  0.2543,  0.5853,  0.7792,  0.1656,  0.8258,  0.0844,  0.5797,  0.1839,
 
   71         { 0.2760,  0.2551,  0.8143,  0.5497,  0.9340,  0.6020,  0.5383,  0.3998,  0.5499,  0.2400,