Vpython Lab

This is what the original code would have looked like if we did not change the code that was given to use my Prof. Mason. It was very simple in practice. The code that made this image is below.

This code is what we made with some help from Mason during class. We changed the direction and the size of each arrow to best represent what we had.

from visual import *

## CONSTANTS

k = 9e9   ## OneOverFourPiEpsilonZero

q1 = 1.6e-19

## OBJECTS

particle = sphere(pos=vector(1e-10, 0, 0), radius = 2e-11, color=color.red)

xaxis = cylinder(pos=(-5e-10,0,0), axis=vector(10e-10,0,0),radius=.2e-11)

yaxis = cylinder(pos=(0,-5e-10,0), axis=vector(0,10e-10,0),radius=.2e-11)

zaxis = cylinder(pos=(0,0,-5e-10), axis=vector(0,0,10e-10),radius=.2e-11)

## the position of the arrow is the observation location:

Earrow1 = arrow(pos=vector(3.1e-10,-2.1e-10,0), axis = vector(1e-10,0,0), color=color.orange)

Earrow2 = arrow(pos=vector(3.1e-10,2.1e-10,0), axis = vector(1e-10,0,0), color=color.orange)

Earrow3 = arrow(pos=vector(-1.1e-10,-2.1e-10,0), axis = vector(1e-10,0,0), color=color.orange)

Earrow4 = arrow(pos=vector(-1.1e-10,2.1e-10,0), axis = vector(1e-10,0,0), color=color.orange)

Earrow5 = arrow(pos=vector(1e-10,0,3e-10), axis = vector(1e-10,0,0), color=color.orange)

Earrow6 = arrow(pos=vector(1e-10,0,-3e-10), axis = vector(1e-10,0,0), color=color.orange)

## CALCULATIONS

R1=Earrow1.pos-particle.pos

R2=Earrow2.pos-particle.pos

R3=Earrow3.pos-particle.pos

R4=Earrow4.pos-particle.pos

R5=Earrow5.pos-particle.pos

R6=Earrow6.pos-particle.pos

Ef1=((k*q1)/(mag(R1)**3))

Ef2=((k*q1)/(mag(R2)**3))

Ef3=((k*q1)/(mag(R3)**3))

Ef4=((k*q1)/(mag(R4)**3))

Ef5=((k*q1)/(mag(R5)**3))

Ef6=((k*q1)/(mag(R6)**3))

## write instructions below to tell the computer how to calculate the correct

## electric field E1 at the observation location (the position of Earrow1):

## change the axis of Earrow1 to point in the direction of the electric field at that location

## and scale it so it looks reasonable

scalefactor= 1e-20

Earrow1.axis= Ef1*R1*scalefactor

Earrow2.axis= Ef2*R2*scalefactor

Earrow3.axis= Ef3*R3*scalefactor

Earrow4.axis= Ef4*R4*scalefactor

Earrow5.axis= Ef5*R5*scalefactor

Earrow6.axis= Ef6*R6*scalefactor

print "Ef1: ", Ef1

print "Ef2: ", Ef2

print "Ef3: ", Ef3

print "Ef4: ", Ef4

print "Ef5: ", Ef5

print "Ef6: ", Ef6

## additional observation locations; do the same thing for each one

The output is such.

Our next time we use Vpython we had a cool picture both in 2D and 3D but we dont have the picture of it right now. However, i do have a prediction of what the forces around the two charges might look like in a 2D space.

 Here is the code for it. We were just representing the motion that a small positive charge would move if it was near a charge. The image is really cool. Just found the picture. It is at the bottom

from __future__ import division from visual import * ## constants k = 9e9 # stands for One Over Four Pi Epsilon-Zero qe = 1.6e-19 # proton charge s = 4e-11 # charge separation R = 3e-10 # display Enet on a circle of radius R V = 2e-10 scalefactor = 1e-20 # for scaling arrows to represent electric field ## objects ## Represent the two charges of the dipole by red and blue spheres: plus = sphere(pos=vector(s/2,0,0), radius=1e-11, color=color.red) qplus = qe # charge of positive particle neg = sphere(pos=vector(-s/2,0,0), radius=1e-11, color=color.blue) qneg = -qplus # charge of negative particle ## calculations theta = 0 while theta < 2*pi: rate(50) # tell computer to go through loop slowly ## Calculate observation location (tail of arrow) using current value of theta: Earrow1 = arrow(pos=R*vector(cos(theta),sin(theta),0), axis=vector(1e-10,0,0), color=color.orange) R1= Earrow1.pos-plus.pos R2= Earrow1.pos-neg.pos Ef1=((k*qplus)/(mag(R1)**3)) Ef2=((k*qneg)/(mag(R2)**3)) Earrow1.axis= ((Ef1*R1)+(Ef2*R2))*scalefactor print "Ef1: ", Ef1 ## write instructions below to tell the computer how to calculate the correct ## net electric field Enet at the observation location (the position of Earrow): ## change the axis of Earrow to point in the direction of the electric field at that location ## and scale it so it looks reasonable ## Assign a new value to theta theta = theta + pi/12 theta = 0 while theta < 2*pi: rate(50) # tell computer to go through loop slowly ## Calculate observation location (tail of arrow) using current value of theta: Earrow1 = arrow(pos=R*vector(cos(theta),0,sin(theta)), axis=vector(1e-10,0,0), color=color.orange) R1= Earrow1.pos-plus.pos R2= Earrow1.pos-neg.pos Ef1=((k*qplus)/(mag(R1)**3)) Ef2=((k*qneg)/(mag(R2)**3)) Earrow1.axis= ((Ef1*R1)+(Ef2*R2))*scalefactor print "Ef1: ", Ef1 ## write instructions below to tell the computer how to calculate the correct ## net electric field Enet at the observation location (the position of Earrow): ## change the axis of Earrow to point in the direction of the electric field at that location ## and scale it so it looks reasonable ## Assign a new value to theta theta = theta + pi/12 theta = 0 while theta < 2*pi: rate(50) # tell computer to go through loop slowly ## Calculate observation location (tail of arrow) using current value of theta: Earrow1 = arrow(pos=R*vector(0,cos(theta),sin(theta)), axis=vector(1e-10,0,0), color=color.orange) R1= Earrow1.pos-plus.pos R2= Earrow1.pos-neg.pos Ef1=((k*qplus)/(mag(R1)**3)) Ef2=((k*qneg)/(mag(R2)**3)) Earrow1.axis= ((Ef1*R1)+(Ef2*R2))*scalefactor print "Ef1: ", Ef1 ## write instructions below to tell the computer how to calculate the correct ## net electric field Enet at the observation location (the position of Earrow): ## change the axis of Earrow to point in the direction of the electric field at that location ## and scale it so it looks reasonable ## Assign a new value to theta theta = theta + pi/12 theta = 0 scalefactor= 1*(10**(-20.5)) while theta < 2*pi: rate(50) # tell computer to go through loop slowly ## Calculate observation location (tail of arrow) using current value of theta: Earrow1 = arrow(pos=V*vector(cos(theta),sin(theta),0), axis=vector(1e-10,0,0), color=color.red) R1= Earrow1.pos-plus.pos R2= Earrow1.pos-neg.pos Ef1=((k*qplus)/(mag(R1)**3)) Ef2=((k*qneg)/(mag(R2)**3)) Earrow1.axis= ((Ef1*R1)+(Ef2*R2))*scalefactor print "Ef1: ", Ef1 ## write instructions below to tell the computer how to calculate the correct ## net electric field Enet at the observation location (the position of Earrow): ## change the axis of Earrow to point in the direction of the electric field at that location ## and scale it so it looks reasonable ## Assign a new value to theta theta = theta + pi/12 theta = 0 while theta < 2*pi: rate(50) # tell computer to go through loop slowly ## Calculate observation location (tail of arrow) using current value of theta: Earrow1 = arrow(pos=V*vector(cos(theta),0,sin(theta)), axis=vector(1e-10,0,0), color=color.red) R1= Earrow1.pos-plus.pos R2= Earrow1.pos-neg.pos Ef1=((k*qplus)/(mag(R1)**3)) Ef2=((k*qneg)/(mag(R2)**3)) Earrow1.axis= ((Ef1*R1)+(Ef2*R2))*scalefactor print "Ef1: ", Ef1 ## write instructions below to tell the computer how to calculate the correct ## net electric field Enet at the observation location (the position of Earrow): ## change the axis of Earrow to point in the direction of the electric field at that location ## and scale it so it looks reasonable ## Assign a new value to theta theta = theta + pi/12 theta = 0 while theta < 2*pi: rate(50) # tell computer to go through loop slowly ## Calculate observation location (tail of arrow) using current value of theta: Earrow1 = arrow(pos=V*vector(0,cos(theta),sin(theta)), axis=vector(1e-10,0,0), color=color.red) R1= Earrow1.pos-plus.pos R2= Earrow1.pos-neg.pos Ef1=((k*qplus)/(mag(R1)**3)) Ef2=((k*qneg)/(mag(R2)**3)) Earrow1.axis= ((Ef1*R1)+(Ef2*R2))*scalefactor print "Ef1: ", Ef1 ## write instructions below to tell the computer how to calculate the correct ## net electric field Enet at the observation location (the position of Earrow): ## change the axis of Earrow to point in the direction of the electric field at that location ## and scale it so it looks reasonable ## Assign a new value to theta theta = theta + pi/
`````````````````````````````````````````````````````````````````````````````````