use Modelica.Math.Matrices.solve in setSat_T solves #16

(although it should do exactly the same as the manual solution)