Нелинейный решатель ОДУ (Рунге-Кутта) ⇐ Python

[Anonymous]

Я пытаюсь решить следующую нелинейную систему (см. прилагаемую картинку), но я действительно запутался.
УРАВНЕНИЯ
Мне нужно создать один метод под названиемsolve_prec, который решает систему. Я не знаю, как использовать scipy.integrate.solve_ivp в этом случае.
Примечание: T (x=0) = 1, а область интегрирования равна (0, -0,25). Кроме того, энергетический метод относится к e, метод давления относится к p. Наконец, «ро» — это den, e0 = энергия(1,1) и p0 = давление(1,1)
Заранее спасибо!!
class eos_gamma_c():

def __init__(self, gamma=5/3.0, M0 = 10, P0 = 0.0001, k = 0.0001):
self.gamma = gamma
self.M0 = M0
self.P0 = P0
self.k = k
self.cv = 1./self.gamma/(self.gamma-1.)
return

def energy(self, tem, den):
return self.cv*tem

def pressure(self, tem, den):
return den*tem/self.gamma

def solve_prec(self, tem, den):
ene0 = self.energy(1,1)
pre0 = self.pressure(1,1)


Подробнее здесь: https://stackoverflow.com/questions/741 ... unge-kutta