Pythonによるモンテカルロ法を用いたIsing模型の数値計算

バージョンによる変更内容を調べたいのであればリリースノートを確認すればよいかと思います。
参考にされた記事では Python3.5.1 を使用していると書かれておりますので、まずはお手元でも同じバージョンで実行してみてはいかがでしょうか？

また、参考にされた記事では T=2.2691 の値は式で算出されております。
この分野には疎いので私の理解が誤っている可能性はございますが、図からはあくまでそのような傾向が読み取れるということで、正確にT=2.7であるということが読み取れますという内容では無いように思いました。