# PHYS 1312 vpython lab 1 shell # Fall 2016 PCS from __future__ import division, print_function from visual 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)