我正在近似数值数学方法的错误,但是在编写了if else语句后,代码似乎会不断在该语句中继续,除非我将其破坏。但我想继续使用下面的代码进行仿真,以完成近似的第二部分。有什么建议吗?
while True:
m1 = h*func1(t,v, Hs)
k1 = h*func2(t,v, Hs)
m2 = h*func1(t+(h/2),v+(m1/2), Hs + (k1/2))
k2= h*func2(t+(h/2),v+(m1/2), Hs +(k1/2))
m3 = h*func1(t+(h/2),v+(m2/2),Hs+(k2/2))
k3= h*func2(t+(h/2),v+(m2/2),Hs+(k2/2))
m4 = h*func1(t+h,v+m3,Hs+k3)
k4= h*func2(t+h,v+m3,Hs+k3)
v= v +(1/6)*(m1+(2*m2)+(2*m3)+m4)
Hs=Hs + (1/6)*(k1+(2*k2)+(2*k3)+k4)
t= t+h
if Hs > 11:
f = f +0.0001
Hs=Hr*d/(L*f+d)
v = math.sqrt(2*g*Hr*d/(L*f+d))
t = 0
else:
break
while (t2 <= 10):
m12 = z*func1(t2,v2,Hs2)
k12 = z*func2(t2,v2,Hs2)
m22 = z*func1(t2+(z/2),v2+(m12/2), Hs2 + (k12/2))
k22= z*func2(t2+(z/2),v2+(m12/2), Hs2 +(k12/2))
m32 = z*func1(t2+(z/2),v2+(m22/2),Hs2+(k22/2))
k32= z*func2(t2+(z/2),v2+(m22/2),Hs2+(k22/2))
m42 = z*func1(t2+z,v2+m32,Hs2+k32)
k42= z*func2(t2+z,v2+m32,Hs2+k32)
v2= v2 +(1/6)*(m12+(2*m22)+(2*m32)+m42)
Hs2=Hs2 + (1/6)*(k12+(2*k22)+(2*k32)+k42)
t2 =t2+z
if t2 == 10:break
在Hs达到<11达一定的f后,程序将继续运行直到它变得无限大。但是,当t达到75后,我希望代码继续执行下面的while循环。
答案 0 :(得分:0)
python中有一个plot(d1, ylim=c(0, .14), main="Intersection point")
lines(d2)
points(intersection.point, col="red", pch=19)
legend("topright", pch=19, col="red", legend="intersection point")
语句。
例如:
continue
如您所见,由于>>> for i in range(10):
>>> if i == 3:
>>> continue
>>> print(i)
0
1
2
4
5
6
7
8
9
在3
为continue
时继续执行i
语句,因此缺少3
。
它并没有停止循环,只是继续运行4
。
我认为这会有所帮助。