Python ода первого порядка, как решить эту проблему с помощью Sympy

Question

Когда я пытаюсь решить эту первую оду, используя Sympy, как показано ниже:

import sympy
y = sympy.Function('y')
t = sympy.Symbol('t')

ode = sympy.Eq(y(t).diff(t),(1/y(t))*sympy.sin(t))
sol = sympy.dsolve(ode,y(t))
csol=sol.subs([(t,0),(y(0),-4)]) # the I.C. is y(0) = 1
ode_sol= sol.subs([(csol.rhs,csol.lhs)])
print(sympy.pprint(ode_sol))

Это дает мне эту ошибку:

Traceback (most recent call last):
File "C:/Users/Mohammed Alotaibi/AppData/Local/Programs/Python/Python35/ODE2.py", line 26, in <module>
csol=sol.subs([(t,0),(y(0),-4)]) # the I.C. is y(0) = 1
AttributeError: 'list' object has no attribute 'subs'

Answer 1

Ваша проблема в том, что у этого ODE нет уникального решения. Таким образом, он возвращает список решений, которые вы можете узнать из сообщения об ошибке и распечатав sol,

Сделать оценку в цикле,

for psol in sol:
    csol = psol.subs([(t,0),(y(0),-4)]);
    ode_sol = psol.subs([(csol.rhs,csol.lhs)]);
    print(sympy.pprint(ode_sol))

чтобы найти следующую ошибку, эта замена не решает для константы. Что работает, чтобы определить C1=sympy.Symbol("C1") и используя

    ode_sol= psol.subs([(C1, sympy.solve(csol)[0])]);

но это все еще кажется счастливым. Или лучше избегать сообщений об ошибках для неразрешимости второго случая:

C1=sympy.Symbol("C1");
for psol in sol:
    csol = psol.subs([(t,0),(y(0),-4)]);
    for cc1 in sympy.solve(csol):
        ode_sol= psol.subs([(C1, cc1)]);
        print(sympy.pprint(ode_sol))