bifurcation pragrame

经验证没有问题，和作者文章图吻合。算法是频闪法duffing——vanderpol。