Add the Octave c_inner_product function.

[octave.git] / c_inner_product.m

function ip = c_inner_product(w, a, b, v1, v2)
  integrand = @(x) w(x)*v1(x)*v2(x)
  ip = quad(integrand, a, b)
end