我在第24行中发现了错误“ TypeError:'int'对象不是子对象”:u [i] = un [i] -c dt / dx (un [i] -un [i-1])
这里是完整的代码:
import numpy
from matplotlib import pyplot
import time, sys
import copy
nx=41
dx=2/( nx-1 )
nt=25
dt=.025
c=1
u=numpy.ones(nx)
u[int(.5/dx):int(1/dx+1)] =2
print(u)
un = numpy.ones(nx)
for n in range(nt):
un(i)=copy.copy(n)
for i in range (1,nx):
u[i]=un[i]-c*dt/dx*(un[i]-un[i-1])
谢谢。