主要是在出现不稳定极限环时应该如何处理。
不稳定极限环在时间演化曲线中自然是不能直接得到的，这里如果想找到他就要利用不稳定极限环作为态空间分界面的性质。在极限环内部，最终趋于稳定点，例如（0，0），而极限环外，很有可能就发散了，或者到了远处的某稳定点或者某稳定极限环。
于是，画分岔图非常不爽地需要两重循环，第一重是初值的循环，利用的是二分法，逐渐逼近分界面；第二重自然是参数。
在做数值计算的时候，很重要的一个问题是，通常在有限的时间内（我们自然不可能算到无穷久以后，时间的设定可以先跑几次看一看，大致设一个比较大的，能保证分得开极限环两侧性质的也就好了～其实如果初值在极限环外，很快值就会增得很大，以至于Matlab不干了，给你抱怨两句，你丫给我这么大的数，还让我算这么久，就算到这儿了，后面的不算了，然后自动跳入下一循环～足以看出Matlab还是一个贤妻受啊，抱怨归抱怨，没撂挑子不是，反正算到那儿也的确足够了，比起极其挑剔的C或C++来，好对付多了，慢一点也可以原谅了），初值在极限环内时，最终并非严格等于零的，这是需要设定一个比较小的数作为标准（不能太小）。
不罗嗦了，上代码。(some method to sovle equation of Lorentz)