■
7月1日の明け方近く。窓外にヨタカの声が聞こえる。今年はホトトギスとヨタカの声が良く聞こえる。気をつけて聞いているせいかもしれない。室内では子供の寝息が聞こえる。
Meteorological Analyses over the Northern Hemisphere の 12Z Tue 30 JUN 2009 の 200mb Streamlines and Isotaches のデータを見ると日本の南海上からアラビア半島にいたる西向きのジェット気流が形成されている。こんなこともあるんだ。へぇ。
sinnx (nは奇数)の不定積分を計算するために sin2n+1x=Σaksin( (2k+1)x ) を計算する。漸化式 sin2n+1x=(1/2)(1-cos(2x))sin2n-1x を下から逐次解いていく。ここまできれいに二項係数で書けるんだ。へぇ。
誤差伝播の計算で有効数字の桁数の概算をテイラー展開を使って実行する方法について覚書をまとめたいとずっと考えていてやっていない。