Added new 'system_jacobian' function for use by IDA (maybe elsewhere?) Refactored the matrix output stuff in IDA. Fixed the index checking in idaanalyse Still a problem with checking rank of small matrices.