blob: f950aa59770016e11b834e5ebbabdbf5a6a03d7d (
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.Special(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, mapVector cumdist x, mapVector gaussianPDF x]
mesh (sombrero 40)
|
