当计数器 < 时间时：w = d*u.copy()对于范围（1，节点 -1）内的 i：对于 np.arange（-1/2, 1/2）内的 j：u[i] = ((w[j+1//2]/dx**2)*(u[i+1]-u[i]))counter += dtif f <= 0.5：IndexErr...

当计数器<时间时：

w = d*u.copy()

for i in range(1, nodes -1):
    for j in np.arange(-1/2, 1/2):

        u[i] = ((w[j+1//2]/dx**2)*(u[i+1]-u[i]))

counter += dt
if f <= 0.5 :

IndexError：只有整数、切片（ : ）、省略号（ ... ）、numpy.newaxis（ None ）和整数或布尔数组才是有效索引

代码程序