Нужна помощь в решении ОДУ второго порядка с заданными условиями

Всем привет! Помогите, пожалуйста, решить следующее ОДУ:

Ниже привела свой код, но графики получаются, мягко говоря, не красивые. Догадываюсь, что есть ошибка. Не понимаю, куда отправить начальное значение pi/4 и что отправить в np.linspace. Заранее спасибо за помощь!

from scipy import integrate
import numpy as np
from matplotlib import pyplot as plt


def f(y, t):
 v = y[0]
 x = y[1]

 f0 = 1/(np.sin(x)**3)-v
 f1 = v
 return [f0, f1]

# массив точек интегрирования
ti = np.linspace(np.pi/4,np.pi/2, 1000)
# начальная координата
x0 = 2.0
# начальная скорость
v0 = 1.0
# список начальных условий
y0 = [x0, v0]
# решение ОДУ
sol = integrate.odeint(f, y0, ti)
# рисунок для построения графиков функций x(t) и v(t)
fig = plt.figure()
# панель для рисования графика координаты
ax1 = fig.add_subplot(121)
# панель для рисования графика скорости
ax2 = fig.add_subplot(122)
ax1.set_xlabel(r'$t$')
ax1.set_ylabel(r'$x$')
ax1.set_title('координата')
# график x(t) - зависимость нулевого столбца вектора
# решения sol от ti
ax1.plot(ti, sol[:, 0])
ax2.set_xlabel(r'$t$')
ax2.set_ylabel(r'$v$')
ax2.set_title('скорость')
ax2.plot(ti, sol[:, 1])
# настройка оптимального расположения панелей
plt.tight_layout()
fig.savefig("ode2.png")