Построение диаграммы MR из уравнений TOV с помощью Python ⇐ Python

[Anonymous]

Я пытаюсь решить модифицированные уравнения TOV dP'(r)/dr = -1,474 ε'(r)M'(r)/r^2 (1+P'(r)/ε'(r)) (1+11,210^-6 r^3 P'(r)/M'(r))[(1-2,948M'(r)/r)]^-1 dM'(r)/dr = bo * ε '(r) где M(r)=M'(r)Mo, ε(r)=ε'(r)εο, εο= 1 МэВфм^-3, GMo/c^2=1,474 Км и 4pi/(Moc^ 2)=0,710^-40 с^2/(кгКм^2),bo=8,910^-7 км^-3 с уравнением состояния ε(P)=15,55P^0,666+76,71P^0,247. Эти модифицированные уравнения наверняка верны, поскольку они взяты из двух разных статей. Проблема, с которой я столкнулся, заключается в том, что радиус не меняется при разных значениях Pc, тогда как M меняется лишь немного. Диапазон использования Pc такой же, как в статье: Pc=1-1200 Мэв*фм^-3. Может ли кто-нибудь мне помочь? Это код, который я использую:
import numpy as np
from scipy.integrate import solve_ivp

# Constants for the modified TOV equations
epsilon_0 = 1 # MeV/fm^3
b0 = 8.9e-7 # km^-3

# Equation of state: ε(P) = 15.55 * P^0.666 + 76.71 * P^0.247
def energy_density(P):
return 15.55 * P**0.666 + 76.71 * P**0.247

# Modified TOV equations: dP/dr and dM/dr
def tov_equations(r, y):
P, M = y
epsilon_r = energy_density(P) * epsilon_0 # Energy density from the equation of state
if P

Подробнее здесь: https://stackoverflow.com/questions/790 ... ith-python