使用scipy.integrate.solve_bvp实现测地线方程的调试求助

• 学习Python的sympy和scipy库时，尝试用scipy.integrate.solve_bvp结合sympy符号计算实现测地线方程求解
• 已完成克里斯托费尔符号的符号化计算，也定义了积分步骤所需的ODE函数
• 目前无法让scipy.integrate模块正常工作：算法应收敛到连接两点且以恒定速度遍历的曲线，但现在既无法匹配初始和终点，曲线速度模长也不恒定
• 已知scipy.integrate.solve_bvp需要初始猜测值，参考了官方示例仍未解决问题，寻求调试帮助

参考资料说明：

1. 测地线方程求解相关理论
2. 椭球测地线实现的Jupyter笔记本
3. SciPy官方文档中scipy.integrate.solve_bvp的使用说明

内容的提问来源于stack exchange，提问作者Usuário 6789