Some major refactoring of refode and refpde, possibly buggy
Some further changes to resode.
Larger class of initial data, still doesn't work for d>5
Testing the initialization values for derivatives when calling ddassl