from __future__ import print_function, division
from visual import *
scene.width = scene.height = 1000
max = 15*0.5
scene.range = max
scene.autoscale = 0
scene.background = color.white

proton = sphere(pos=(0,0,0),radius=0.2, color=color.red)
rad = 0.5
arr = []
for theta in arange(0,2*pi,pi/8):
    for n in arange(1,10,1):
        obs = vector(n*rad*cos(theta), n*rad*sin(theta),0)
        rhat = norm(obs - proton.pos)
        aa = arrow(pos=obs, axis=rad*rhat, color=(1,.5,0))
        arr.append(aa)
        
scene.mouse.getclick()
for theta in arange(0,2*pi,pi/8):
    cc = curve(pos=[proton.pos,(25*rad*cos(theta),25*rad*sin(theta),0)],
               color=(1,.5,0), radius=0.03)
for aa in arr:
    aa.visible = 0