[Python] Explicit Euler

아래의 ODE를 Euler method를 통해 풀어보자. 앞서 포스팅한 EE의 stability analysis에 따르면 h는 4보다 작거나 같아야 stable 하다. h가 1일때와 h가 4.2일 때로 나누어 그래프에 표시해보자. 코드는 아래와 같다. import numpy as np import matplotlib.pyplot as plt # 정확한 해를 계산하는 함수 def y(t): return np.exp(-0.5 * t) # t의 구간 정의 a = 0 b = 20 # Explicit Euler 메소드를 적용하여 수치적인 해를 계산하는 함수 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_values.append(y_next...

원문링크 : [Python] Explicit Euler