【发布时间】:2015-03-01 07:58:15
【问题描述】:
我的代码编写如下,但由于索引超出范围而无法运行,有人可以告诉我哪里做错了吗?
def potential(a,b,c,d):
energylist = []
for x in np.arange(-5,10,0.1):
for y in np.arange(-5,10,0.1):
expa = np.exp(-(x-a)**2)
expb = np.exp(-(y-b)**2)
expc = np.exp(-(x-c)**2)
expd = np.exp(-(y-d)**2)
if x-y != 0:
integ = expa*expb*(1/np.abs(x-y))*expc*expd
energylist.append(integ)
val = sum(energylist)
return val
def htable(*myions):
k = len(myions)
matrix = np.zeros((k^2,k^2),float)
for i in np.arange(0,k^2-1,1):
for j in np.arange(0,k^2-1,1):
m = myions[i/k]
n = myions[i%k]
o = myions[i/k]
p = myions[i%k]
matrix[i,j] = potential(m,n,o,p)
print matrix
return matrix
htable(1,3)
【问题讨论】: