We compute correlation functions, specifically 1-point and 2-point functions, in holographic boundary conformal field theory (BCFT) using geodesic approximation. The holographic model consists of a massive scalar field coupled to a Karch-Randall bran
e -- a rigid boundary in the bulk AdS space. Geodesic approximation requires the inclusion of paths reflecting off of this brane, which we show in detail. For the 1-point function, we find agreement between geodesic approximation and the harder $Delta$-exact calculation, and we give a novel derivation of boundary entropy using the result. For the 2-point function, we find a factorization phase transition and a mysterious set of anomalous boundary-localized BCFT operators. We also discuss some puzzles concerning these operators.
