临界值𝝉是怎么计算的

我看了您的《双时滞四维捕食网络的分析【基于matlab的动力学模型学习笔记_6】》，讲的真好，但是有个地方没能看懂，就是时滞的临界值𝝉是怎么计算的呀？有相应代码能参考一下吗？