root finding with jacobian

1 files changed, 14 insertions, 0 deletions

diff --git a/examples/root.hs b/examples/root.hs
index 69db243..2a24f0f 100644
--- a/examples/root.hs
+++ b/examples/root.hs
@@ -11,6 +11,15 @@ test method = do
     print s -- solution
     disp p -- evolution of the algorithm
 
+jacobian a b [x,y] = [ [-a    , 0]
+                     , [-2*b*x, b] ]
+
+testJ method = do
+    print method
+    let (s,p) = rootJ method 1E-7 30 (rosenbrock 1 10) (jacobian 1 10) [-10,-5]
+    print s
+    disp p
+
 disp = putStrLn . format "  " (printf "%.3f")
 
 main = do
@@ -18,3 +27,8 @@ main = do
     test Hybrid
     test DNewton
     test Broyden
+
+    testJ HybridsJ
+    testJ HybridJ
+    testJ Newton
+    testJ GNewton