Clean up
- removed mat files
- updated initpath to be the gateway file for new 
future clones
- moved fht to base level

To do:
- update the analytical solution files to be less 
confusing and cluttered.
- write a method to generate Green's function 
response for some given relative tolerance.
Improved analytical green's function calculation
Added features:
- new function endcor.m written to evaluate
hankel integral via quadrature and to compute the
end correction term for qfht
- combined analy_sol functions into a single 
function called analy_green.m
- Added error estimation to analy_green, by 
comparing qfht result with quadrature results 
generated using clenshaw-curtis quadrature.
Quadrature hankel transform + fixes
quadht.m:
- Implemented a Clenshaw-Curtis quadrature version 
of the hankel transform. Nodes and weights are 
generated using a fft method, allowing millions of
points to be simultaneously quadratures.
- The quadrature also gives an error estimate.

test_error.m:
- tests the error of the fdm against the analytical
solution, gives a relative error estimate.

Fixes:
- analy_green: moved the 1e-7 factor to the end of
the method.
- fdm_run: moved the 1e-7 factor after the fdm.

To do:
- Implement new version of the fdm
Quadrature hankel transform + fixes
quadht.m:
- Implemented a Clenshaw-Curtis quadrature version 
of the hankel transform. Nodes and weights are 
generated using a fft method, allowing millions of
points to be simultaneously quadratures.
- The quadrature also gives an error estimate.

test_error.m:
- tests the error of the fdm against the analytical
solution, gives a relative error estimate.

Fixes:
- analy_green: moved the 1e-7 factor to the end of
the method.
- fdm_run: moved the 1e-7 factor after the fdm.

To do:
- Implement new version of the fdm
Merge commit 'd3c61309341675da68d14529f01d6e477032bdac' into newfdm

Conflicts:
	green/test_error.m
Changed boundary condition merge philosophy:
- Instead of making square matrices in each domain,
rectangular domain matrices are used, where two lines
are left for the boundary condition on each side.
- Boundary condition enforcement modified to allow this
- Dirichlet condition explicitly enforced at the top
and bottom.
Implemented z-directional non-uniform grids
- Added interpolation within the coil layer
- Allowed each layer to have its own non-uniform grid
- Implemented non-uniform grid in each domain
Implemented automatic non-uniform meshing
Enabled caching in test_error
Completed new fdm code.
- Uses FDCOEFFF.M for calculation of fd coeffs.
- Solves d/dz + r d2/dr2 + d2/dz2 problem, thus 
overcoming several numerical formulation problems.
- Updated documentation.
Merge pull request #1 from ryz-mit/newfdm

Newfdm