X-Git-Url: http://gitweb.michael.orlitzky.com/?p=octave.git;a=blobdiff_plain;f=divided_difference_coefficients.m;h=92e991575b54460cbf3475ebb92e68d7faa8253f;hp=53d5799e922af57d70c74604cca5814d4d4b7b4d;hb=1a6f56b0dd6750649725b2fd07edb3fe0850a886;hpb=c2dbf1f8827b39b5f5916675f7d319ec1e728e96 diff --git a/divided_difference_coefficients.m b/divided_difference_coefficients.m index 53d5799..92e9915 100644 --- a/divided_difference_coefficients.m +++ b/divided_difference_coefficients.m @@ -1,23 +1,22 @@ function coefficients = divided_difference_coefficients(xs) - ## Compute divided difference coefficients of `f` at points `xs`. - ## - ## INPUTS: - ## - ## * ``xs`` - A vector containing x-coordinates. - ## - ## OUTPUTS: - ## - ## * ``coefficients`` - The vector of coefficients such that - ## dot(coefficients, f(xs)) == dd. Used to solve linear systems. - ## + % Compute divided difference coefficients of `f` at points `xs`. + % + % INPUTS: + % + % * ``xs`` - A vector containing x-coordinates. + % + % OUTPUTS: + % + % * ``coefficients`` - The vector of coefficients such that + % 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