Dank dir für die Auszeichnung. Deshalb hab noch mal drüber geschaut und festgestellt, dass der Normalisierungsfaktor im Richtungsvektor abhanden gekommen ist. u muss lauten auf
u=1+tan2(α)r(cos(β),sin(β),tan(α))
ich korrigiere das in der Antwort oben...