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.