Dann komm ich auf
A={0,-a}
Z2={e,b}
Z1=r {cos(t),sin(t)
solve({0,-a} + d (r {cos(t),sin(t)}-{0,-a}) - {e,b},{d,e})
\( \left\{ d = \frac{ao + bo}{ro \; \operatorname{sin} \left( to \right) + ao}, e = ro \; \left(ao + bo \right) \; \frac{\operatorname{cos} \left( to \right)}{ro \; \operatorname{sin} \left( to \right) + ao} \right\} \)
h(t)= ((r * (a + b)) * cos(t) / ((r * sin(t)) + a))
das ist das, was Du auch hast
h(t)= ( (a+b)rcost ) / ( a+rsint )