c_inner_product.m

   1 function ip = c_inner_product(w, a, b, v1, v2)
   2   ##
   3   ## The usual inner product defined on the space of continuous
   4   ## functions over the interval [a,b].
   5   ##
   6   ## INPUT:
   7   ##
   8   ##   * ``w`` -- The weight function.
   9   ##
  10   ##   * ``a`` -- The left endpoint of the interval.
  11   ##
  12   ##   * ``b`` -- The right endpoint of the interval.
  13   ##
  14   ##   * ``v1`` -- The first vector.
  15   ##
  16   ##   * ``v2`` -- The second vector.
  17   ##
  18   ## OUTPUT:
  19   ##
  20   ## The inner product <v1, v2>.
  21   ##
  22   integrand = @(x) w(x)*v1(x)*v2(x);
  23   ip = quad(integrand, a, b);
  24 end