EOS

任意の状態方程式で、Kerr metric 上の平衡形状をつくりたい。その上で各タイムステップで のマップをつくらないかんということは、分かった。各タイムステップで、というところ。これはとても時間がかりそう。GRB がどうのと言えるレベルというかタイムスケールじゃないなあ。
然るにとりあえずフラクションの時間変化をとめて、計算できる形まで変形して、コードを書いた。が、全然上手くいかない。高温高密度でフラクション下がらない。むしろ一様分布、みたいな。