python - スプリングに逆らう物理を遅くする

Question

私は物理クラスのチュートリアルを行っています。

バネで動かされる箱の上にカートがあり、停止するまで減速するプログラムを作成しましたが、実行すると、バネがカートを加速しているように見えます (どちらを否定しても)。

vPython の数値の丸めの問題で加速が発生する可能性があると聞きましたが、これが本当なら、すべての数値を 1000 倍大きくして修正できますか?

ありがとう！

from visual import *
from visual.graph import *

length=1.0

track=box(pos=vector(0,-0.05,0),
          size=(length, 0.05, 0.10),
          material=materials.bricks,)

# creates the track which is "length"meters in the
# x direction, 0.05m tall, and 1m deep
start=-0.5*length+0.05

cart=box(pos=vector(start+0.01,0,0),
         size=(0.1,0.05,0.1),
         color=color.green)
k=-4
#spring constant

sprL=(start-0.05)-0.1

#sets position of left end of spring
spring=helix(pos=(sprL,0,0),
             axis=((cart.x-0.05)-sprL,0,0),
             radius=0.02,
             color=color.yellow)

cart.m=0.70
#mass of cart
cart.vel=vector(0,0,0)
#initial velocity of cart
cart.force = k*(cart.x)*vector(1,0,0) 
#force of the spring
cart.accel=cart.force/cart.m
#acceleration of the cart taking into account the fan
t=0
deltat=0.01
end=0.5*length-0.05
#defining the end of the track
gdisplay(x=100,
         y=500,
         xtitle='time (sec)',
         ytitle='X (cyan), Px (red)')

xcurve = gcurve(color=color.cyan)
pcurve= gcurve (color=color.red)

while cart.x<end+0.01 and (cart.x>(start-0.01)):
    #we include -0.01 so the cart does not fail immediately upon start...
    cart.pos = cart.pos + cart.vel*deltat+(0.5)*(cart.accel)*deltat**2
    #x equals x naught plus v times delta t plus one half a delta t squared
    #note that ** means "to the power of"
    xcurve.plot(pos=(t,cart.x))
    pcurve.plot(pos=(t,cart.vel.x))
    cart.vel=cart.vel+cart.accel*deltat
    #new velocity is old velocity plus acceleration times time
    cart.force=k*(cart.x)*vector(1,0,0)
    cart.accel=cart.force/cart.m
    spring.axis=((cart.x-0.05)-sprL,0,0)    
    t=t+deltat
    #increments time
    rate(100)
    #rate means no more than 100 loops per second