Add Newton's method and some tests.

[octave.git] / divided_difference.m

diff --git a/divided_difference.m b/divided_difference.m
index a9039b10c042caa773c4b162f0d6ee94e95734fe..cd4b0181b10ed4eb7c460798a9a55c711ee813ab 100644 (file)
--- a/divided_difference.m
+++ b/divided_difference.m
@@ -15,7 +15,11 @@ function dd = divided_difference(f, xs)
   ## OUTPUTS:
   ##
   ##   * ``dd`` - The divided difference f[xs(1), xs(2),...]
-  ## 
+  ##
+  if (exist('../homework1/src', 'dir'))
+    addpath('../homework1/src');
+  end
+
   order = length(xs) - 1;
   
   if (order < 0)
@@ -23,10 +27,10 @@ function dd = divided_difference(f, xs)
     dd = NA;
   elseif (order == 0)
     ## Our base case.
-    dd = f(xs(1))
+    dd = f(xs(1));
   else
     ## Order >= 1.
-    cs = divided_difference_coefficients(xs)
+    cs = divided_difference_coefficients(xs);
     dd = dot(cs, f(xs));
   end
 end