[groom03]

Hi, I've got to model how a charge will move away from a stationary charge in python, i need the acceleration to be recalculated based on the force between the charges because the moving charge will decelerate, i've been up all night trying to get this to work but i just cannot figure out whats wrong

Code:

Code:
#Problem 4a
from visual import*

#mass is assumed to be proton rest mass

k=8.987551787362e09 #coulomb constant

#particle data

#Stationary particle
alpha=sphere(pos=(0,0,0), radius=0.01,color=color.blue)
qa=0.02
m=1.67262158e-27

#Moving particle
beta=sphere(pos=(1e-03,0,0), radius=0.01, color=color.red)
#particle seperations have been magnified by 1000
qb=0.02

#particle separation
dist=beta.pos-alpha.pos

a=f/m

#Lets make it move
while beta.pos.x>15e-03: #the moving particle will have moved 15mm
    rate(20000)
    dt=0.1
    f=(k*qa*qb)/dist**2
    beta.velocity=vector(1,0,0)
    #beta.velocity=a*dt
    betadelta_s=beta.velocity
    beta.pos=beta.pos+betadelta_s
    print a
    print f
    print beta.velocity

any help greatly appreciated

p.s i know nothing about computer programming