findings

a standard approach is implemented. 2x2 and 3x3 systems can be
solved faster with a simpler approach without loops.
While rather big systems can be solved, it should be noted
that special approaches exist to solve sparse systems and
for those with bands.