Do not use the provided "intial dt" for startup steps in higher order methods

Because the startup steps are not adaptive this was causing problems if dt
was set high.

Instead always use a very small dt.
Add skew function to create skew symmetric matrices

and self tests.
Merge branch 'master' of github.com:davidshepherd7/Landau-Lifshitz-Gilbert-ODE-model