ODE solver for BVProblem never calls the provided function for boundary conditions

OK, I see what’s going wrong, from your error stacktrace, the Radau solver you are trying to use is the RadauIIA5 solver from OrdinaryDiffEq.jl, RadauIIa5 in BoundaryValueDiffEq.jl should be the lowercase RadauIIa5.