stochastic-0.1.1.0: test/Normal3.hs

import Data.Stochastic
import Data.Stochastic.Types
import Data.Stochastic.Internal

import System.Random

normalInit :: StochProcess
normalInit = normalProcess 10 3

f :: Double -> StochProcess
f d = if d < 7 then normalProcess (0.9 * d) (0.3 * d)
  else if d > 13 then normalProcess (1.1 * d) (0.3 * d)
  else normalProcess d (0.3 * d)

normal3 :: StochProcess
normal3 = composeProcess 3 normalInit f

main = do
    gen <- newStdGen
    print "normal3"
    print $ runProcessN 100000 normal3 gen