2016-04-02 9 views

答えて

0

:しかし、この部分は次のようなものになるだろう。

ynew = [dt * u[i] + y[i] for i in xrange(1,n+1)] 
unew = [dt * (k/m) * (y[i+1] + y[i-1]) for i in xrange(1,n+1)] 
- ここ

はPython的な方法です

0

これは大きなプログラムの一部にすぎないと思われますので、私はあなたが表示していないものを記入することはできません。私は境界エラーのうちのインデックスを取得しないように、私はuynbeforeを使用し、nafterのインデックス作成には注意でしょう

for n in range (1, N + 1): 
    nbefore = n - 1 
    nafter = n + 1 
    ynew[n] = dt * u[n] + y[n] 
    unew[n] = dt * (k/m) * (y[nafter] + y[nbefore]) 
関連する問題