summaryrefslogtreecommitdiff
path: root/examples/plot.hs
blob: e4025fd1019c0d9341ea9ed045bab5233a5cc86f (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
import Numeric.LinearAlgebra
import Graphics.Plot
import Numeric.GSL(erf_Z, erf)

sombrero n = f x y where 
    (x,y) = meshdom range range
    range = linspace n (-2,2)
    f x y = exp (-r2) * cos (2*r2) where 
        r2 = x*x+y*y

f x =  sin x + 0.5 * sin (5*x)

gaussianPDF = erf_Z
cumdist x = 0.5 * (1+ erf (x/sqrt 2))

main = do
    let x = linspace 1000 (-4,4)
    mplot [f x]
    mplot [x, liftVector cumdist x,  liftVector gaussianPDF x]
    mesh (sombrero 40)