solve_bvp value error for coupled second order ode

Question

I am new on Python and trying to solve a BVP problem, basing myself on the code from scipy.integrate.solve_ode page. Here is what I did:

mu=0
eta