TSIRlyap function

Function to do Lyapunov exponent calculations from a TSIR simulation