Sympy의 ODE 솔버 초기 조건은

Question

내가 초기 값 문제지만 초기 조건 ics를 해결하기 위해 dsolve를 사용하려면

아무런 영향이없는 것으로 나타 영향을주지 않습니다. 이 예에서 dsolve은 동일한 결과를 제공합니다. 두 경우 모두

from sympy import * 
x = symbols('x') 
f = Function('f') 

eq = Eq(Derivative(f(x), x), f(x)) 

# No initial condition 
pprint(dsolve(eq, f(x))) 

# With initial condition f(0) = 1 
pprint(dsolve(eq, f(x), ics={f(0):1})) 

나는

  x 
f(x) = C1 ℯ 

C1

도 ics으로 값 1로 대체되지 않습니다 동일한 솔루션을 얻을. Sympy second order ode에 sympy issue 4720이 언급되어 있지만 이제는 해당 문제가 종결되었습니다. SymPy 1.1.1을 사용하고 있습니다.

Answer 1

PR implementing the initial conditions는 SymPy 1.1.1의 일부는 아니지만 1.2이어야한다. 당신은 SymPy의 버전이 (가에서 오류가 발생하는 경우, 당신은 ICS 전원 시리즈를 제외하고 지원이없는)

from sympy.solvers.ode import solve_ics 

을 실행하여이 있는지 여부를 테스트 할 수 있습니다. 당신이 언급 한 후뿐만 아니라, 내 answer 찾아 정수를 대체 할 (다른 상황에서 서로 다른 유형을 반환하는 습관이있다) solve를 사용하는 방법에 대한 간단한 설명이 있습니다. 요청을 당겨

eq = Eq(Derivative(f(x), x), f(x)) 
sol = dsolve(eq, f(x)).rhs 
constants = solve([sol.subs(x,0)-1], dict=True) 
print(sol.subs(constants[0]))