Moved all of the driver-type code into gamma/phi

The initial driver code which tested correctness has moved from main
to the respective files.

Additionally, the code to create Figure 1 in the paper is in gamma

Lastly, main is now a while-loop with a menu system to allow the user
to select which function they want to display or test.
ve=i
Changed "test_all" for gamma and phi so no data file is saved.
Added calculation of nesting B-spline coefficients.

This addition affected main.cpp, phi.cpp, and phi.h. Some care
has been taken in the calculation of the coefficients to avoid use
of the factorial function. Instead, a recurrence is used to calculate
the coefficients.
Added asserts to gamma&phi, started nested phi plot.

Some dynamic memory allocates were not asserted and so I've added them
where I noticed in gamma and phi code.

Main was modified by adding two new entries for Figures 4 and 5 in the
paper. These will be the nested B-splines and semi-nested C1 functions.

NOTE: There appears to be a problem with the current nested B-spline
calculation b/c the sum of the fine grid values do not equal the course
grid value.
Fixed sum of fine grids, wrote code for command file.

The fine grid values needed to be multiplied by the coefficients which
were calculated in a previous commit. This fixed the sum of the fine grids
matching the coarse grid.

I played with different gnuplot options in the command file. It seems that
dashed lines on the windows terminal is not easy?