Added a function that writes the concentration at two different times, t1 and t2, to file, one file for each method. The solution methods now operate with (n-2 x n-2) matrices (boundary conditions) and the final solution is the calculated solution pluss the steady state solution.
The explicit method works, have not yet tested the implicit one. Not yet started the CN method. Changed the formalism of the explicit method to vector notation.