Oh, der Winkel im Bild ist natürlich falschherum
. Das Ergebnis ist der Winkel der roten Gerade zur Horizontalen, also von einem Punkt rechts von P1 aus rechtsherum zur roten Gerade.
Die Funktion selbst sollte aber stimmen, habe ich aus meinem Refraction-Programm kopiert.