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