TSIRllyap function

Function to calculate the local Lyapunov exponents for the TSIR