2006-09-01から1ヶ月間の記事一覧

p + k * p * (1 - p)  蛇足1

Plot[population[k, 51], {k, 1.5, 3.0}]

p + k * p * (1 - p)  補足

Mathematica 5.1の「ヘルプ」→「デモ」→「サウンドギャラリー」→「ロジスティックマップ」のコードを加工。 population[k_Real, n_Integer] := Module[{f, x}, f = Compile[{x}, Evaluate[N[x + k x (1 - x)]]]]; Nest[f, 0.3, n]] Plot[{population[k,…

p + k * p * (1 - p)

その後、OCamlを囓っている。数行の小物は組めるようになった。 let rec population n k = match n with 0 -> 0.3 | n -> let p = population (n - 1) k in p +. k *. p *. (1. -. p);; let rec range n m = if n > m then [] else n :: (range (n + 1) m);;…

Mathlink on IronPython 1.0

先日公開されたIronPython 1.0は、予想以上に優れモノで、.Netランタイム用に作られたモジュール(マネージコード)と全くシームレスに連携できる(ようだ)。対するオリジナルPython(+ctypesパッケージ)による既存のネイティブコードとの連携では、相互…