q is a function of u

I changed q to be a function of u, rather than to be a set function of
phi_0.  I thought this made more sense since we don't have incident
light, we have a dipole oriented at some angle.
changed to pettersson approach

both approaches are in the file now.  Approach 1 (currently in use) is
from petterson, approach 2 is from neyts