[Python] ODE_Trapezoidal method

TR을 이용하여 미분방정식을 풀어보자. 주어진 식에 대해 몇 가지 타입 스텝별로 해의 거동을 살펴보겠다. import numpy as np import matplotlib.pyplot as plt from scipy.integrate import odeint from scipy.optimize import newton # 미분 방정식 정의 def f(y, t): dydt = -((3 * np.tanh(3 * t) + 5) * y) return dydt # Trapezoidal Method로 한 단계 전진을 위한 함수 def trapezoidal_step(y_n, t_n, dt): # newton을 사용하여 비선형 방정식 풀이 implicit_function = lambda y_next: y_n + 0.5 * dt * (f(y_n, t_n) + f(y_next, t_n + dt)) - y_next y_next = newton(implicit_function, y_n) return y_...

원문링크 : [Python] ODE_Trapezoidal method