## OUTPUTS:
##
## * ``dd`` - The divided difference f[xs(1), xs(2),...]
- ##
+ ##
+
order = length(xs) - 1;
-
+
if (order < 0)
## Can't do anything here. Return nothing.
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