fint calculated correctly

solver does not yet converge
kinematic force calculated correctly

LAPack function DGESV was modifying the Mass matrix while solving and
outputting the modified mass matrix to replace the one which was
correctly calculated