# PHYS 1312 vpython lab 1 shell
# Fall 2022 PCS
from __future__ import division, print_function
from vpython import *

oofpez=9.0e9
qe=3e-19
sf=3e-16

source = sphere(pos = vector(0,0,0), radius = 1.e-9, color = color.red)

q_source = qe
R = 1e-8

r_obs = vector(R,0,0)
r = r_obs - source.pos
rhat = r/mag(r)