from __future__ import print_function, division #requested per VPython
import time
from visual import * #make sure this is asterisk char, copy/paste no go
## constants
mzofp = 1e-7
L = 2.5
q = 2*1.6e-19
scalefactor = 7e-9
deltat = 9e-20
#objects & initial conditions
particle = sphere(pos=vector(0,0,-4e-10), radius=1e-11, color=color.red)
velocity = vector(0,0,4e4)
barrow=arrow(pos=vector(-8e-11,0,0), axis=(0,0,0), color=color.cyan)
time.sleep(5) #delays the animation so I have enough time to hit record
scene.autoscale = 0 #starts the scene
while particle.y > -5e-10:
rate(10000) #ensures that only 10000 calculations are made per second
particle.pos = particle.pos + velocity*deltat
r = barrow.pos - particle.pos
rmag = sqrt(r.x**2+r.y**2+r.z**2)
rhat = r/rmag
B = mzofp*q*cross(velocity,rhat)/rmag**2
barrow.axis = B*scalefactor
well, this code runs well. I have a question - What will be the direction of the arrow representing the magnetic field? I was assuming based on the arrow position, that it's in the -y direction. Just need confirmation if either t or f.
Also another question - When the arrow representing magnetic field reaches its maximum length, what will the location of the particle be? I'm not certain how to calculate this part. Confused by this question. Any help appreciated.
