Проверь коэффициенты перевода единиц:

• Масса: В геометрических единицах (G=c=1) 1 солнечная масса (M☉) соответствует ~1.477 км. Чтобы получить массу в M☉, раздели результат на этот коэффициент.

• Уравнение состояния (EOS):

  • Убедись, что функция eos(p) корректна для твоих единиц. Например, если EOS предполагает определённые коэффициенты (вроде 20 в твоем коде), проверь их соответствие единицам давления и плотности энергии.

• Начальные условия и параметры интегрирования:

  • Начальное давление 1e-13 может быть слишком малым. Проверь физически реалистичные значения для нейтронных звезд (например, ~1e35 Па в ядре).

• Проверка условий остановки:

  • Условие y[0] > 1e-17 может преждевременно останавливать интегрирование. Попробуй уменьшить этот порог или использовать условие на радиус.

• Для уравнения состояния с известным решением (например, e = 3p) убедись, что код даёт ожидаемые массу и радиус.

• Сравни результаты с литературными данными, учитывая коэффициенты перевода единиц.