python - Python - 钟摆的 Runge-Kutta 方法 - 怎么做 o_O

Question

我是 python 新手，我想在这段代码的末尾添加 runge-kutta 方法：

GlowScript 2.9 VPython

a0 = float(input('initial swing [A0]: '))
b0 = float(input('initial speed [B0]: '))
g0 = float(input('Earth acceleration [G0]: ')) 
cycle = 500
t = 0
l = 10
dt = 0.1

ball=sphere(pos=vector(10,2,0),radius=0.5,color=color.blue)
pivot=vector(0,20,0)
rod=cylinder(pos=pivot,axis=ball.pos-pivot,radius=0.1,color=color.red)
 


while (t<100):
 rate(100) 
 acc=-g0/l*sin(a0) 
 a0=a0+b0*dt 
 b0=b0+acc*dt 
 ball.pos=vector(l*sin(a0),pivot.y-l*cos(a0),0)
 rod.axis=ball.pos-rod.pos 
 t=t+dt

你能帮我吗？