scipy.integrate.solve_bvp(fun, bc, x, y, p=None, S=None, fun_jac=None, bc_jac=None, tol=0.001, max_nodes=1000, verbose=0, bc_tol=None)# 求解ODE 系统的边值问题。 此函数在数值上求解服从 two-point 边界条件的 ODE 的一阶系统: dy / dx = f