Naive Berechnung der unvollständigen Beta-Funktion

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.

Βx1(a,b)
Βx2(a,b)
Β1(a,b)
Ix1(a,b)
Ix2(a,b)
min(|summands|)
max(|summands|)
max / min