dasslSolve now returns temporal derivatives in addition to function values
Added a support for the `maxorder` parameter and improved the test cases