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を閉じないで気軽に実行できるのがいいね。

グラフも簡単に書けるし!

ジャジャガッチ | python | 22:48 | comments(0) | trackbacks(0) |
1/1PAGES | |

03
--
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
--
>>
<<
--
PR
RECOMMEND
RECENT COMMENT
MOBILE
qrcode
OTHERS
Since 2013/09/17
LATEST ENTRY
CATEGORY
ARCHIVE
LINKS
PROFILE
SEARCH