Add function to generate skew matrix of a vector
Use solution fpt jacobian in ode problems
Add ode llg jacobian function