Trial of new calculation in which coherence 21 is completely out of the game
+ remove upperdiag (no improvement)
+ Runge-Kutta reintegration