Не работает sympy.solve_univariate_inequality

@TvigMaster · Регистрация: 14.04.2023

Студворк — интернет-сервис помощи студентам

Надо поработать с пределом функции fc, но sympy выдаёт ошибку cannot be solved using
solve_univariate_inequality. Я знаю, что иногда помогает задание переменных целых чисел с точкой, но со степенями это не работает. Пробовал так-же объявить n как положительное через n = sp.Symbol('n', positive=True), но тоже не помогло. Пробовал сначала найти lim, а потом вручную вписать его в функцию, но безрезультатно.

Python
1
2
3
4
5
6
7
8
9
10
import sympy as sp
from sympy.abc import n
from sympy import oo
fc = (n**3 + n**2)**(1/3) - (n**3 - n**2)**(1/3)
lim = sp.limit(fc, n, oo)
a = sp.solve_univariate_inequality(sp.Abs(fc - lim) < 0.01, n)
b = sp.solve_univariate_inequality(n > 0, n)
c = sp.Intersection(a, b)
 
print(c)

@Royal_X · 27.10.2024, 14:10

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

Python
1
lim = sp.limit(fc, n, oo)

вот у sympy получается вычислить и он выдает 2/3

Но вот, видимо, неравенство

$https://www.cyberforum.ru/cgi-bin/latex.cgi?\left| \sqrt[3]{n^3+n^2}-\sqrt[3]{n^3-n^2}-\frac{2}{3}\right| <\frac{1}{100}$

слишком сложное для sympy, даже если попытаться решать только в вещественных числах или если попытаться решить как уравнение, а не как неравенство

Однако, можно решить как уравнение численными методами

Python
1
2
3
4
5
6
import sympy as sp
n = sp.symbols('n')
fc = (n**3 + n**2)**(1/3) - (n**3 - n**2)**(1/3)
lim = sp.limit(fc, n, sp.oo)
num_sol = sp.nsolve((n**3 + n**2)**(1/3) - (n**3 - n**2)**(1/3) - lim - 1/100, n, 1)
print(num_sol)

Выводит 3.58331494152949

Т.е. можно сказать, что решением неравенства является n > 3.58331494152949

Если же нужно точное решение, то тогда нужно использовать что-нибудь посерьёзнее. Например, такое неравенство без проблем решается аналитически с помощью Wolfram Alpha или Wolfram Mathematica:

Code
1
Reduce[Abs[-(-n^2 + n^3)^(1/3) + (n^2 + n^3)^(1/3) - 2/3] < 1/100, n, Reals] // ToRadicals

Получим

$https://www.cyberforum.ru/cgi-bin/latex.cgi?n>\frac{203}{900} \sqrt{\frac{203 \left(2^{2/3} \sqrt[3]{260829273107671+74181681 i \sqrt{20225818319}}+\frac{5146170800}{\sqrt[3]{\frac{1}{2} \left(260829273107671+74181681 i \sqrt{20225818319}\right)}}+100000\right)}{243618}}$

Численно это приблизительно равно

Code
1
n > 3.5833149415294972810278287494320883614564201591411050945639726473039788255636312608728960163392219436307797619158170749645672130922437806550633400869452242681302097623512274453514122759786851713230...

Новые блоги и статьи Все статьи Все блоги /
Отчёт о спецтехнике находящейся в ремонте Maks 20.04.2026 Отчёт из решения ниже размещен в конфигурации КА2. Задача: отобразить спецтехнику, которая на данный момент находится в ремонте. Есть нетиповой документ "Заявка на ремонт спецтехники" который. . .	Памятка для бота и "визитка" для читателей "Semantic Universe Layer (Слой семантической вселенной)" Hrethgir 19.04.2026 Сгенерировано для краткого описания по случаю сборки и компиляции скелета серверного приложения. И пусть после этого скажут, что статьи сгенерированные AI - туфта и не интересно. И это не реклама -. . .	Запрет удаления строк ТЧ документа при определенном условии Maks 19.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "Аккумуляторы", разработанного в конфигурации КА2. У данного документа есть ТЧ, в которой в зависимости от прав доступа. . .	Модель заражения группы наркоманов alhaos 17.04.2026 Условия задачи сформулированы тут Суть: - Группа наркоманов из 10 человек. - Только один инфицирован ВИЧ. - Колются одной иглой. - Колются раз в день. - Колются последовательно через. . .
Мысли в слух. Про "навсегда". kumehtar 16.04.2026 Подумалось тут, что наверное очень глупо использовать во всяких своих установках понятие "навсегда". Это очень сильное понятие, и я только начинаю понимать край его смысла, не смотря на то что давно. . .	My Business CRM MaGz GoLd 16.04.2026 Всем привет, недавно возникла потребность создать CRM, для личных нужд. Собственно программа предоставляет из себя базу данных клиентов, в которой можно фиксировать звонки, стадии сделки, а также. . .	Знаешь почему 90% людей редко бывают счастливыми? kumehtar 14.04.2026 Потому что они ждут. Ждут выходных, ждут отпуска, ждут удачного момента. . . а удачный момент так и не приходит.	Фиксация колонок в отчете СКД Maks 14.04.2026 Фиксация колонок в СКД отчета типа Таблица. Задача: зафиксировать три левых колонки в отчете. Процедура ПриКомпоновкеРезультата(ДокументРезультат, ДанныеРасшифровки, СтандартнаяОбработка) / / . . .

Опции темы

Не работает sympy.solve_univariate_inequality

Решение