Numerical methods for TMF cycle modeling