Add utility functions that help to define the system and Jacobian #46

Add:

• Print/save to a file the mass/Jacobian matrix
• Print/save to a file the vector function
• Print/save to a file the solution and initial state vectors
• Compare analytic - automatic Jacobians