Added function for source term.
Problem with algoritm.

After testing solving algorithm and plotting using the python script I
have found that there must be a mistake somewhere. The result that the
Thomas algorithm now produces looks mistakingly much like the source
term as described in ClosedForm.py.

Need to look over the algorithm and check that it does what it's
supposed to do.