Added a main file for unit tests.
Added solve_ndiag function.  Updated unit tests.
Fix loop terminating condition error.  Was running one extra time before.
Add first draft of redone MB3 acceleration. Problems converging when dx is small.  Need to study/fix.
Fixed Sth = sigma_t / h typo.  Should be sig_t * h.  Fix two for loop bounds for Qhat.  Answer now symmetric.
Fix three bugs with accelerate_MB3().  Still converges to a different solution than nonaccelerated MB3.

Three bugs:
1.  zeta_1, zeta_3 not initialized before using.
2.  L_tilde on edges was wrong.
3.  Sign errors in first and last rows of A.
Improved definition of L-tilde on edges.  A matrix modified for new def.  Fix one bug with L-tilde.

The one bug was that I was accidentally multiplying zeta_1 by 2 in the definition of L_tilde in the non-boundary edges.

There is still a slight discrepancy causing the accelerated scheme to converge to a slightly different solution.
Added Utilities function print_dmatrix_matlab to print out matrices for easy input to Matlab.
Acceleration now working after bug fixes.  Spectral radius seems like it matches Fourier predictions.

L_tilde in non-boundary edges misdefined -- this is fixed now.
Boundary cases for Qhat acceleration were redone.
Removed some comments and printouts.  Added # in front of printouts for better C&P compatibility with Python.
Added a new main file from 20141117
Add main file for epsilon problem for MC paper.
Minor style bug in print_dmatrix_matlab.
Made MB-3 acceleration functional for non-homogeneous case.

1.  Replaced material variables with edge material variables in the auxiliary equation.
I did so by changing the leftIteration and rightIteration formulas and by making an
exception for the def. of source_edge
2.  Fixed for-loop bounds for _avg variables.
3.  Fixed evaluations of material properties in acceleration method.