Eu escrevi um programa em Fortran e agora estou tentando escrevê-lo em Python. Contudo, está dando um erro:
D[i][j+1] = - u_med[i+1]/(delta_r[i]*delta_r[i+1])
IndexError: list index out of range
Sou iniciante em Python, o que pode está dando errado? Este é o código:
D = []
u_med = []
delta_r = []
for i in range(2,n):
for j in range(2,n):
if i == j:
D[i][j+1] = - u_med[i+1]/(delta_r[i]*delta_r[i+1])
t1 = u_med[i+1]/(delta_r[i]*delta_r[i])
t2 = u_med[i]/(delta_r[i]*delta_r[i])
D[i][j]= t1 + t2 + V[i]
D[i+1][j]= - u_med[i+1]/(delta_r[i]*delta_r[i+1])
t1 = u_med[2]/(delta_r[1]*delta_r[1])
t2 = 0
D[1][1]= t1 + t2 + V[1]
D[1][2]= - u_med[2]/(delta_r[1]*delta_r[2])
D[2,1]= - u_med[2]/(delta_r[2]*delta_r[1])
t1 = 0
t2 = u_med[n]/(delta_r[n]*delta_r[n])
D[n][n]= t1 + t2 + V[n]
Edit: Eu estou trabalhando com a equação de Schrödinger unidimensional num grid variável... Enfim, estou tetando escrevê-la na forma matricial, com uma matriz N x N tridiagonal simétrica.
Mesmo código no Fortran:
do i=2,n-1
do j=2,n-1
if (i.eq.j) then
D(i,j+1)=-u_med(i+1)/(delta_r(i)*delta_r(i+1))
t1 =u_med(i+1)/(delta_r(i)*delta_r(i))
t2 = u_med(i)/(delta_r(i)*delta_r(i))
D(i,j)= t1 + t2 + V(i)
D(i+1,j)=-u_med(i+1)/(delta_r(i)*delta_r(i+1))
end if
end do
end do
t1 =u_med(2)/(delta_r(1)*delta_r(1))
t2 = 0.d0
D(1,1)= t1 + t2 + V(1)
D(1,2)= -u_med(2)/(delta_r(1)*delta_r(2))
D(2,1)= -u_med(2)/(delta_r(2)*delta_r(1))
t1 =0.d0
t2 = u_med(n)/(delta_r(n)*delta_r(n))
D(n,n)= t1 + t2 + V(n)