研究は無次元化など．いつものことだがこれをやってると頭が混乱してくる．質量降着率が大きいから Msun どうしようみたいな．それと最初から無次元化してる量の扱いとか．TAのやつは Msun=1 にして軌道半径 1 AU のケプラー回転を考えれば全て上手くゆく．…

汝の父祖の遺せしものを，おのれのものとすべく，自ら獲得せよ

やっと半袖で出歩く季節になった．今年の雪駄買わないと． 本コードの前処理した．単純作業．フェルミ積分の計算がかなり重いのでもう少し粘る．xn/(1+ex) は，ζ は出てくるものの解析的計算できるのにと何度も思う．Fermi 統計も含め，統計物理の勉強しなけ…