物理の駅 Physics station by 現役研究者

テクノロジーは共有されてこそ栄える

Pythonでエアリー関数を使う

from scipy import special
import numpy as np
import matplotlib.pyplot as plt

x = np.linspace(-15, 5, 201)
(ai, aip, bi, bip) = special.airy(x)

plt.plot(x, ai, label="ai")
plt.plot(x, bi, label="bi")
plt.ylim(-0.5, 1.0)
plt.legend()
plt.show()

f:id:onsanai:20200107200735p:plain

解説はWikipedia

ja.wikipedia.org