Форум программистов, компьютерный форум, киберфорум
Python для начинающих
Войти
Регистрация
Восстановить пароль
Блоги Сообщество Поиск Заказать работу  
 
Рейтинг 4.80/5: Рейтинг темы: голосов - 5, средняя оценка - 4.80
1 / 1 / 0
Регистрация: 26.12.2022
Сообщений: 2

SymPy неправильно решает дифференциальное уравнение

26.12.2022, 21:55. Показов 1175. Ответов 3

Студворк — интернет-сервис помощи студентам
Здравствуйте!
Решаю с помощью Google Colab задачку в SymPy. Там нужно решить дифференциальное уравнение:
y * y' = 3x - 2 , y(0) = 1
Когда решала сама на бумажке, а потом перепроверяла в калькуляторах в интернете - выходит другое решение. SymPy же решает неверно! Не могу найти, где я что ей недодаю.

Итак, ход решения на бумажке:

y'(x) = dy / dx, отсюда

y * dy / dx = 3x - 2

Дальше по правилам решения таких уравнений нужно и левую, и правую часть умножить на знаменатель, т.е. на дифференциал dx, получаем:

y * dy = (3x - 2) * dx

Это уравнение с разделяющимися переменными, его решаем, интегрируя обе части, получаем в левой части y ^ 2 / 2, а в правой - 3x ^ 2 / 2 - 2x + C, откуда вычисляем y = корень квадратный из 2 * (3x ^ 2 / 2 - 2x + C)

В SymPy пишу так:

Python
1
2
3
4
5
6
from sympy import *
import numpy as np 
x, y, t, C = symbols("x y t C")
dydx = Derivative(y, x)
mq = Eq(y * dydx , 3 * x - 2)
sol = solve(mq, dydx)[0]
Уже на этом шаге SymPy даёт мне неверный ответ: (3x - 2) / y , т.е. она не умножала на dx, а всю производную выразила через y и x.

Ну и дальше, понятно, все шаги идут уже неверными, поэтому мои письменные результаты с SymPy не совпадают.
Подскажите, пожалуйста, я вторые сутки бьюсь и зашла в тупик - что я делаю не так?
И почему не подхватывается функция dsolve ? Может, в этом корень проблемы?

Спасибо заранее за помощь!
0
Лучшие ответы (1)
IT_Exp
Эксперт
34794 / 4073 / 2104
Регистрация: 17.06.2006
Сообщений: 32,602
Блог
26.12.2022, 21:55
Ответы с готовыми решениями:

Не решает дифференциальное уравнение методом рунге кутты
Вроде бы верно написал сам метод, а решать не хочет пишет "Значение должно быть вектором" Проблема найдена

Программа решает уравнение неправильно
var a,b,x:real; begin readln(a,b,x); if a = 0 then if b = 0 then write ('любое число') else write ('нет решений') else if b...

Код неправильно решает уравнение
Помогите пожалуйста решить уравнение. Нам дано x,y,z нужно решить уравнение. Нужно чтобы сошлось с ответом w. Вот код, помогите исправить,...

3
 Аватар для Aviz__
2744 / 2053 / 507
Регистрация: 17.02.2014
Сообщений: 9,472
27.12.2022, 11:13
Лучший ответ Сообщение было отмечено salangha как решение

Решение

salangha, читай внимательно тут https://pythonru.com/biblioteki/sympy-v-python
твое выражение такое:
eqn =Eq(f(x)*f(x).diff(x), 3*x - 2)
dsolve(eqn, f(x))
4
 Аватар для Aviz__
2744 / 2053 / 507
Регистрация: 17.02.2014
Сообщений: 9,472
27.12.2022, 11:19
решение в картинке)
Миниатюры
SymPy неправильно решает дифференциальное уравнение  
3
1 / 1 / 0
Регистрация: 26.12.2022
Сообщений: 2
27.12.2022, 21:52  [ТС]
Йес! Ребятки! Всё заработало, спасибо огромное!! И отдельное спасибо за ресурс по ссылке, я до него не добралась сама.
1
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
BasicMan
Эксперт
29316 / 5623 / 2384
Регистрация: 17.02.2009
Сообщений: 30,364
Блог
27.12.2022, 21:52
Помогаю со студенческими работами здесь

Программа неправильно решает уравнение
Она решает уравнение вида x ^ 3 + ax ^ 2 + bx + c = 0 но дает непонятную ответ может исправить чтобы все давало правельно и если можно...

Программа неправильно решает квадратное уравнение
Не пойму где ошибка в проге. Работает но не считает правильно. Что не правильно сделал. ...

Написать программу которая решает дифференциальное уравнение методом Эйлера и Эйлера-Коши
Написать программу которая решает дифференциальное уравнение методом Эйлера и Эйлера-Коши y' =x+cos(y/3) a=1.6 , b=2.6 Шаг h=0.1 ...

Написать программу, которая решает уравнение с одним неизвестным и выводит в консоль значение неизвестного. Уравнение посимвольно вводится с клавиатур
Заупутался :-(

18. Написать программу, которая решает уравнение с одним неизвестным и выводит в консоль значение неизвестного. Уравнение посимвольно вводится с клави
#include <iostream> using namespace std; void main() { int urov1 = 0; int urov2 = 0; int urov3 = 0; float stad = 0; ...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
4
Ответ Создать тему
Новые блоги и статьи
Символьное дифференцирование
igorrr37 13.02.2026
/ * Программа принимает математическое выражение в виде строки и выдаёт его производную в виде строки и вычисляет значение производной при заданном х Логарифм записывается как: (x-2)log(x^2+2) -. . .
Камера Toupcam IUA500KMA
Eddy_Em 12.02.2026
Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .
И ясному Солнцу
zbw 12.02.2026
И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.
«Знание-Сила»
zbw 12.02.2026
«Знание-Сила» «Время-Деньги» «Деньги -Пуля»
SDL3 для Web (WebAssembly): Подключение Box2D v3, физика и отрисовка коллайдеров
8Observer8 12.02.2026
Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами и вызывать обработчики событий столкновения. . . .
SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL_LoadPNG (без SDL3_image)
8Observer8 11.02.2026
Содержание блога Библиотека SDL3 содержит встроенные инструменты для базовой работы с изображениями - без использования библиотеки SDL3_image. Пошагово создадим проект для загрузки изображения. . .
SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL3_image
8Observer8 10.02.2026
Содержание блога Библиотека SDL3_image содержит инструменты для расширенной работы с изображениями. Пошагово создадим проект для загрузки изображения формата PNG с альфа-каналом (с прозрачным. . .
Установка Qt-версии Lazarus IDE в Debian Trixie Xfce
volvo 10.02.2026
В общем, достали меня глюки IDE Лазаруса, собранной с использованием набора виджетов Gtk2 (конкретно: если набирать текст в редакторе и вызвать подсказку через Ctrl+Space, то после закрытия окошка. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru