backup of work so far. Finally, I have solved the fourway-problem. JUmp-times turned uot to be stored incorrectly in simulations of fourway coalescents
Started implementing the recursion-equiations in multiple-precision arithmetics