 geodesic propagates from the turning point to the edge of the sphere, we need to solve how the geodesic propagates outside the sphere, and how the geodesic solution inside is matched to a solution outside. The geodesic equations outside are given by $$\begin{aligned}
\frac{f'\,\dot{r}\,\dot{t}}{f} + \ddot{t} & = 0\\
-\frac{f'}{2\,f}\dot{r}^2 + \frac{1}{2}f\,f'\,\dot{t}^2 -r\,f\,\dot{\phi}^2 + \ddot{r} & = 0\\
\frac{2\,\dot{r}\,\dot{\phi}}{r}+\ddot{\phi} & = 0 \;.\end{aligned}$$ These can be integrated to give $$\begin{aligned}
\label{eq:outgeoint1} f\,\dot{t} &=E\\
\label{eq:outgeoint2} r^2\dot{\phi} & = L\\
\label{eq:outgeoint3} \dot{r}^2 & = f\left( 1-\frac{L^2}{r^2}\right)+E^2 \;.\end{aligned}$$ To relate the constants of motion outside to the trajectory of the geodesic inside, we construct a coordinate system $(\tau,\lambda)$ which is continuous at the edge of the sphere and its close surroundings and use that to transform the geodesic derivatives in the inside coordinate system to those in the outside coordinate system. This computation is described in more detail in Appendix \[sec:continuous\], but its result is given by $$\begin{aligned}
\label{eq:matchingOS1} E &= \sqrt{\frac{\bar{a}^2}{a^2}-1}\sqrt{1+\rho_0^2} + \frac{\bar{a}}{a}\sqrt{1-\frac{\bar{\rho}^2}{\rho^2}}\dot{r}_s\\
\label{eq:matchingOS2} L &= \bar{L} = \pm \bar{a}\bar{\rho} \;.\end{aligned}$$ The procedure to compute the geodesic length for all times and angles is then as follows: Given the initial conditions $\bar{a}$ and $\bar{\rho}$, integrate equations (\[eq:OSingeoint1\] - \[eq:OSingeoint3\]) until the geodesic reaches the edge of the sphere
geodesic propagates from the turning point to the edge of the sphere, we necessitate to resolve how the geodesic propagates outside the sphere, and how the geodetic solution inside is matched to a solution outside. The geodesic equation outside are given by $ $ \begin{aligned } 
 \frac{f'\,\dot{r}\,\dot{t}}{f } + \ddot{t } & = 0\\ 
 -\frac{f'}{2\,f}\dot{r}^2 + \frac{1}{2}f\,f'\,\dot{t}^2 -r\,f\,\dot{\phi}^2 + \ddot{r } & = 0\\ 
 \frac{2\,\dot{r}\,\dot{\phi}}{r}+\ddot{\phi } & = 0 \;.\end{aligned}$$ These can be integrated to move over $ $ \begin{aligned } 
 \label{eq: outgeoint1 } f\,\dot{t } & = E\\ 
 \label{eq: outgeoint2 } r^2\dot{\phi } & = L\\ 
 \label{eq: outgeoint3 } \dot{r}^2 & = f\left (1-\frac{L^2}{r^2}\right)+E^2 \;.\end{aligned}$$ To relate the constants of apparent motion outside to the trajectory of the geodesic inside, we construct a coordinate arrangement $ (\tau,\lambda)$ which is continuous at the edge of the sphere and its close surroundings and use that to translate the geodesic derivatives in the inside coordinate system to those in the outdoor coordinate system. This computation is describe in more detail in Appendix \[sec: continuous\ ], but its consequence is given by $ $ \begin{aligned } 
 \label{eq: matchingOS1 } E & = \sqrt{\frac{\bar{a}^2}{a^2}-1}\sqrt{1+\rho_0 ^ 2 } + \frac{\bar{a}}{a}\sqrt{1-\frac{\bar{\rho}^2}{\rho^2}}\dot{r}_s\\ 
 \label{eq: matchingOS2 } L & = \bar{L } = \pm \bar{a}\bar{\rho } \;.\end{aligned}$$ The procedure to compute the geodesic length for all time and angles is then as follows: Given the initial conditions $ \bar{a}$ and $ \bar{\rho}$, integrate equation (\[eq: OSingeoint1\ ] - \[eq: OSingeoint3\ ]) until the geodesic reaches the edge of the sphere
 geldesic propagates from tme turning point to the edge kf the sohere, we need to solve how tie gwodesuc propagates outside ghe spherv, and how the teodesic solution luside ls maccied to a solutipn outside. The geodesic aqjacions outside are given by $$\begin{aligged}
\frac{g'\,\dlt{r}\,\dot{t}}{f} + \ddot{j} & = 0\\
-\fgas{f'}{2\,f}\dkn{r}^2 + \frac{1}{2}f\,f'\,\dot{t}^2 -r\,f\,\dot{\phi}^2 + \ddot{r} & = 0\\
\frac{2\,\dmt{r}\,\dot{\phi}}{r}+\ddoy{\phi} & = 0 \;.\end{aligned}$$ These cwn bf integrated to gige $$\begin{alitned}
\jqbel{eq:outgeoknt1} f\,\dot{t} &=T\\
\lcbel{eq:outgekint2} r^2\dot{\phi} & = L\\
\label{eq:outgeoing3} \dot{x}^2 & = f\left( 1-\ftcx{L^2}{r^2}\ghght)+E^2 \;.\end{almgned}$$ No relate the constands of mption outside bo thx trqjectory of the geodevic inside, we consjruct a comrbinate system $(\tau,\lambea)$ whicv is congunuuus av tge edgf oh the sphers and its coose surroundings amd lxe that to tdansfowm the geodesic derivatives in the inside cokrdinate system to thosw in the outside coorfinate syftem. This computation is described in more detail in A'pdndnw \[rwc:fontinuous\], but its result is given by $$\begin{alyfntd}
\lsbel{eq:matchingJS1} E &= \sqrt{\ftaf{\bst{a}^2}{a^2}-1}\sqrt{1+\rho_0^2} + \frxc{\bar{a}}{c}\aqdt{1-\frac{\bar{\rho}^2}{\rho^2}}\dot{g}_s\\
\label{gq:matcyingOS2} L &= \bar{K} = \pm \bar{a}\bar{\rho} \;.\end{aligned}$$ The proceduge ti compute the geodzsic length yor alk timrs and angles is then ar fomlows: Given the initjxl conditions $\baf{a}$ sng $\bar{\rho}$, integrate equatiogs (\[eq:OSinjeoinc1\] - \[eq:OSivgeoont3\]) uneil the geldesig reaches the edge lf thg sphese
geodesic propagates from the turning point to of sphere, we to solve how sphere, how the geodesic inside is matched a solution outside. The geodesic equations are given by $$\begin{aligned} \frac{f'\,\dot{r}\,\dot{t}}{f} + \ddot{t} & = 0\\ -\frac{f'}{2\,f}\dot{r}^2 + \frac{1}{2}f\,f'\,\dot{t}^2 + \ddot{r} & = 0\\ \frac{2\,\dot{r}\,\dot{\phi}}{r}+\ddot{\phi} & = 0 \;.\end{aligned}$$ These can be to $$\begin{aligned} f\,\dot{t} \label{eq:outgeoint2} r^2\dot{\phi} & = L\\ \label{eq:outgeoint3} \dot{r}^2 & = f\left( 1-\frac{L^2}{r^2}\right)+E^2 \;.\end{aligned}$$ To relate the constants motion outside to the trajectory of the geodesic we construct a coordinate $(\tau,\lambda)$ which is continuous at edge the sphere its surroundings use that to the geodesic derivatives in the inside coordinate system to those in the outside coordinate system. This computation described in in Appendix but result given by $$\begin{aligned} &= \sqrt{\frac{\bar{a}^2}{a^2}-1}\sqrt{1+\rho_0^2} + \frac{\bar{a}}{a}\sqrt{1-\frac{\bar{\rho}^2}{\rho^2}}\dot{r}_s\\ \label{eq:matchingOS2} L \pm \bar{a}\bar{\rho} \;.\end{aligned}$$ The procedure to compute the length for times and angles is then as Given the initial conditions $\bar{a}$ and $\bar{\rho}$, integrate (\[eq:OSingeoint1\] - \[eq:OSingeoint3\]) until the geodesic reaches the edge of the sphere
 geodesic propagates from the Turning poiNt to tHe eDge Of The sPherE, we need to solve HOw thE geodesic propagates outSide tHe SPherE, AnD how tHe geodeSIc SOLutIoN iNsiDe IS mAtcheD to A solutiOn outside. THe gEoDesic equatioNS oUtside are gIveN by $$\begin{aligNed}
\Frac{f'\,\dOt{R}\,\doT{T}}{f} + \ddoT{t} & = 0\\
-\fRac{f'}{2\,f}\Dot{r}^2 + \frAC{1}{2}f\,f'\,\dot{T}^2 -r\,f\,\dot{\phi}^2 + \DdOT{r} & = 0\\
\frac{2\,\DOt{r}\,\dot{\pHI}}{R}+\dDot{\pHi} & = 0 \;.\end{aligned}$$ These CAn BE integrated to gIve $$\begIn{ALiGNEd}
\lAbeL{eq:outgeoiNt1} F\,\dot{t} &=e\\
\Label{eq:OUtGEOInt2} R^2\Dot{\phi} & = L\\
\label{eQ:outgeoint3} \dOT{r}^2 & = f\Left( 1-\frAc{l^2}{r^2}\rIGht)+E^2 \;.\enD{aligNeD}$$ to rElate the conStanTs of motioN outsiDE to the tRAjectorY of the GeoDesIc inSIdE, wE coNsTRucT A cOorDInaTe system $(\TaU,\lAmbda)$ WhicH IS COntiNuoUs at The edGe of the sphere And Its cLOse SurroUndinGs anD uSe thaT to traNsforM tHe geodesic derivAtivEs in the inSidE cOorDiNate sYStem to ThoSe iN the outSide cooRDinAtE SYStEm. This computation iS dESCrIbed in moRe detaIL iN APPendix \[seC:cOntInuoUS\], But itS resULt Is given bY $$\begin{ALiGnEd}
\label{Eq:MatchiNgoS1} E &= \SqrT{\frac{\BAr{a}^2}{a^2}-1}\Sqrt{1+\rhO_0^2} + \frac{\bar{A}}{a}\sqrT{1-\Frac{\bar{\rho}^2}{\rho^2}}\dOT{r}_s\\
\label{eq:matCHiNGoS2} l &= \Bar{L} = \Pm \bAr{a}\bar{\rho} \;.\enD{aliGNed}$$ THe prOCeDurE To comPute tHe GEoDEsic length for all timEs And angLes is Then as follows: given the inITIAl conditIons $\BAr{A}$ And $\bar{\rho}$, integRate eQuations (\[eq:osingeoinT1\] - \[eq:OSIngeoint3\]) Until the gEODesic reaCheS thE edGe oF THe Sphere
 geodesic propagates fromthe turnin g poi nttoth e ed ge o f the sphere,w e ne ed to solve how the ge odesi cp ropa g at es ou tside t h es p her e,  a ndho w  t he ge ode sic sol ution insi deis  matched toa  s olution ou tsi de. The geod esi c equa ti ons outsi deare g iven b y  $$\be gin{align ed } 
\frac { f'\,\do t { r} \,\d ot{t}}{f} + \ddot { t} & = 0\\
-\frac {f'}{2 \, f }\ d o t{r }^2  + \frac{1 }{ 2}f\, f '\,\dot { t} ^ 2 -r\ , f\,\dot{\phi} ^2 + \ddot{ r } &  = 0\\ 
\ fra c {2\,\d ot{r} \, \ dot {\phi}}{r}+ \ddo t{\phi} &  = 0 \ ; .\end{a l igned}$ $ Thes e c anbe i n te gr ate dt o g i ve  $$ \ beg in{align ed }\labe l{eq : o u t geoi nt1 } f\ ,\dot {t} &=E\\
\la bel {eq: o utg eoint 2} r^ 2\do t{ \phi}  & = L \\
\l ab el{eq:outgeoint 3} \ dot{r}^2& =  f \le ft ( 1-\ f rac{L^ 2}{ r^2 }\right )+E^2 \ ; .\e nd { a l ig ned}$$ To relate t he c on stants o f moti o nou t side toth e t raje c t ory o f th e  g eodesicinside ,  w econstru ct  a coo rd ina tesyste m  $(\ tau,\l ambda)$which is continuousa t the edge of th e sp h ereand  its closesurr o undi ngsa nd  us e  that  to t ra n sf o rm the geodesic der iv atives  in t he inside coo rdinate sy s t e m to tho se i n  t h e outside coor dinat e system.T his comp utati on is de scribed i n more det ail  in  Ap pen d i x\[sec:continu o u s\],  b ut itsres ult isgiv enby$$\ be gin{align ed}
\lab el {e q: ma tch ingOS 1 } E &= \ sq rt{ \f rac {\bar { a}^2}{ a^2}- 1}\s qr t{ 1 +\r ho_0^2} +\ f rac{ \b ar {a}} {a} \s qrt{1 -\fr a c{\ bar{\rh o}^2}{\rh o^2 } }\do t{ r} _s\\
\l abel{eq:match in gOS2} L &=  \ bar {L} =\ p m \bar{a }\bar{\rho} \;.\end{ali g ned}$$The  proc edur e to comp ute  the g eod e sic le ngth f or al ltim e s  anda n gl esis  then as f o l low s: Gi ve n th e initi al conditions $\ba r {a} $ and $\bar{\ rho }$,i n te gra t ee qua ti o ns( \ [eq:OSingeoint1 \] - \[eq: OS i ng eoint3\])u nti lthe geo desic r eache s  the ed ge of the  sphere
 geodesic_propagates from_the turning point to_the edge_of_the sphere,_we_need to solve_how the geodesic_propagates outside the sphere,_and how the_geodesic_solution inside is matched to a solution outside. The geodesic equations outside are given_by_$$\begin{aligned}
\frac{f'\,\dot{r}\,\dot{t}}{f} +_\ddot{t}_&_= 0\\
-\frac{f'}{2\,f}\dot{r}^2 + \frac{1}{2}f\,f'\,\dot{t}^2 -r\,f\,\dot{\phi}^2_+ \ddot{r} & = 0\\
\frac{2\,\dot{r}\,\dot{\phi}}{r}+\ddot{\phi}_& =_0 \;.\end{aligned}$$ These can be integrated to give_$$\begin{aligned}
\label{eq:outgeoint1}_f\,\dot{t} &=E\\
\label{eq:outgeoint2} r^2\dot{\phi}_& = L\\
\label{eq:outgeoint3} \dot{r}^2 & = f\left( 1-\frac{L^2}{r^2}\right)+E^2 \;.\end{aligned}$$_To relate the constants of motion_outside to the_trajectory_of_the geodesic inside, we_construct a coordinate system $(\tau,\lambda)$ which_is continuous at the edge of_the sphere and its close surroundings and_use that to transform the geodesic_derivatives in the inside coordinate_system to_those in the outside coordinate_system. This computation_is described_in more detail_in Appendix \[sec:continuous\], but its result_is given by_$$\begin{aligned}
\label{eq:matchingOS1} E &= \sqrt{\frac{\bar{a}^2}{a^2}-1}\sqrt{1+\rho_0^2} + \frac{\bar{a}}{a}\sqrt{1-\frac{\bar{\rho}^2}{\rho^2}}\dot{r}_s\\
\label{eq:matchingOS2}_L_&= \bar{L} =_\pm_\bar{a}\bar{\rho}_\;.\end{aligned}$$ The_procedure to compute_the_geodesic length_for_all times and angles is then_as_follows: Given the initial conditions $\bar{a}$ and_$\bar{\rho}$, integrate equations (\[eq:OSingeoint1\]_-_\[eq:OSingeoint3\]) until the geodesic_reaches the edge of the_sphere