wvogel日記

自分用の技術備忘録が多めです.

2011-09-20から1日間の記事一覧

線形合同法

線形合同法による疑似乱数を作って遊んでたら面白いものができた。 linearRand :: Integral a => a -> a linearRand x = (456*x+759) `mod` 24 何を渡しても 15 が返ってくる笑 定数w 値が1違うだけで全然違いますね rand_list :: Integral a => [a] rand_l…

Monte Carlo

(?) :: (Ord a,Num a) => (a,a) -> Bool (?) (x,y) = 1 > x^2 + y^2 pi' :: Fractional a => a pi' = 4.0*fromIntegral (length $ filter (?) dots) / fromIntegral (length dots) dots :: (Enum a,Fractional a) => [(a,a)] dots = [(x,y) | x <- [0.01,0.0…