Bei Werten ab etwa 10 für a und b führt diese naive Implementierung
mit double-Genauigkeit zu unakzeptablen bis absurden Rechenfehlern,
verursacht durch extreme Unterschiede der Größenordnungen der Summanden.
Nicht nachmachen!
Weil die Werte von Βx(a,b) rasch sehr klein werden, wird auch der normalisierte Wert Ix(a,b)=Βx(a,b)/Β1(a,b) angezeigt.