the solution is obtained by gaussian elimination with partial pivoting, with 2/3n³ + 2n² floating point operations where "n" is 100, the order of the dense matrix "a" that defines the problem. its small size and the lack of software flexibility doesn't allow most modern computers to reach their performance limits. however, it can still be useful to predict performances in numerically intensive user written code using compiler optimization.