X-Git-Url: http://gitweb.michael.orlitzky.com/?p=octave.git;a=blobdiff_plain;f=divided_difference.m;h=4705f0a3d8d1ff3cf699cca2dd21c04d30a3c2fa;hp=a9039b10c042caa773c4b162f0d6ee94e95734fe;hb=84b8fb9002d091f84d0205e923c3989d0138ec9e;hpb=c2dbf1f8827b39b5f5916675f7d319ec1e728e96 diff --git a/divided_difference.m b/divided_difference.m index a9039b1..4705f0a 100644 --- a/divided_difference.m +++ b/divided_difference.m @@ -15,18 +15,19 @@ function dd = divided_difference(f, xs) ## 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