4章の後半くらいから理解があやしくなった。

カーネル多変量解析―非線形データ解析の新しい展開 (シリーズ確率と情報の科学)
赤穂 昭太郎
岩波書店 / ￥ 3,675 ()
在庫あり。

何度か読まないといけないな。

なんというか燗つけるにはまだちょっと暑いし、かといって華やかな酒が吞みたいわけではないしなぁと、久々の焼酎、しかも球磨。

たまにはいいかな。

娘とお菓子を作ってきた。

• さつまいものごまビスケット
• ヨーグルトゼリー

黙々と蕎麦打ちすんのも好きだけど、子供とお菓子を作るのも結構楽しい。

ジャガイモが伸びてきた。

島おくらも収穫開始。

NHK 趣味の園芸ビギナーズ & やさいの時間 2009年 09月号 [雑誌]

日本放送出版協会 / ￥ 650 ()
在庫あり。

まわりは白菜植えたりキャベツ植えたりしているのに、自分の畑は空きがなくて出遅れてる感があって焦る。

モロヘイヤとつるなしインゲンのレーンをどけて冬野菜に備えた。 そういえばこの前植えたつるなしインゲンは固くて食えたもんじゃなかった。(実が大きくならずに枯れはじめたし)

大根はちょっと無理そうなのでほうれん草と小松菜を植えようと思っている。

オプションで指定する。

sudo port install swig +ocaml

でも、なんかエラーが出てて、いまのとこうまいこと使えてない。

floatしかべき乗の演算子がないのは何故？

# 2.0 ** 3.0;;
- : float = 8.

int

# 2 ** 3;;
Characters 0-1:
  2 ** 3;;
  ^
Error: This expression has type int but is here used with type float

let rec power m = function
    0 -> 1
  | n -> m * power m (n-1)
;;

let rec factorial = function
    0 -> 1
  | n -> n * factorial (n-1)
;;

let rec calc_pi = function
    0 -> 3.0
  | n -> (float_of_int (3 * factorial(2*n))) /.
 (float_of_int ((power 2 (4*n)) * (power (factorial n) 2) * (2*n+1))) +. calc_pi (n-1)
;;

これだとcalc_pi 6ぐらいでneg_infinityになってしまう。あと型で結構つまづく。

普通に求めるヤツ

let mc () =
let x = Random.float 1.0 in
let y = Random.float 1.0 in
if sqrt ((x ** 2.0) +. (y ** 2.0)) <= 1.0 then 1
else 0
;;

let rec count n res = match n with
  0 -> res
| _ -> count (n-1) (res + mc())
;;

let calc_pi n =
  let res = count n 0 in
  (float_of_int res)*.4.0 /. (float_of_int n)
;;

vector loversは良いが、アルバム持っているからな。octogenのもたついてるっぽいリズムがメローでよかった。mirror musicってsilicon soulの別ユニットなのか。

SOMA COMA PART.3
オムニバス,オクトゲン,ベクター・ラバーズ,スラム,ミラー・ミュージック,パーシー・X vs ブラッドシュガー,ザ・ブラック・ドッグ,デヴィッド・ホルムス&スチュアート・マクミラン,アレックス・スモーク,スキントレード
SOMA/OCTAVE-LAB / ￥ 2,678 (2009-08-12)
在庫あり。

SOMAが毎年夏にリリースするエレクトリック・チルアウト/ダウンテンポコンピレーション第3弾が登場 !! 野外レイヴシーズンにリリースされヨーロッパでも人気の高いクラブミュージックファンにフィットする気持ちイイダウンテンポのサマー・テクノチューンを収録。グラスゴーテクノシーンの至宝、アレックス・スモーク、SOMAオーナーのスラムを始め、ブラック・ドッグからアンドリュー・ウェザオールの変名プロジェクト、ブラッドシュガーの90年代の名曲、デヴィッド・ホルムズまでUKテクノシーンのメロウサイドを切り取った注目作 !

ビュフォンの針でとおもったんだけど、なぜかarcsinのテーラー展開が目に留まったので。

def factorial(i):
   return 1 if i ==0 else i * factorial(i-1)

def calc_pi(i):
   pi = 0
   for n in range(i+1):
       pi += 3.0*factorial(2*n)/(2**(4*n)*((factorial(n))**2)*(2*n+1))
       print pi
   return pi


if __name__ == "__main__":
   print calc_pi(16)

「無限と連続」の数学―微分積分学の基礎理論案内
瀬山 士郎
東京図書 / ￥ 2,625 ()
在庫あり。

itunesいじっててmetallicaにはまった。

nothing else matters と and justice for allが好き。metallicaといえば大体この二枚。

Metallica
Metallica
Universal/Mercury / ￥ 1,940 (1991-08-12)
在庫あり。

メタル・ジャスティス
メタリカ
UNIVERSAL INTERNATIONAL(P)(M) / ￥ 2,500 (2008-09-12)
通常10～14日以内に発送

そういえば、初めて買ったアルバムがAppetite for Destructionなのだよね。

Appetite for Destruction
Guns N' Roses
Geffen / ￥ 1,429 (1990-10-25)
在庫あり。

patienceも良いですね。

G N' R Lies
Guns N' Roses
Universal/Polygram / ￥ 1,020 (1990-10-25)
在庫あり。