The existence of trapped edge waves in a rotating stratified fluid with non-constant topography is studied using asymptotic and numerical techniques. A refinement of the classical WKBJ method is employed that is uniform at both the shoreline and caustic, where the classical approximation is singular, and is also uniform over long distances from the shore. This approach requires the use of comparison equations and it is shown that the two used previously in the literature are asymptotically equivalent in terms of the wave amplitude, but have small differences in the predicted wave frequencies. These asymptotic results, and results using shallow-water theory, are then compared to results from a careful numerical study of the nonlinear differential eigenvalue problem, allowing their range of practical applicability to be assessed. This numerical approach is also used to investigate whether trapping occurs in non-trivial and realistic geometries in the internal gravity wave band, which has been an open question for some time.