diff options
Diffstat (limited to 'packages/hmatrix/examples/root.hs')
-rw-r--r-- | packages/hmatrix/examples/root.hs | 31 |
1 files changed, 0 insertions, 31 deletions
diff --git a/packages/hmatrix/examples/root.hs b/packages/hmatrix/examples/root.hs deleted file mode 100644 index 8546ff5..0000000 --- a/packages/hmatrix/examples/root.hs +++ /dev/null | |||
@@ -1,31 +0,0 @@ | |||
1 | -- root finding examples | ||
2 | import Numeric.GSL | ||
3 | import Numeric.LinearAlgebra | ||
4 | import Text.Printf(printf) | ||
5 | |||
6 | rosenbrock a b [x,y] = [ a*(1-x), b*(y-x^2) ] | ||
7 | |||
8 | test method = do | ||
9 | print method | ||
10 | let (s,p) = root method 1E-7 30 (rosenbrock 1 10) [-10,-5] | ||
11 | print s -- solution | ||
12 | disp p -- evolution of the algorithm | ||
13 | |||
14 | jacobian a b [x,y] = [ [-a , 0] | ||
15 | , [-2*b*x, b] ] | ||
16 | |||
17 | testJ method = do | ||
18 | print method | ||
19 | let (s,p) = rootJ method 1E-7 30 (rosenbrock 1 10) (jacobian 1 10) [-10,-5] | ||
20 | print s | ||
21 | disp p | ||
22 | |||
23 | disp = putStrLn . format " " (printf "%.3f") | ||
24 | |||
25 | main = do | ||
26 | test Hybrids | ||
27 | test Hybrid | ||
28 | test DNewton | ||
29 | test Broyden | ||
30 | |||
31 | mapM_ testJ [HybridsJ .. GNewton] | ||