Dirk2

I added the comments on which methods are used.
I also implemented a different DIRK scheme that has more stages and is L-stable, which resulted in an error that is a 100 times smaller.