- #1
groom_03
- 1
- 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
any help greatly appreciated
p.s i know nothing about computer programming
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