jacobiメソッドを使ってlaplace方程式を解くPYTHON

Question

私はPythonにはかなり慣れていて、ラプラス方程式とjacobi方程式を使って金属箱の電位を作り直そうとしています。私は最初に動作するように見えるコードを書いたが、私はエラーが出ている：IndexError：インデックス8は、軸7の軸0の範囲外であり、理由を理解できない。どんな助けも素晴らしいだろう！

from visual import* 
from visual.graph import* 
import numpy as np 

lenx = leny = 7 
delta = 2 

vtop = [-1,-.67,-.33,.00,.33,.67,1] 
vbottom = [-1,-.67,-.33,.00,.33,.67,1] 
vleft = -1 
vright = 1 

vguess= 0 

x,y = np.meshgrid(np.arange(0,lenx), np.arange(0,leny)) 

v = np.empty((lenx,leny)) 
v.fill(vguess) 

v[(leny-1):,:] = vtop 
v [:1,:] = vbottom 
v[:,(lenx-1):] = vright 
v[:,:1] = vleft 

maxit = 500 

for iteration in range (0,maxit): 
    for i in range(1,lenx): 
     for j in range(1,leny-1): 
      v[i,j] = .25*(v[i+i][j] + v[i-1][j] + v[i][j+1] + v[i][j-1]) 
      print v 

Answer 1

ちょうどあなたのコードでチラッから、それはインデックスエラーがこの部分で起こっているし、それに応じて変更することができているかのように思える：あなたは、単に追加

# you had v[i+i][j] instead if v[i+1][j] 
v[i,j] = .25*(v[i+1][j] + v[i-1][j] + v[i][j+1] + v[i][j-1]) 

と余分なiあなたのインデックスにどの範囲外だったでしょうか？