「今はもうない」を読んだ。

久々に森作品を読んだ。これは全巻読んでいて、人間関係を全部覚えている人がいたら最初のほうでわかるお話ですね。私はそこまで読んでいないので。。。ネタバレとしては、主人公は誰でしょう、ということですかね（トリックのネタバレではないが）。

長池の伝説＠長池駅

大蛇の伝説があるそうだけど、この文面だけみるとそんなに大蛇が悪そうに見えないんだけれど。。。埋め立てたほうがたたられそう。。。

京都街道と奈良街道の分岐点

アルプラザ城陽の近くで撮影。あまりにひっそりしていて気づかれない。。。

文化パルク城陽の前を通りかかった。

変わった形の建物があるなあ、と近づいていくと。。。

文化パルク城陽でした。

http://www.city.joyo.kyoto.jp/parc/parc-indexnew.htm

丸いダイエーのパビリオン＠ポートピアのような部分はプラネタリウムらしいです。

残念ながら、この日は休館日でした。

カロリ。の新CMソングは懐かしいVillage peopleのIn the navyというよりピンクレディーのピンク・タイフーンか。。。

佐々木希ちゃんがかわいらしいカロリ。の新CM。

http://www.suntory.co.jp/rtd/calori/cm/index.html

かかっている曲はヴィレッジ・ピープルのイン・ザ・ネイビーですよ。

でもこれって海軍に入ろう！という曲なんですが、実際はゲイの人（ヴィレッジ・ピープルはそういう売り）が海軍に入って男たちのなかでうほうほとかそういうことでしょ。。。

で、それをピンクレディーがカバーして、ピンク・タイフーンにした。

しかし歌詞のインザネイビーがピンクレイディーになったのには度肝を抜かれた。

CMではさらにシンカローリーになってますが、どっちかというとピンクレイディーに近いなあ。

やっとプラズマテレビを買った（Panasonic Viera 42型）

もういい加減にテレビを買い換えようと思い立って、パナソニックのビエラ42型を買ってきた。プラズマで、品番はTH-P42R1。

いままでこのテレビデオだったんで、格段の違いだ。

ハードディスク200GBｙte内蔵で45時間くらい録画できるし、その録画も番組表をボタンひとつで呼び出して録画設定するだけ。超簡単。マニュアルにはかかれてないけど、ダブルチューナーなんで録画と視聴は別々にできますよ。

さっそくBSもWOWOWも登録した。今まで何でデジタルテレビ買わなかったんだろうと後悔。

さて、テレビも新しくなったし、ゲーム機も買い換えるか（続く）。

「姑獲鳥の夏」を読んだ。

だいぶ前の作品だけれど、最近の脳ブーム？からみても古くないですね。

語り手が見てない、というのは反則っぽいけど。。。でも民俗とかの話が面白いからいいです。

映画にもなったようですが、見てない。。。映像化はかなり難しかったでしょうが、また見てみよう。

SH-03Aが故障！→修理に出してやっと帰ってきた。

一年を待たずして壊れたよ。症状は、メールの送信ができなくなるというもの。受信はできるのだが、送信BOXを開こうとすると待ち受け画面に戻るし、送信しようとしても「送信に失敗しました」エラーがでる。全部メモリー消そうとしてもエラー。ドコモショップに持っていったら全設定消去をされたのだが、それでも動かず。もうどうしょうもない。

で、修理にだしたら、ガワは全部変わって、かつ中身の基板も変わって、

。。。てもう絶対新しい端末おくってきただけやんーーーー。

こういう症状の人って検索しても1件しかなかった。皆さん大丈夫？？？

平成22年センター試験 数学II・数学B 第6問をkeisan.casio.jpにUP!

毎年恒例のBASIC移植シリーズ。問題は、

一般に、三つの正の数について、どの二つの数の和も残りの数より大きければ、それらを三辺の長さとする三角形が存在する。逆に、すべての三角形において、どの二辺の長さの和も残りの一辺の長さより大きい。
この事実を用いて、自然数Nを三角形の三辺の長さとなり得る三つの自然数a,b,c (a≦b≦c)の和として表す方法をすべて列挙し、その総数を求める。

でリンクはこれ↓

平成22年センター試験　数学II・数学B　第6問

プログラムは死ぬほど簡単で、

x=1;
for (a = 1; a <= int(n/3.); a=a+1) {
    for (b = a ; b <= int((n-a)/2) ; b = b+1) {
        c = n-a-b;
        if (c < (a+b)) {
     println (x, a, b, c);
        x = x + 1;
        }
    }
}

でした。

「アバター」（3D版）をようやく見てきた。

想像していたよりずっとちゃんと論理的なお話でびっくりした。以下ネタバレあり。

木々が根で情報ネットワークを組んでいる、いわばバイオコンピュータが星全体のようなところが舞台のお話（「銀河ヒッチコックガイド」の地球みたいなものか）。そういう星で生物が進化していくなら、それにアクセスする手段としての器官ができるのが自然。なので、ナヴィをはじめ、すべての生物には触手のような神経ネットワークをつなぎ合わせる器官がある。それを木々にではなく、同一/他種目間のコミュニケーションに使えるようになっていくことも自然ですね（全部同一インターフェースのはずだし）。で、最後のシーンも何も考えないと超自然的？な感じに見えるかもしれないが、人間の脳の情報をバイオコンピュータを介してすべてナヴィに送り込む、つまりグレッグ・イーガンの小説でよくでてくる脳の情報を全部コンピュータに移したようなものですね。移した結果、基の個体は死ななくてはならないというところも似てるな。

。。。とはいうものの全体的な印象はすごくナウシカっぽいというか。触手うねうねと森の意志ってのががそう見えるのか。

あとはやっぱり神経接続とガラス張りの部屋で暴走を見てるシーンは死ぬほどエヴァンゲリオンっぽいし、実写版ダグラムは動いてるし、どう考えても日本のアニメの影響はあると思うけどなあ。

あー、3Dはやっぱり面白かったですよ。でも一番驚いたのは本編の前のメガネの説明のところで。。。あとはなんか飛び出す絵本っぽい。立体というより奥行きをつけて平面を並べているような感じ。でもガス弾で思わず体をよけたのは内緒だ。

しかし飛行機のパイロットの女性はすごく私のタイプなんですが（肉感的で）、残念なことになってしまって。。。

遅延微分方程式（酵素反応論）をルンゲクッタ8次でExcel VBAで。

次は酵素反応論からの例。方程式は

y1'(x) = I - z*y1(x)

y2'(x) = z*y1(x) - y2(x)

y3'(x) = y2(x) - y3(x)

y4'(x) = y3(x) - 0.5*y4(x)

ここで

z=1/(1+0.0005*y4(x-4)^3)

で遅延している。

計算結果はこちら。

 

この文献の結果を再現してますね。だいぶ慣れてきた。次いってみよう。

遅延カーマック・マッケンドリック方程式をルンゲクッタ8次でExcel VBAで。

だいぶ前に伝染病の方程式をやってみたけど、その遅延方程式ヴァージョンを考える。

免疫を持つ人が一定期間τ1(=10)で再び感染してしまい、かつ潜伏期間τ2(=1)がある場合のモデルは

y1'(x) = -y1(x)y2(x-1) + y2(x-10)

y2'(x) = y1(x)y2(x-1) - y2(x)

y3'(x) = y2(x) - y2(x-10)

となる。ここでy1は感染可能な人、y2は感染者、y3は治癒した人を表す。計算結果はこちら。

なるほど、周期的突発があるんだな。

参考文献はこちら（続く）。

遅延微分方程式をルンゲクッタ8次でExcel VBAで（Mackey-Glass方程式編）

遅延微分方程式を高次ルンゲクッタで解くのに、以前は相当複雑な方法でやってみた。この常微分方程式の数値解法にもっと簡単な方法があったのでやってみる。

方法は、過去のxの変わりにgをつかうというもの。

 

あー、前は過去のxをわざわざ補完してつかってたのだ。これのほうが圧倒的に簡単。

でマッキーグラス方程式を計算してみた。

計算結果はこれ。前は非常に複雑だったのがすこぶる簡単に計算できた。

これからしばらく遅延微分方程式を解いていってみよう（続く）。

「犬と主人問題」をExcel VBA (Dormand&Prince8次）で計算。

主人がｙ軸方向に速度vで走っていて、犬がwの速度で主人に向かって走っていくときの犬の軌跡を考える。

微分方程式にすると、

xy'' =　√(1+(y')^2)

となる。実は厳密解が変数分離で解けるけど、ここはルンゲクッタ8次で計算してみよう。

主人が遅いとほぼまっすぐ追いつくけど、早いとだいぶ犬はがんばらないといけないな。おもしろいな。

Billy Bat 2巻をかった。

下山事件、白州次郎とか面白そうなねたで1巻が始まったBilly Batだけれど。。。

キリストとユダの話とか、だんだん風呂敷が大きくなってきた。大丈夫かなあ。面白いというより心配になってきた。

1000均ショップって。。。

目を疑った。一桁違う。。。向島付近のスーパー内にて。

銀河力学をルンゲクッタ8次でExcelで計算。

「常微分方程式の数値解法I」に出ていた例。ハミルトニアンを

H = (p1^2 + p2^2 +p3^2)/2 + Ω*(p1*q2 - p2*q1) + A*ln(C + q1^2/a^2 + q2^2/b^2+q3^2/c^2)

でH=2としてDormand-Prince8次で計算してみた。

三次元的に書いたのがこれ。

q2=0のポアンカレ写像がこれ。h=1/10でもっともらしい値がでてる。

空間微分を含むFitzHugh-南雲方程式をExcelで高次ルンゲクッタで計算

だいぶ前に空間微分のない常微分方程式を計算してみたけれど、今度は

τu_t = u (u-α) (1-u) - v + ε²u_xx

v_t = u - γv

を計算してみよう。

空間微分はいつものように2次精度の差分で近似して、得られた常微分方程式系を5次のルンゲクッタを使って計算する。

α=0,τ=0.018,γ=1,ε=0.05

です。

ちゃんとパルスが伝播していく様子がよくわかる。

「ファウンデーション」（アイザック・アシモフ）を読んだ。

いつか読もうと思っていて、シリーズが長いのでちょっとしり込みしていたのを意を決して読んでみました。確かにすごく面白い。

心理歴史学によって未来を予想したハリ・セルダンが要所要所の危機で現れるさまや、危機をその時々の人物が暴力を使わずに解決していくさまは本当に面白い。でもまだ物語は始まったばかりでこの巻は終わっちゃう。次も楽しみ。

ただ、宇宙旅行（というかワープのような航法をもっている）ができるのに、原子力の技術がないというのはいまいち納得しかねるなあ。そこだけはちょっと不満。だけどそんなことを吹き飛ばすくらいのパワーがあります。

あー、あとガンダム00でもこうやって死後もホログラムで出てくるという博士がでてくるとか。パクリか？（いや、ほとんど見たことないので全然的外れかもしれませんが）

「摩天楼の怪人」を読んだ。

島田荘司さんの摩天楼の怪人を読みました。創元推理文庫って初めて読んだ気がする。

ネタバレかも知れませんが、高所恐怖症の私は最後のほうでくらくらした。私なら絶対落ちてる。年代がポイントとずっと言っているのでそういうお話ですが、時計台の針はルパン三世のカリオストロの城のようだ。。。

しかし（最大のネタバレ）今ならgoogle earthがあるので一発でばれますね。google恐るべし。

ExcelでCML (Coupled Map Lattice) を計算。

こういうのはExcelにぴったりの計算だと思いますが。

x^k+1_n = (1-g) * f(x^k_n) + g (f(x^k_n+1) + f(x^k_n-1) ）/2

でf(x) = 1-a*x*x

とする。いろんなa,gのパラメータで書いたのがこちら。

しかし最初、この

http://en.wikipedia.org/wiki/Coupled_map_lattice

で書かれたようなパターンにぜんぜんならなくて、なぜかなあと思っていたら2つずつ時間ステップ飛ばして図示してるんだって。。。最初から言ってくれよ。

この本も会社の図書室にありました。

が、こういうのは昔ちょっと流行って驚くほど急激にすたれたな。一発屋？

「はじめての現代数学」を読んだ。

ハヤカワ・ノンフィクション文庫で出ていて気になっていたのを読んで見た。

もともとは1988年出版ですが、2009年に文庫化されたということで、ちゃんとポアンカレ予想とかフェルマー予想が解決されたことまで付け加えられています。

専門家向けではなくトポロジーや論理学の話がわかりやすく書かれていて、こういう本ってなかなかないのでお勧めですよ。

しかしこの当時はカタストロフィー理論が流行っていたようなんですが、今はこれはどうなってるのかな。私が大学生の時代ですでに全然聞かなかったし（本は持ってますが）。

非局所結合TDGL方程式(1次元)を高次精度ルンゲクッタでExcelで。

空間微分が非局所結合（積分）に変わったTDGL方程式を見てみよう。式はこちら。

で計算結果。C0=0,C1=-2,C2=2,K=1.1で計算。

うーむ、あってるのかどうかわからん。参考文献があんまりないからなあ。私が見たことあるのは、

「数理科学1997年6月号　結合振動子系のダイナミクス」と、このリズム現象の数理くらい。

TDGL方程式(1次元)を高次精度ルンゲクッタでExcelで。

TDGL方程式、

の空間微分を2次で差分化し、常微分方程式化して5次精度のルンゲクッタ（Cash-Karp)で計算してみた。

パラメータはC0=-4/3,C1=2,C2=-1。

あーこれはホール解なのか？結構面白いパターンが書けた。でもう少し拡張して非局所版を計算してみよう（続く）。

福知山のかき末さんに今年も行って来た。

毎年恒例です。かきをたらふく食べるために、福知山にあるかき末さんにいってきました。

4500円の定食を選択。上定食はこれにお刺身がつくだけなんでパス。

つきだしと。。。

かき鍋です。4人前でこんなにある。

さらにかきフライも熱々。

最後にかきご飯がつくけど、もうおなかいっぱいで全部は食べきれず、折り詰めにしてもらった。これはほんとうにおなかいっぱいになりますよ。大満足。

最近の容疑者はみんな同じパーカーを着てる。

市橋容疑者で有名になった紫っぽいパーカーだけど、元光ゲンジの赤坂容疑者も同じのを着ていてあれ？と思った。でまた女性を殺した容疑のやつが同じやつを着ていてあれれ？と。こんなやつ↓（絵心0なので。。。）警察が準備してるんだろうね。標準品なんだろうか。

(1/31追記）

さっき「バンキシャ！」を見てたらこの話題やっていた。千葉県警の特注品の被疑者護送用パーカーだって。なんと袖がなく、本当に護送専用。顔の隠し方も調整可能とのこと。9年前から使ってるそうですが、最近千葉で凶悪事件が多いから目立つそうです。

京都グランヴィアの五山望に行ってきた。

高級鉄板焼です。コックさんが目の前にいるので写真撮りづらく、玄関のみ。

3万円のコースなり。

フォアグラ/海老/蟹/あわび/サーロイン＆フィレ/いくら・蟹釜飯とかもう限界まで高級食材をおなかいっぱい。しかしもう”ごちになります”より設定金額高くなってないか？

5人で20万円ほどでした。

でも高いだけのことはありましたよ。

田中神社（伏見）に行ってきた。

馬の神様らしいです。近くに京都競馬場もあるしね。

虫の神様とも。

立派な馬の像がわきにあったけど、一部が立派過ぎて撮影は自粛。。。

大手筋商店街のSOLAR MOON

伏見の大手筋商店街には太陽光発電がついてますね。

でもソーラームーン？セーラームーンなら知ってるけど。。。

http://www.otesuji.jp/solar.html

「御手洗潔のメロディ」を読んだ。

もはやだんだんミステリーものでもなくなって御手洗さんスーパースター小説のようになってきてるなあ。これはこれでファンが多いのはわからんでもないですが。。。

マイルス・デイビスまで登場しちゃってるよ。脳科学でノーベル賞までとりそうだし。もはや初期のおどろおどろした雰囲気は皆無ですが、逆に女性ファンがふえていっていそう。やおい小説？

カドムツェフ・ペトリアシビリ方程式をExcelで計算

2次元のkdV方程式、カドムツェフ・ペトビアシビリ（Kadomtsev-Petviashivili)方程式

の2ソリトン解はよくこうやって書かれている。

が、これ微分計算するのがめちゃくちゃめんどくさい。そこで、Maximaで計算させると次のように長い式が得られる。

これをそのままExcelのVBAに持ってきて図示したのがこれ。

なるほど、交差しながら動いていく。

殺人方程式（切断された死体の問題）を読んだ。

新年からこういうタイトルの本はどうかと思いますが、館シリーズ以外のものをはじめて読みました。

トリックはさすがに理学部物理学科卒業の私なので、最初からわからないと恥ずかしいということでわかりましたが、犯人はあー、そういうことかと。ただ館シリーズのおどろおどろしい感じが好きな私としては双子の兄弟はさわやかすぎかな。やっぱり次は4巻もある暗黒館の殺人を読むしかないのか。。。

あけましておめでとうございます。

さて、今年はブログで何をやろうかな。Excelで計算シリーズでは、昨年は8次のルンゲクッタやシンプレクティックで常微分方程式を解くシリーズからスペクトル法、2次元の偏微分方程式（反応拡散系）くらいまでいった。

そろそろねた切れ感があるけれど。。。。

やっぱり3次元まで上げるか。。。でもうちの非力なPCではちょっと計算するだけで数時間かかりそう。あとは確率微分方程式かなあ。それとルンゲクッタ以外の方法とか。

また、おなじみのあのCMに出てる人は誰？あの曲は何？シリーズとおっさんホイホイ（特撮＆アイドル）は継続。

グルメ記事は増やしたいなあ。経済力が続く限り。