Fixed several errors, it works!!!!!

The sigma's needed to be index by region and not by j.  Made discret,
J, N, bc, and the sigmas a class variable. Also added functions for
finite difference and updating phi/calculating the source.