这是我的代码。我将不胜感激任何帮助。我在“while”行遇到语法错误。 它应该做一个Verlet近似u''(theta)= u(theta)-sqrtu(theta)然后绘制u(theta)。
tlist=[]
ulist=[]
uold=uo
uo= uold-(0.5)[uold+sqrt(uold)]*dt*dt
t=0
tmax=10
dt=.001
while t < tmax
u,uold=2.0*u-uold-(u+sqrt(u))*dt*dt,u
t=t+dt
ulist.append(u)
tlist.append(t)
plot(tlist,ulist)'
答案 0 :(得分:1)
你在while语句后缺少一个冒号。
应该是:while t < tmax:
答案 1 :(得分:0)
正如@Matthias指出的那样,您在while
声明的末尾错过了:
此外,您可能希望在最后一次致电'
后失去plot
。