在做碳税模拟，怎么都做不出来了，希望可以找一篇有详细代码的文章参考。