debugging colddae_ + newton method

colddae_ uses a self written newton method to compute the shifted time
points

I also changed the selection strategy for the regular system, before we
had U2 = orth(...), but then Example 5 in the noncausal test set caused
some errors.