Решить диффур питон

from sympy import *

Никогда так не делай в рабочем коде.

я учусь

Для учёбы можно.

Задача именно диффур решить или сделать это на питоне?

Если первое, то Maxima однозначно твой выбор.

на питоне

Сначала написать дифур на бумашке в математических терминах и показать людям

https://ibb.co/n0CqRTm

Вроде правильно

Наверное надо взять что то вместо dsolve

Говорят pdsolve но он не решает, может из-за того что я не ввел 0=<x=<1

С этими cas сплошная молока. Да, им все надо разжевывать

Спросил у GPT бота.

Вот его ответ, я не знаю насколько он прав или сочиняет:

Начало ответа бота"

Ошибка возникает из-за того, что dsolve из библиотеки SymPy не может решать уравнения, содержащие функции нескольких переменных, в данном случае u(x,t).

Чтобы решить эту проблему, можно заменить функцию u(x, t) одной переменной u(z), где z - это новая переменная, зависящая от x и t. Это можно сделать с помощью замены z = x + c*t, где c - произвольная константа. Тогда x и t могут быть выражены через z как x = z - c*t и t = (x - z) / (-c).

Применяя эту замену, можно переписать уравнение в следующей форме:

u_z = diff(u(z), z)
u_zz = diff(u_z, z)
eq = u_z + c*u_zz - 1 - 2*(z - c*t)*t - 5*t*sin(3*pi*(z - c*t)/2)

Здесь u_z и u_zz - это первая и вторая производные функции u(z), соответственно. eq - это новое уравнение с замененной функцией u(z).

Также нужно изменить начальные условия в соответствии с новой переменной z:

con = {u(0): t, u_z.subs(z, 1): t*t, u(z).subs(z, 0): 3*z}

Здесь u_z.subs(z, 1) заменяет x на 1 и t на (1 - z) / (-c), а u(z).subs(z, 0) заменяет x на 0 и t на z / c.

Измененный код будет выглядеть следующим образом:

from sympy import *

if __name__ == '__main__':
    z, c = symbols('z c')
    u = Function('u')
    u_z = diff(u(z), z)
    u_zz = diff(u_z, z)
    eq = u_z + c*u_zz - 1 - 2*(z - c*t)*t - 5*t*sin(3*pi*(z - c*t)/2)
    con = {u(0): t, u_z.subs(z, 1): t*t, u(z).subs(z, 0): 3*z}
    y = dsolve(eq, ics=con)
    print(eq)
    print(simplify(y))

Здесь c является произвольной константой замены, которую нужно выбрать самостоятельно.

"конец ответа бота.

А почему не sage? Он же на питоне.

Ну во-первых Maxima привычнее.

Во-вторых мне Lisp кажется более интересным.

в натуре четко, хотя он похоже 9 потерял у второй производной и знак поменял

Да, за ботом нужно перепроверять.

Я попросил уточнить бота:

Мой вопрос: «Не потерял ли ты 9 у второй производной? И не ошибся ли со знаком, можешь ли перепроверить?»

Его ответ:

Вы правы, при замене переменных я допустил ошибку и потерял коэффициент 9 перед второй производной. Исправленный уравнение должно выглядеть следующим образом:

u_z = diff(u(z), z)
u_zz = diff(u_z, z)
eq = diff(u(z), t) - 9*u_zz - 1 - 2*z*t - 5*t*sin(3*pi*z/2)

Что касается знака перед sin в уравнении, я также допустил ошибку и указал неправильный знак. Верный знак должен быть плюс, то есть уравнение должно выглядеть так:

u_z = diff(u(z), z)
u_zz = diff(u_z, z)
eq = diff(u(z), t) - 9*u_zz - 1 - 2*z*t + 5*t*sin(3*pi*z/2)

Извините за путаницу. Спасибо, что обратили на это внимание!

теперь замену переменных потерял

Новый ответ бота:

Да, извините, я допустил ошибку в предыдущем ответе. Действительно, в уравнении необходимо использовать замену переменных, чтобы привести его к уравнению вида u_zz + p(z)*u_z + q(z)*u = f(z, t).

Исправленный код должен выглядеть следующим образом:

from sympy import *

if __name__ == '__main__':
    x, t = symbols('x t')
    u = Function('u')
    z = x + c*t
    u_z = diff(u(z), z)
    u_zz = diff(u_z, z)
    p = simplify(diff(z, x))
    q = simplify(diff(z, x, x))
    f = -1 - 2*x*t - 5*t*sin(3*pi*x/2)
    eq = u_zz + p*u_z + q*u - f
    con = {u(0): t, u_z.subs(z, 1): t*t, u(z).subs(z, 0): 3*x}
    y = dsolve(eq, ics=con)
    print(eq)
    print(simplify(y))

Здесь p и q - это функции, соответствующие коэффициентам уравнения при первой и второй производных функции u(z). Переменная f - это правая часть уравнения.

Извините за путаницу, и благодарю за замечание!

А почему LISP, а не какой-нибудь Haskell?

Ну записывай недостатки хаскеля:
1) не гомоиконный
2) не поддерживает императивное программирование

1 не понял
2 А как же монады? В том числе ST

не вдаваясь в код скажу лишь что в нем опять ошибки - он не запускается

Код отформатируй, а то прочитать сложно

В VSCode ctrl+shift+i

Почти во всех туториалах для новичков так почем то делают

Ну всё, наступила эпоха, и ЛОР не нужен теперь…