diff options
Diffstat (limited to 'examples/integrate.hs')
-rw-r--r-- | examples/integrate.hs | 24 |
1 files changed, 0 insertions, 24 deletions
diff --git a/examples/integrate.hs b/examples/integrate.hs deleted file mode 100644 index 3a03da6..0000000 --- a/examples/integrate.hs +++ /dev/null | |||
@@ -1,24 +0,0 @@ | |||
1 | -- Numerical integration | ||
2 | import Numeric.GSL | ||
3 | |||
4 | quad f a b = fst $ integrateQAGS 1E-9 100 f a b | ||
5 | |||
6 | -- A multiple integral can be easily defined using partial application | ||
7 | quad2 f y1 y2 g1 g2 = quad h y1 y2 | ||
8 | where | ||
9 | h y = quad (flip f y) (g1 y) (g2 y) | ||
10 | |||
11 | volSphere r = 8 * quad2 (\x y -> sqrt (r*r-x*x-y*y)) | ||
12 | 0 r (const 0) (\x->sqrt (r*r-x*x)) | ||
13 | |||
14 | -- wikipedia example | ||
15 | exw = quad2 f 7 10 (const 11) (const 14) | ||
16 | where | ||
17 | f x y = x**2 + 4*y | ||
18 | |||
19 | main = do | ||
20 | print $ quad (\x -> 4/(x^2+1)) 0 1 | ||
21 | print pi | ||
22 | print $ volSphere 2.5 | ||
23 | print $ 4/3*pi*2.5**3 | ||
24 | print $ exw | ||