function ip = c_inner_product(w, a, b, v1, v2) % % The usual inner product defined on the space of continuous % functions over the interval [a,b]. % % INPUT: % % * ``w`` -- The weight function. % % * ``a`` -- The left endpoint of the interval. % % * ``b`` -- The right endpoint of the interval. % % * ``v1`` -- The first vector. % % * ``v2`` -- The second vector. % % OUTPUT: % % The inner product . % integrand = @(x) w(x)*v1(x)*v2(x); ip = quad(integrand, a, b); end