This is part of the code of the master thesis of Eugeniu Zaicanu, a student of mine back at Seeburger AG