atcoder の D - Remainder Reminder のような，剰余 a%b を調べる問題で a と b のどちらを固定するかで議論の難しさがとても違ったので（自分は最初 a を固定して考えていて分からなくなった），この際，剰余の関数グラフをプロットして見てみる． 特別な結…

無限体積の2次元イジングモデルの厳密解であるオンサーガー解を数値的に計算する． 1スピン当たりの自由エネルギー \begin{align} -\beta f &= \log(2\cosh 2K) + \frac{1}{2\pi ^ 2} \int ^ \pi _ 0 d\omega _ 1 \int ^ \pi _ 0 d\omega _ 2 \log(1 - k ^ 2…

numpy.gradientで配列を微分するとき，配列の境界部分（始端終端）の値だけ目的の正しい微分値が出ないことがあった． すぐに解決したのでメモです． 境界部分での微分はedge_orderというオプションで精度を選べる． numpy.gradient — NumPy v1.16 Manual 何…

python3.6のチュートリアルでループ制御でのelse節について知ったのでメモです． docs.python.jp ループのelse節の動作についてはチュートリアルに書いてある通りでたぶん理解できたと思います． チュートリアルにあったサンプルプログラムについて考えてみ…