2015

finished debugging and fixing the dense output code. I tested the polynomial interpolation code with the problem y'(t)= -sin(t) and got the correct expected results
now working on optimizing the dense output code and testing it against DOP853 and ODEX-P(12) (but haven’t really figure out yet how to characterize the error in the dense output case)
added interpolation error control and tested the whole code. Everything looks good :)