* add DtypeError with informative error message

* use float64 in augment system because of scipy

* set dtypes explicitly in all places

* fix incorrect sensitivity matrix column order
closes #3676
+ order of sensitivity columns now beginning with y0
+ tests were aligned
+ order of y0, theta now everywhere

* more detailed docstrings

* re-run after the fix

* update & re-run after the fix