[Python] 2nd ODE with numerical method(1)

이번에는 2차 미분 방정식을 파이썬을 통해 나타내보자. 주어진 식은 아래와 같다. 식을 풀어보면 y=cos(ωt) 라는 것을 알 수 있지만, 식을 풀지 않고 해에 대한 그래프를 그려보자. import numpy as np from scipy.integrate import odeint import matplotlib.pyplot as plt # 주어진 미분 방정식 # y[1]=y', y[0]=y # dydt = [y', y''], y''=-omega**2 * y def model(y, t): dydt = [y[1], -omega**2 * y[0]] return dydt # 초기 조건 y0 = 1.0 # 초기 위치 v0 = 0.0 # 초기 속도 omega = 4.0 # 감쇠율 # 시간 범위 및 간격. exact 해를 찾는 거니까 촘촘하게 t = np.linspace(0, 6, 1000) # odeint를 사용하여 미분 방정식을 푸는 함수 호출하여 solution 변수에 입력 # odei...

원문링크 : [Python] 2nd ODE with numerical method(1)