y1的一阶导数等于y2； y2的一阶导数等于（-5/6）y2根号下(y22+y42); y3的一阶导数等于y4; y4的一阶导数等于-10-（5/6）y4根号下(y22+y42); 求解这四个非线性微分方程，用python来求解

可以使用 Python 中的 SciPy 库中的odeint函数来求解这个非线性微分方程组。以下是求解的代码：

from scipy.integrate import odeint
import numpy as np

# 定义微分方程组
def f(y, t):
    y1, y2, y3, y4 = y
    dydt = [y2, -5/6 * y2 * np.sqrt(y2**2 + y4**2), y4, -10 - 5/6 * y4 * np.sqrt(y2**2 + y4**2)]
    return dydt

# 初始值
y0 = [1, 0, 1, 0]

# 时间点
t = np.linspace(0, 10, 101)

# 求解微分方程组
sol = odeint(f, y0, t)

# 输出结果
print(sol)

输出结果为一个数组，其中每一行表示在对应时间点的 y1、y2、y3 和 y4 的值。