1

Почему у меня не считает дифференциальное уравнение?
Это Бернулли. Вот само уравнение: xy′ − 2y − x2 * sqrt(y) = 0
Вот мой код:

from sympy import *
from scipy import *
import numpy as np

x = symbols('x')
y = symbols('y', cls=Function)
result = dsolve(Eq(y(x).diff(x), x * y(x).diff(x) - 2 * y(x) - x**2 * np.sqrt(y(x))), 0)
print('1.4.1')
pprint(result)

Вот traceback:
traceback

1 ответ 1

3

Ругается оно явно на np.sqrt, на простой корень sqrt заменить не помогает, остаётся ошибка, хоть и другая уже. Но можно для квадратного корня использовать такое написание в питоне **0.5 и так уже не ругается:

result = dsolve(Eq(y(x).diff(x), x * y(x).diff(x) - 2 * y(x) - x**2 * y(x)**0.5), 0)
6
  • к сожалению, не помогло, вот ошибка: sympy.core.function.PoleError: Cannot expand gamma(_x - 1.0) around 0
    – Iangyl
    1 дек 2020 в 20:27
  • @Iangyl Странно, у меня в Google Colab этот код прекрасно работает
    – CrazyElf
    1 дек 2020 в 21:19
  • Действительно всё работает в Colab. Я так же попробовал pow использовать, но оно же по сути тоже самое выходит. И поэтому и ошибка та же была. Я в VS Code работаю.
    – Iangyl
    2 дек 2020 в 12:30
  • @Iangyl Скорее всего в колабе другие версии библиотек, сверьте со своими версиями
    – CrazyElf
    2 дек 2020 в 12:53
  • Я обновил версии библиотек до последних в VS Code - ошибка та же.
    – Iangyl
    2 дек 2020 в 18:25

Ваш ответ

By clicking “Отправить ответ”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Всё ещё ищете ответ? Посмотрите другие вопросы с метками или задайте свой вопрос.