6:00に起きて妻が購入を予定している実験器具の相談に乗ってから，子どもたちを起こして一緒に食事をし，正月休みの間に溜まった可燃ゴミ２袋を出してから家を発った。相変わらず小雪が舞っていた。長野電鉄で長野駅へ出て，あさま504号には７分間の乗り換え時間があり，ちょうど良かった。この行き方は楽だし時間もかからないが，やはり雪でなければ自転車で行きたいと思う。
今日は，数値演算ライブラリを探索する予定である。シミュレーションなどのプログラムをＣで書くとき，既に確立した数値演算ルーチンについては，車輪の再発明をするのも嫌なので，ライブラリを利用したいと思うわけだが，フリーで安心できるものはなかなか見つからない。後で役に立つだろうと思い，検索結果をまとめておくことにした。今回探しているのは，尤度関数が微分不可能なときの最尤法のルーチン，あるいは２変数関数の最小値を与えるパラメータを求める滑降シンプレックス法かPowellの方法のルーチンでもいいのだが，実は今月届くはずの人口学研究に掲載される結果では試行錯誤して得た尤度曲面を目で見て求めたパラメータを，自動的に求めるプログラムを作って公開するためのものである。公開に当たっては，データの形式や使い方も書かねばならないが，小集団での出産間隔データからの出生力の分析には役に立つと思われるので，作る価値はあるだろう。ついでにカプラン＝マイヤ法もやらせれば（グリーンウッドの公式で標準誤差も求める），もっと役に立つだろう。学生に実習させるため，15日の講義までには完成させたいと思っているが，間に合うかどうかは微妙なところかもしれない。以下，探索結果を列挙する。
商用はたくさんあるし，Fortranで書くならnetlibでいいのだが，Cで書かれたライセンスフリーのソース公開のものはほとんど見つからなかったので，馬鹿馬鹿しいが，車輪の再発明みたいなことをしようと思う。といっても，アルゴリズムはより単純なものということで滑降シンプレックス法を使おうと思っているので，やるのはインプリメントだけだから，たいしたことはないが。
関係ないが，今夜はジュラシックパークの続編，ロストワールドのノーカット版が放映されるらしいので，録画しておくといいんじゃないかなと思っている。放映される時間帯には新幹線に乗っている確率が一番高いので，電話で録画を頼まねばならないかもしれない。ジュラシックパークが放映されたときは子どもたちが熱中してみていたので，録画しておけば子守りのときに役に立ちそうな気がするし，きっとやってくれるだろう。さらに関係ないが，昨夜はテレビ東京系でコンタクトが放送されたらしいのを見逃したのだが，先日上野のCDショップでDVDが2000円だったので衝動買いしてしまっており，別に悔しくはない。
夜までかかって最適化以外の部分を１つの関数に書き直すことはできたのだが，最適化のインプリメントができないままに終わった。新年会にも出ずにやったのにケリがつかなかったのは，今ひとつ気分が良くない。長野電鉄に乗るため，今日も21:30の直前にホームに着くように大学を出た。今日はダイヤ通りだった。藤原正彦「心は孤独な数学者」（新潮文庫）を読み進め，読了。本書はニュートン，ハミルトン，ラマヌジャンという３人の大数学者を生み出した風土を数学者藤原正彦が訪ねる旅である。優れた評伝となっているばかりでなく，紀行文としてもすばらしい。