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 you have read our privacy policy.

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