X-Git-Url: http://gitweb.michael.orlitzky.com/?p=octave.git;a=blobdiff_plain;f=divided_difference_coefficients.m;h=a40f5597d6b47712ad7e0b397585d71641b80ec8;hp=53d5799e922af57d70c74604cca5814d4d4b7b4d;hb=6cd387b2e42c3273c6cd48ad2e7e551bddd86474;hpb=c2dbf1f8827b39b5f5916675f7d319ec1e728e96 diff --git a/divided_difference_coefficients.m b/divided_difference_coefficients.m index 53d5799..a40f559 100644 --- a/divided_difference_coefficients.m +++ b/divided_difference_coefficients.m @@ -8,16 +8,15 @@ function coefficients = divided_difference_coefficients(xs) ## OUTPUTS: ## ## * ``coefficients`` - The vector of coefficients such that - ## dot(coefficients, f(xs)) == dd. Used to solve linear systems. + ## dot(coefficients, f(xs)) == f[xs]. Used to solve linear systems. ## coefficients = []; - + for xj = xs this_coeff = 1; for xi = xs if (xi != xj) - ## Append (xj - xi) to the vector of coefficients. this_coeff = this_coeff * (1 / (xj - xi)); end end