Given two points $p_1,p_2$ on a circular billard table. I want to know all billard trajectories from $p_1$ to $p_2$ hitting the boundary precisely once.
Model of the circular billard:
Denote the orgin of $\mathbb{R}^2$ by $M:=(0,0)$ and assume the billard table has radius one. Thus the boundary of the table is given by $\lbrace (x,y) \in \mathbb{R}^2 : x^2 +y^2 =1 \rbrace$. One can assume that the point $p_1$ is given by $p_1=(0,a)$ for some $a \in (-1,0)$. Notice that the problem is trivial if $p_1$ is either the midpoint $M$ or is lying on the boundary (i.e. $a=0$ or $a=-1$). The situation above can be achieved by rotation.
The point $p_2$ can after rotation and reflection assumed to be given by $p_2=(b,c)$ for some $(b,c) \in (0,1) \times (-1,1)$ satisfying $b^2 + c^2 \leq 1$.
Now the problem is to find for given $a,b,c$ a point $(x,y)$ such that the billard tratjectory from $p_1=(0,a)$ to $(x,y)$ is going through $p_2=(b,c)$. (Billard trajectory means that the angle of incidence is equal to the angle of reflection).
First approach:
Denote by $B$ the point at the boundary of the billard table such that the trajectory from $p_1$ to $p_2$ goes through $B$. For given $p_1=(0,a), p_2=(b,c)$ the distance from $p_1$ to $B$ and the distance from $B$ to $p_2$ should be minimal. Therefore to find $B$ define the function
$$f: S_1((0,0)) \rightarrow \mathbb{R}; x \mapsto \Vert p_1 - x \Vert + \Vert p_2 - x \Vert$$
Thus the problem above reduces to finding the minimum of $f$. This problem can be attacked using Lagrange multipliers.
Thus one can assume $f$ is defined an all of $\mathbb{R}^2$. Its gradient is then given by
$$ \nabla f(x,y) = \begin{pmatrix} \frac{x}{\sqrt{x^2 + (y-a)^2}} + \frac{x-b}{\sqrt{(x-b)^ 2+ (y-c)^2}}\\ \frac{y-a}{\sqrt{x^2 + (y-a)^2}} + \frac{y-c}{\sqrt{(x-b)^2 + (y-c)^2}} \end{pmatrix}. $$
At this point I get stuck, since calculations seem to get to difficult.