used fitsio to read data
supressed plotaxis output on plot, but devauc plot still looks like delta function
trying to check if modelflux looks like deVauc 2D
looked over Hogg plot again still can figure out why; understand more about FITS data now, going to check  MODEL_FLUX values and whether it fits de vauc  later
commit before major change
got the 1D plot working, trying a larger sample size and seeing if I could use the r from the luv and original dev (aopt=1, aopt=2) plots
able to plot them on same plot, but had to restrict the r values
added the approximation by adding together 6,8,10 Gaussian (Mdev=6,8,10)
looking into Hogg formula for deVauc