#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