dra = 1! - ra - 2! * DA '- Diese Konstanten sparen spÑter Rechenzeit
  drb = 1! - rB - 2! * db '- sie mÅssen dann nicht immer wieder neu berechnet
  drc = 1! - rc - 2! * DC '- werden

olddecaydiffA = (1 - r_a - 2*D_a) * a + D_a * (a + 1)
aq = s * b * (a^2 + ba) / 
na =