[Python] Implicit (Backward) Euler

Implicit Euler 방법으로 구한 해를 표시해보자. 문제는 Explicit Euler에서의 문제와 동일하게 하고, EE로 구한 해를 같이 표시해보겠다. import numpy as np import matplotlib.pyplot as plt from scipy.optimize import fsolve # 정확한 해를 계산하는 함수 def y(t): return np.exp(-0.5 * t) # t의 구간 정의 a = 0 b = 20 # Explicit Euler method 적용하여 수치적인 해를 계산하는 함수 def EE(h): t_values = np.arange(a, b + h, h) y_values = [1.0] # 초기 조건 y(0) = 1 for t in t_values[:-1]: y_next = y_values[-1] - 0.5 * h * y_values[-1] # y_n+1 = y_n - 0.5 * h * y_n y_values.append(y_next) re...

원문링크 : [Python] Implicit (Backward) Euler