Tuesday, April 21, 2015

Vpython assignment 2

#look lines of charge!  : D  ////below is my code

from visual import *

k = 8.988e9

loc1= sphere(pos=(0,0,0), radius = 0.1, color=color.green)
loc2= sphere(pos=(1,0,0), radius = 0.1,color=color.green)
loc3= sphere(pos=(0,2,0), radius = 0.1,color=color.green)
loc4= sphere(pos=(3,0,0), radius = 0.1,color=color.green)
loc5= sphere(pos=(1,3,7), radius = 0.1,color=color.green)

#Charge A
i=0
VNetA=0
VNetB=0
VNetC=0
VNetD=0
VNetE=0

while i < 200:
   
    x=(0.1*i)

    c1 = sphere(pos=(2,-10+x,0), radius=0.2,color=color.red)

    q1= 0.02*(1e-9)      #charge of charge 2


    #loc1

    rca= mag(loc1.pos-c1.pos)

    Vca= k*q1/rca

    VNetA = VNetA + Vca


    #loc2

    rcb= mag(loc2.pos-c1.pos)

    Vcb= k*q1/rcb

    VNetB = VNetB + Vcb


    #loc3

    rcc= mag(loc3.pos-c1.pos)

    Vcc= k*q1/rcc

    VNetC = VNetC + Vcc


    #loc4

    rcd= mag(loc4.pos-c1.pos)

    Vcd= k*q1/rcd

    VNetD = VNetD + Vcd
   

    #loc5

    rce= mag(loc5.pos-c1.pos)

    Vce= k*q1/rce

    VNetE = VNetE + Vce


    i= i+1

i=0
x=0

while i < 200:
   
    x=(0.1*i)

    c2 = sphere(pos=(-2,-10+x,0), radius=0.2,color=color.blue)

    q2= -0.02*(1e-9)      #charge of charge 2


    #loc1

    rca= mag(loc1.pos-c2.pos)

    Vca= k*q2/rca

    VNetA = VNetA + Vca


    #loc2

    rcb= mag(loc2.pos-c2.pos)

    Vcb= k*q2/rcb

    VNetB = VNetB + Vcb


    #loc3

    rcc= mag(loc3.pos-c2.pos)

    Vcc= k*q2/rcc

    VNetC = VNetC + Vcc


    #loc4

    rcd= mag(loc4.pos-c2.pos)

    Vcd= k*q2/rcd

    VNetD = VNetD + Vcd
   

    #loc5

    rce= mag(loc5.pos-c2.pos)

    Vce= k*q2/rce

    VNetE = VNetE + Vce


    i= i+1



#Observation Locations


L1 = label(pos=loc1.pos,xoffset=20,yoffset=20,text= "V=%1.2f" % VNetA)
L2 = label(pos=loc2.pos,xoffset=20,yoffset=20,text= "V=%1.2f" % VNetB)
L3 = label(pos=loc3.pos,xoffset=20,yoffset=20,text= "V=%1.2f" % VNetC)
L4 = label(pos=loc4.pos,xoffset=20,yoffset=20,text= "V=%1.2f" % VNetD)
L5 = label(pos=loc5.pos,xoffset=20,yoffset=20,text= "V=%1.2f" % VNetE)




No comments:

Post a Comment