如何让Python求解这个二阶非线性常微分方程？

import numpy as np
from scipy.integrate import solve_ivp

def equation(t, y):
    # 定义二阶非线性常微分方程
    dydt = np.zeros_like(y)
    dydt[0] = y[1]
    dydt[1] = -np.sin(y[0])
    return dydt

# 定义初始条件和时间范围
y0 = [0, 1]  # 初始条件
t_span = [0, 10]  # 时间范围

# 使用solve_ivp函数求解微分方程
sol = solve_ivp(equation, t_span, y0)

# 打印结果
print(sol.y)