2015

added the n-body IVP to the test cases
tested parallel vs serial on the n-body problem
noticed that we’re slowed down by concatenating of ndarrays when constructing ts and ys arrays, and figured out later that we don’t need keep intermediate values of the computed solution, and only need y(tf). I changed the code across all the files to fix this.