2017.08.22 Tuesday
pythonでオイラー法
ここんとこちょこちょこpython書いてて、少しずつ自由が効くようになってきた。まだまだだけど。
書かないと忘れちゃうのでオイラー法でも。
import numpy as np
import matplotlib.pyplot as plt
f = 1.0
N = 1000
x = 0.0
dx = 0.001
for n in range(N):
if n%100==0:
plt.plot(x,f,'.',c='blue')
#print(f)
f = (1.0+dx)*f
x += dx
X = np.arange(0,1,0.01)
Y = np.exp(X)
plt.plot(X,Y)
plt.show()
コンパイルがいらないのでvimを閉じないで気軽に実行できるのがいいね。
グラフも簡単に書けるし!