What could be wrong with my charge motion model in Python?

  • Context: Python 
  • Thread starter Thread starter groom_03
  • Start date Start date
  • Tags Tags
    Work
Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
Sign up Log in
1 reply · 2K views
groom_03
Messages
1
Reaction score
0
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 what's wrong

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
 
Physics news on Phys.org
What are some of the error messages that you're getting?
Also, use print statements at each stage to check that the values seem sane.

Also, far as I can tell you don't define f for
Code: 
a=f/m