Ich sehe gerade, dass die beiden Winkel falsch sind.
Phi 1 = 133.5668293844902°
Phi 2 = 61.46435809214322°
Da hab ich mich tatsächlich in der Zuordnung der Ausgabe der Werte vertan. Richtig sind die Werte für die Winkel wie in der ausführlichen Liste der Zwischenwerte angegeben.
Ahh - das hat sich mit deiner Antwort überschnitten. OK, dann suche ich erst mal nicht weiter, ob ich einen Fehler im Programm eingebaut habe und warte erst mal ab.