Посчитать дифференциал

@Iangyl · Регистрация: 31.03.2018

Здравствуйте!

Подскажите, пожалуйста, почему не считает дифференциал?

Вот код:

Python

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))), y(x))
print('1.4.1')
pprint(result)

Вот traceback:

@dondublon · 02.12.2020, 11:42

Похоже, вы смешали scipy и sympy. Не надо, это разные вещи. scipy работает с настоящими числами, sympy - со своими символьными объектами.

@Iangyl · 02.12.2020, 12:25 **[ТС]**

dondublon, там нигде не используются методы scipy, потому что таковых в нём нет.

Добавлено через 37 секунд

Сообщение от Iangyl

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))), y(x))

это всё sympy

@dondublon · 02.12.2020, 12:30

А мне показалось, что dsolve оттуда. Не видно, когда import *, да ещё и из двух модулей.
В любом случае, тот ufunc, про который говорится в тексте ошибки - из numpy.

@Iangyl · 02.12.2020, 15:26 **[ТС]**

dondublon, я сделал по другому, вот-так:

Python

1	result = dsolve(Eq(y(x).diff(x), x * y(x).diff(x) - 2 * y(x) - x*2 y(x)**0.5), 0)

В VS Code мне это принесло такую ошибку: sympy.core.function.PoleError: Cannot expand gamma(_x - 1.0) around 0

А в Google Colab, всё заработало:

@dondublon · 02.12.2020, 16:16

Iangyl, интересный случай. Браво, не знал, что так можно.

@dondublon 4387 / 1825 / 338 Регистрация: 17.03.2012 Сообщений: 9,505 Записей в блоге: 5
	02.12.2020, 11:42	2
	Похоже, вы смешали scipy и sympy. Не надо, это разные вещи. scipy работает с настоящими числами, sympy - со своими символьными объектами. 0

@Iangyl 23 / 27 / 18 Регистрация: 31.03.2018 Сообщений: 426
	02.12.2020, 12:25 [ТС]	3
	dondublon, там нигде не используются методы scipy, потому что таковых в нём нет. Добавлено через 37 секунд Сообщение от Iangyl 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))), y(x)) это всё sympy 0

@dondublon 4387 / 1825 / 338 Регистрация: 17.03.2012 Сообщений: 9,505 Записей в блоге: 5
	02.12.2020, 12:30	4
	А мне показалось, что dsolve оттуда. Не видно, когда import *, да ещё и из двух модулей. В любом случае, тот ufunc, про который говорится в тексте ошибки - из numpy. 0

@dondublon 4387 / 1825 / 338 Регистрация: 17.03.2012 Сообщений: 9,505 Записей в блоге: 5
	02.12.2020, 16:16	6
	Iangyl, интересный случай. Браво, не знал, что так можно. 0

Опции темы